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The Supersonic Conical Wing of Minimum 
Drag’ 


S. H. TSIENT 
Cornell Unwersity 


ABSTRACT 


A theoretical investigation is made of thin lifting wings 
of the conical family, with subsonic leading edges, cambered to 
give minimum pressure drag at supersonic speeds. The method 
is an approximate one based on the linear small-perturbation 
theory and an approximation to the lift distribution in series 
form. Two different assumptions are made concerning leading- 
edge suction: (1) full theoretical value of this thrust component, 
and (2) no thrust component, representing an extreme case of 
flow separation at the leading edge. 

The results show that in case (1) the minimum drag coefficient 
of the conical wing is almost identically the same as that of a flat 
triangular wing, at the same lift coefficient, except when the 
leading edge lies close behind the Mach cone of the apex. Thus 
little improvement over a flat wing can be achieved, within 
the conical family, if full leading-edge suction is obtained. In 
case (2), on the other hand, substantial improvement over the 
flat delta wing can be achieved by suitable camber; in fact, the 
drag can be reduced to a value that is again approximately the 
same as that of the flat wing with full leading-edge suction. 

The camber distributions and lift distributions of these wings 
are worked out and are presented in charts. 

In a brief supplementary investigation, the effect of chordwise 
camber of parabolic shape is also determined. It is found that 
for cases with leading-edge suction there is negligible improve- 
ment in minimum drag-due-to-lift as compared to the conical 
family. For cases with suction omitted, this type of camber 
affords appreciable improvement in drag, at least for wings 
whose leading edges are swept back well inside the Mach cone. 


NOTATION 


c = root chord 
Cp coefficient of drag due to lift 
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Cop, coefficient of thrust due to leading-edge singularities 
i.c., coefficient of total leading-edge suction force- 
D, 
defined by — 
” wing area X dynamic pressure 
lift coefficient 
drag due to lift 
thrust due to leading-edge singularities 
leading-edge suction force 
elliptical integral of the second kind 
analytic complex function (of ¢), whose real part is u 
analytic complex function (of 4), whose real part is u 
analytic complex function (of ¢), whose real part is w 
complete elliptical integral of the first kind. Also, 
constant of integration in the evaluation of wing 
camber 
free-stream Mach Number 
leading-edge slope in xy-plane 
modulus of complex variable ¢ 


rT 2 
the integral : sin™ @ \/ 1 — m? sin? 6 dé 
0 


modulus of complex variable « 

complex variable 

free-stream velocity 

nondimensional perturbed velocities in x, y, and 2 
directions, respectively (perturbed velocity /free- 
stream velocity ) 

V/ M?— 1 

complex variable 

y/x, also complex variable 

sin~! (¢/m) 

Vy? + 2?/x 


argument of complex variable « 


<4 Dd co“ a ®D 


argument of complex variable ¢ 


INTRODUCTION 


be» DETERMINATION OF THE OPTIMUM configuration 
for minimum pressure drag of a three-dimensional 
supersonic wing is in general a complex problem, even 
within the approximation of linear theory in which 
thickness drag and lifting drag are independent. 
Jones! has established a general criterion for the opti- 
mum wing loading, but a conclusive solution like the 
subsonic case based on lifting-line theory is still lacking. 
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MACH CONE 


Fic. 1. Coordinate system. 


However, the problem has been made accessible in 
certain investigations by (1) confining attention to 
slender wings,” * (2) confining attention to fixed plan 
forms with one or more parameters held constant,‘*~® 
and (3) by successive adjustments.’ Graham’s drag- 
reduction method’ requires fewer restrictions on the 
parameters (camber, twist, and thickness) and theo- 
retically leads to the exact solution for a wing with 
fixed plan form, but unlike the investigations referred 
to in (1) and (2), it is not one that permits a direct ap- 
proach. The present paper deals with optimum lifting 
conical wings of zero thickness symmetric about the 
streamwise axes. Attempts are made to find the 
theoretically obtainable minimum values of Cp/C,” 
and the corresponding shapes of such wings at various 
aspect ratios, based on linear theory. 

The “‘conical’’ limitation brings great simplifications 
so as to make possible a direct approach by available 
mathematical methods and to dispose of the limitation 
of slenderness. The accessibility is similar to those 
in category (2) mentioned above, but instead of vary- 
ing only one parameter (such as, for the nonlifting case, 
the position of the ridge line of a delta wing with con- 
stant‘ or linearly-varied® thickness-chord ratio), here 
the variables consist of camber and twist combining 
so as to give conical shapes, but otherwise arbitrary. 
The restriction of generality involved, therefore, is 
just the restriction to a function of one independent 
variable (the conical variable), whereas, for complete 
generality the shape would be described by a function 
of two independent variables. <A similar restriction 
was made in reference 6, where the rectangular wing 
was permitted to have arbitrary camber but no twist. 
Thus, this investigation also leads to supersonic wings 
of “restricted minimum” drag and indicates some of 
the possibilities of drag-reduction compared with the 
flat plate. 


DECEMBER, 1955 


Another restriction has been imposed—i.e., the lead 
ing edges are subsonic. The superiority of subsonic 
over supersonic leading edges for such wings has long 
been recognized.* The effect of leading-edge suction 
and the accompanying phenomenon of possible flow 
separation is thus involved, as the total leading-edge 
thrust is not of negligible order for most aspect ratios, 

The method of series representation for the velocity 
components used in this paper was first applied to 
conical flow in a paper by Baldwin’ for evaluating 
camber shapes and drags corresponding to certain 
specified lift distributions with no leading-edge singu- 
larities. 

For the sake of simplicity, the analysis is carried 
out at Mach Number v2. 
verted to other Mach Numbers by simple similarity 


Results are easily con- 
relations. 


ConrIcAL FLOW THEORY AND DRAG DUE To LIFT 


The wing investigated is sketched in Fig. 1, with the 
wing surface lying approximately in the xy-plane and 
the x-axis in the free-stream direction. The conical 
wing creates inside the Mach cone a conical flow field 
of which the velocity components can be represented 
by the two-dimensional Laplace equation after the 
Chaplygin transformation, 

e = se” = [p/(1 ++ V1 — pe" — i! , 
(1) 
with p= (%* yy? + 2°)/x and r = tan—'2/y | 


After another well-known transformation, viz., 
¢ = 2e/(1 + e?) (2) 


the region inside the unit circle of the Chaplygin plane 
(see Fig. 2) is mapped into the entire ¢-plane (see Fig. 
3) with the advantage that the wing coordinate is not 
distorted—i.e., when z ~ 0, ¢ ~ y/x. In these com- 
plex coordinates the existing conical theory" '! shows 
that the perturbation velocity components are related 











Fic. 2. Chaplygin plane 
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SUPERSONIC CONICAL 


by compatibility relations such that when w is given, 
one can evaluate w by the following integral: 


fi v1i-2,, 
g(C) = —1 : df(c) (3 


e 4 


where /(¢) = an analytic complex function whose 
real part is wu and g(¢ an analytic complex function 
whose real part 1S W. 

Also, within the linear theory, one can obtain the 
following expressions for coefficients of lift and drag 


due to this lift) in the same coordinates: 


| ein 
m/v 


(p= — [ u(C)w(C)deé — Cp (5) 
m Jo 
in which u(¢), w(¢) are values on the upper surface of 
the wing ( = +0, ¢ = y/x) and Cp, is the coefficient 
of total thrust due to the leading-edge singularities, 
commonly designated as the leading-edge suction force. 
The presence of leading-edge thrust has been recog- 
nized in thin-airfoil theory'? as a result of the square- 
root type singularity of pressure distribution for sub- 
sonic two-dimensional wings. The same concept can 
be extended to conical wings of subsonic leading edges 
in view of the fact that the pressure distribution in the 
neighborhood of the leading edge is identical in the 
limit with that for a two-dimensional wing in subsonic 
flow. Under this reasonable assumption, it has been 
found that for a flat delta wing of subsonic leading 
edges, | 
Cp, = [C.? (4am)| V1 — Bm? (6) 
where 6 = V M? — 1. As it will be shown later that 
the lift distribution of a conical wing may be repre- 
sented by a linear combination of components including 
one which has the flat-plate distribution, Eq. (6) can 
be extended to such wings, if the value of C, is inter- 
preted as the total lift contributed by that particular 
component. To avoid confusion, one may denote this 
portion of lift by C,’, to distinguish it from the total 
lift coefficient C,, and say, for conical wings 


Cy, = [CL’?2/(4am)] V1 — B*m? (CL’ < CL) (7) 


SERIES REPRESENTATION FOR LIFT DISTRIBUTION 


In the present approach, the minimum-drag problem 
will be solved by the determination of the optimum lift 
distribution, which will be represented in a form of a 
series. The constituent terms are linearly combined 
with arbitrary coefficients, the values of which will be 
The necessary conditions that the series 
(1) consistency with the boundary 


determined. 
should satisfy are: 
conditions (i.e., antisymmetry about the real axis, 
symmetry about the center line, and vanishing on the 
Mach cone and on the real axis off the wing), (2) con- 
vergence (so that a finite number of terms may be 


) 


taken for sufficient approximation), (3) correct singu- 
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WING MACH CONE 


-co- —| -m (e) +m +! ++eo 








Fic. 3. ¢-plane 


larities acceptable from physical considerations, and 
(4) no loss of generality. 

It is found that the following infinite series, which 
is an extended form of Baldwin's,’ satisfies all the above 


requirements: 


j : (+m : m 
u(¢) = R.P. ya 1 log ~ +a + 
de Ver — = 
$ > 9\2 1 
~ © is Vo? — m*\" { 
E bus ($= VE 8) 
n=1 m 


in which dp, a’, and »,_; are the arbitrary coefficients 
The first and second terms are respectively 
15 


(real). 
the ‘“‘constant lift’’'* and ‘‘constant upwash’ 
ponents, while the remaining part as a whole repre- 


com- 


sents components giving a lift distribution of arbitrary 
shape but with zero intensity at the leading edge. 
That the above series satisfies conditions (3) and (4) 
is obvious. Thus it need only be shown that the in 
finite series satisfies conditions (1) and (2). This is 
easily seen after a Joukowski transformation 


¢ = (m/2) [t+ (1/t)] (9) 
or 
t= (¢ + Ve? — m*)/m = Re™ 10) 


such that the ¢-plane is mapped into the interior of a 
unit circle in the ¢-plane with the surface of the wing 


lying on the circle |¢| = 1 (upper surface on lower 
semicircle) and the Mach cone on real axis extending 
from (1 — V1 — m?*)/m to (-—1 + V1 — m*)/m 


(see Fig. 4). 
Thus, 


er {3 ma [e(@— VE — mT 
m - 


n=1 
R.P. (i D ben—1t?" ) 
n” 1 


(R <1) (11) 


ll 


=— pe bo, R?2"—! sin (2n — 1)¥ 
n=1 


1) is 


Equation (11) clearly indicates that condition 


satisfied. For condition (2), one needs to know the 
nature of do,_;. Putting 
¢ = y/x = msin 0 (12) 


one can observe from Eq. (8) the lift distribution on the 
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top surface of the wing is 


u = amr+a’ sec6+ 
DY  den-1(—1)"—! cos (2n — 1) (13) 


n=1 
or 
u — a’ sec 0 = dor + 


>» ben-1(—1)"—! cos (2n — 1)6 (14) 


n=1 


It is immediately seen that dp, be,-1(—1)"~! are co- 
efficients of the Fourier cosine series representing an 
arbitrary curve u — a’ sec 6, bounded in the interval 


0 <@<2-7. Hence, the convergence of the infinite 
series is proved. 
Assume that the curve u — a’ sec @ represented by 


Eq. (14) is bounded, continuous and with finite number 
of maxima and minima. Then the series is uniformly 
convergent. This permits one to use for approxima- 
tion a finite number of terms to represent the lift dis- 


tribution wu. In fact, if « — a’ sec 6 has merely a 
bounded first derivative, }:,_; will decrease with in- 
creasing like (2n — 1)~?, which indicates a fairly 
rapid convergence. A four-term approximation for 
u — a’ sec 6 gives for the lift distribution on the top 
surface of the wing: 


Uu = apo + a’u’ + by, + dguzs + bsus (15) 


wow 


where 
uy = 7 (16) 
u’ = m/Vm? — & = sec 6 (17) 
Un = R.P. fil(¢ — Ve? — m?*)/m)}-"! 
= (—1)*"*-1 cos (2n — 1)0 (18) 


Or, one may split 7; such that 


/ 9 , Ma 9 
uy = Vm? — &?/m = [m/(2V m? — &?)) + 


[(m? — 2¢2)/(2m~ m? — ¢7)] (19) 
and 

U = aAotteo + by'uy’ + B,"uy” + dgug + bsus (20) 

where 
uu,’ = ™ (2V m? — {?) (21) 
uy” = (m? — 2¢7)/(2mNv ‘m? — ¢) (22) 
b,’ = 2a’ + by (23) 
b,” = (24) 


Substituting Eq. (20) into Eq. (4), one obtains the total 
lift coefficient, 


CL = mw (4a9 + by’) (25 


instead of (4a) + 2a’ + b,) which would be obtained 
if Eq. (15) were used. This facilitates the succeeding 
manipulations a great deal, including the avoidance of 
certain integration difficulties. 

The upwash distribution compatible with the lift 
can be evaluated from Eq. (3)—viz., 


W = ApWo + byw,’ + b,"w,” + byw; + Dsws (26 


ee 1+ met 1 — m¢é 2 I = 
where W = bh Naat. | cosh . : + cosh~! ; | lee, ital (24 
m m+ ¢ m= > m ¢ 
w,’ = —[1/(2m)] E(V 1 — m?) (28 
l / : 3 ee. é 29 
Wy" = re V1 — ¢? — cosh"! + K (V1 — m?) - - E(V 1 — m?) (<9) 
l m* + 4 ‘ 
= * | E 1(V1 — ¢)3 — 38m V1 — &2 + Bm cosh-'!| ~| + — 5m K(V 1 — m) | (30) 
m? é m 
5 {16 , WM» : | 
Yo — i: (V1 — £7) + | 4m? — 3 (V1 — £7)? + mtv 1 — & — m‘ cosh"! iz 
m° ) k c 


l 
15 


[((32 — 52m? — 3m4)E(V 1 — m*) + (—16m? + 39m4)K(V 1 — mt 


(31) 





The ex 
Equ 
With § 


Cb _ 


where . 
are fur 


leading 


Hence, 
The pai 


m 
[ UjWo 


ry 


where 


(The: 





Airly 
for 
top 


(15) 


(16) 


19) 


20) 


9) 





m P ] aie dap T 2 
| uwdt = 24 (1 — V1 — m?*) log m — log (1 + V1 — m?) + c — m*) sin~! m+ cos>!m — “lt (35) 
0 1 
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The expressions for w), w;’, w,”, and ws are obtained, directly or with modifications, from a previous paper.° 
m 
Equations (20) and (26) suffice for the evaluation of f u(¢)w(¢)dt in Eq. (5) for an arbitrary lift distribution. 
0 


With Eqs. (7) and (25), it can be shown that 

+ ‘(= Boy rt { s T / h,’2 4 B Boy us 4 b,’E » 
= . om £ ‘y’ 4 — 2 tes ne = aes — 2 i] ) 
Cp m W416 4 mo 16 ] m 1 "1 ' S Vi m* JO, Oy + 


By: Bos 
(B.. = ") by’bs + (4. = \ ) by’bs + (Ava + = V1l- m?) by”? + (Byrs) 6:%b3 + (Brr5)b1"b5 + 
- ) 


P i Avo Bow C. , Bow C, 
(A 33)b3? + (Bys) bbs + (A 55)b5” - =e T = b, + b,” + 
8 4 T 4 T 
me (S)}n [Be(e)]ue (2) om 
4 T 4 T 16\nr/$ 


m m m m 
where A,; and B,; stand for [ u,w,dé and (f u,wdge + used and A;-; stands for [ u,/w,d¢, etc., which 
0 0 ZU ~~ 


are functions of m. In the above expression, Cp, is included. As a’u’ is the only component of u which produces 
leading-edge thrust [see Eq. (15) ], it comes naturally that 


C,’ = 2xa’ = x(b,’ — 5,") (33) 
Hence, 


Co, = [(#V 1 — m?*)/(4m)] (b)’ — 3”)? (34) 


m m 
The parameters f u,w,de¢ and f u,w,d¢ are tabulated as follows: 
0 0 


n 2 
ain V1 — mm? 7 
Uy 'wde = I, — Iz (36) 
where 
a/2 9 
2 sd 2 sin? sl = 37) 
=2 log (V1 — m? + V1 — m? sin? 0)dé + = log a (34 
we T 2 —_ 
= log (1 + V1 — m? sin? 6)d@ + > log : (38) 
0 é ? 
(The integrals in J; and J] may be expanded into series. ) 
m ~ 
[ uw, de = — E(V 1 — m?) (39) 
[ u,"wde = —E(m) (40) 
~/ 
m af{[m ew ae ts 
[ uyw,"dt = Ek V1 — m* — m ilog (: le at * — 1m + cos~! m + 
J0 m S m é 
: m ~ ra 
aE(V 1 —- ww) - =~ Vi-— 8) —-s% (41) 
m 
f U3swWod ¢ = 2[E(m) — So]* (42) 
0 


* Sx = f sin @ ‘V1 — m? sin? @d0. For evaluation, see Appendix II. 
0 
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=m 
j | 
U pwd = FT 


0 | 2m? 


ms 
i! ‘ i+ Vl — ‘ . ] 

—- mV 1 — m* + 2m log = Sart" ne 2cos-'m + mr] 
2m nl a 





| 2m 1 — m*)? + 3m V1 — m2? 4+ 3 sin-'m 


m* + 4 ; { 
— 2 3 m - JA(V 1 m*) ¢ 
m-* 
ey, »») 16 
521 ¢ 2E(m) 4 ee S 
P 3 3 
j F5 my | zi lee 
My Ws ¢ lis [ Sm (V1 — m*)? + 10m(V 1.— m?)? Lom Vl — m? + 15 sin m | { 
Jo m°? 3 


- ) 
(Wi 

( rs [ 200 (V1 — m*)* + 3m V1 — m? + 3 sin m | 
” 


m+ : 1+ Vil — m-? ) : 
> |m V1 — m* — 2m log + Si" 9 = 2 Cos * 7 | 
m 


m ; : | 
[ (32 — j2m°" 3m4)E(V 1 — m2) + (—16m?2 + 39m4)A(V | m- 


mo . 
uw, de = —-— E(V 1 — m?) 
0 S 


mn 


| uy" de = 0 


=m 


make l : T , i ) 
u, Ww, de = “= } E(m) = = [K(. l— we) —- (V1 - m?) | - Ey 


om 


U3 de = () 


« 


u— ; ly | . x| (/m?> + 4 : ? ‘ ) 
W,'w;de = 4 — _ + 3) E(m) + = E(V 1 — m*) — 5K(V 1 — m2?) | + 482 + 3lee 
Jo 2I m° 2 m* f 


=m 
| usw, 'd¢ = 0 


ing 5 416 16 
| u,'wde = [E(m) — 2m?S. + m'4S;) + (4m: — ) [E(m) — m?S2] + m*|E(m) — [2] + 
v7 ” 


Imi t5 


T : : ( 
~ 30) [ (32 52m? — 3m4)E(V 1 — m*) + (—16m? + 39m1)K(V1 — m:) |e 
m ; = 
uy"w,"de = 2 [—E(m) + S>] 
f uzw,"de = E(m) — 2S.+ S; 
ii ” 3 . - ‘ . 
Wy’ w3de = = [E(m) — 5Se + 454] 
. ” ‘ l . . ‘ ‘ 
usw,’ de = — : [BE(m) — 14S. + 19.8; — 855] 
«7 ” 
an ” oS) . . . sal 
u,"wde = — = [E(m) — 13S: + 28S; — 16.86] 
f usw3de = — 3[E(m) — 6S. + 9S, — 456] 


f uswede = 5S[E(m) — 14S. + 41S; — 44.8, + 16°] 


+) 


50 


(52) 


(59 





of the 
with re 
lineal 
propria 
All the 
merical 
ing lift 
for any 
In th 
that th 
tance 
in the ¢ 
this thr 
total di 
siderati 
flow sey 
tions ol 
Cy 'C? 
f¢ rees ¢ 
sharp Ie 
been vé 
the pre: 
to the 
nomeno 
plained’ 
ring at : 
rates al 
edge. 
Atten 
fore leac 
based o 
suction. 


are inve 


(a) Suc 


Full t 
assumed 
used for 


CC MITES] 0) 


(b) Suc 


The si 
Results 
the sam 
plies an 
be appre 
no sepal 
in the r 
the follo 





i+) 





SUPERSONIC CONICAL 


u;w3d¢ = 3E(m 


, - ao . 
U;wW;a¢ = - [3le(m) — 50S. + 


9 
« ”) 


Most of the integrations are lengthy, but about half of them 


in Appe ndix I. 


DETERMINATION OF MINIMUM DRAG AT CONSTANT LIFT 


Equation (42) shows that Cp is a quadratic function 
of the coefficients },’, b,”, 63, and b;. Differentiation 
with respect to each of these coefficients leads to four 
linear simultaneous equations that determine the ap- 
propriate set of the coefficients giving minimum Cp. 
All these coefficients are multiples of C Thus, nu 
merical values of minimum Cp /C;? and the correspond- 
ing lift distributions (7) in terms of C; can be obtained 
for any particular m. 

In the determination of minimum Cp C 


that the role of leading-edge suction force is of impor 


, it appears 


tance. This is expected, since it has been illustrated 
in the case of a flat delta wing that the magnitude of 
this thrust may amount to as much as one-half of the 
total drag for very small aspect ratios.*. Some con- 
sideration is thus given to the question of leading-edge 
flow separation that may greatly reduce the contribu- 
tions of leading-edge force, and hence, the values of 
Cp C2. In fact, the theoretical effect of suction 
forces can hardly be taken into account for a wing with 
sharp leading edges moving in a real fluid. This has 
been verified by experiments’ " that indicate that 
the pressure peak at the leading edge is smoothed due 
to the presence of local flow separation. This phe- 
nomenon occurs for sharp leading edges and is ex- 
plained'® as the result of the ‘‘thin-airfoil stall,’’ occur- 
ring at small angles of attack, in which the flow sepa- 
rates and reattaches to the surface near the leading 
edge. 

Attention to the effects of flow separation will there- 
fore lead one to different answers for minimum Cp /C,?, 
based on various assumptions concerning leading-edge 
suction. In the present analysis, the following cases 
are investigated. 


(a) Suction Force Included 


Full theoretical effect of the leading-edge suction is 
assumed to be present. Thus, Eq. (32) may be directly 
used for the evaluation of minimum Cp/C,? and the 
corresponding lift distribution. 


(b) Suction Force Omitted 


The suction force is assumed to be totally ineffective. 
Results are obtained by omitting Cp, and carrying out 
the same calculations as in (a). This, however, im- 
plies an assumption that the actual lift distribution can 
be approximated by the theoretical one (as if there were 
no separation) without introducing appreciable error 
in the results. This approximation is motivated by 
the following observations: (1) experiment" indicates 


26S. + 75S; 
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S48, + 32Ss 60 


235.58; — 460.8, + 4008, 128.Si0 | 61 


can be obtained from a general formula which is derived 


that at sufficiently small angles of attack, the lift dis- 
tribution agrees well with the theoretical curve except 
for a small region very near to the leading edge and 
Dh , 


(2) for this case, the appropriate value of a’ | | ' 
b,”)|, the only component of « which shows a pressure 
peak at leading edge, happens to be very small in com 


parison with the same in case (a 


(c) No Suction Force 


The lift distribution is so restricted that its intensity 
at the leading edge is finite. Thus, theoretically, there 
will be no suction force developed Results can be 
obtained by putting 0,’ b,” b; and going through 
the same calculations. Values of minimum Cp/C,? 
obtained for this case, if compared with those of (a 
will give an idea to what extent the suction force 
actually contributes to drag reduction, as compared 


to the best available distributions without such forces. 


DETERMINATION OF OPTIMUM CONICAL WING SHAPES 


For any arbitrary w({), one may obtain the wing 
shape z x = h(¢) from the boundary condition 02 Ox = 


w(¢)* which is equivalent to 


h(c) — th’(c = WIC (62 


with ¢ = y/x. Thus, the optimum shape that pro- 


duces minimum Cp/C,? is readily known: 


X/ opt ( "i f 


where [w(f)].pt. is obtained by substituting the ap- 
propriate coefficients determined in the last section 
into Eq. (26), and A is an arbitrary constant. The 
presence of A in the equation indicates that a unique 
solution cannot be obtained by the linear theory. 
This, in fact, is true, since, according to linear theory, 
a flat-plate delta wing with dihedral at zero angle of 
attack creates no disturbance velocities. Neverthe- 
less, A should be such that the assumption of small z 
everywhere on the wing is not violated, since it is one 
upon which Eq. (63) is based. Using similar nota- 
tions as those for u and w, one may write, for con- 


venience, 


2/X)opt. = @o(So/x) + by'(a"/x) + 
by” (2," x) + bs(23/x) + bs(25/x) — K¢ (64) 
where do, },’, b:", 63, and b; are the appropriate set of 


coefficients and 


* The method of deriving wing shapes from w(¢) has appeared 
in reference 9, with slightly different approach. 
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Equations (65), (67), (68), and (69) all contain a singu- 
lar term, cosh “1/1 tae It is expected that (2/x)opt. 
will have an infinite value at the centerline of the wing 
due to the use of finite number of terms for approxi- 
mation. This undesirable feature, however, can be 
avoided by adjusting one of the appropriate coefficients, 


say bs, so that the terms containing cosh~! |1/¢| will 


add up to zero. Calculation shows that the error 
introduced in the value of minimum Cp/C,? by such 
an adjustment is negligible. In Fig. 7 are plotted to- 
gether the optimum shapes with and without the singu- 
lar terms. The values of K are so chosen that the ele- 
vation of the leading edges is the same as that of the 
centerline for those shapes with the singularities re- 
moved and that the elevation of the leading edges is 
unchanged from the previous case for those shapes with 


the singularities kept. 


It is to be emphasized that the optimum shapes ob- 
tained from the five-term approximation give values of 
minimum Cp/C,? which are sufficiently accurate for 





id 


5 oe 
(—16m? + 39m)K( V1 — m) lt -—e cosh—! | — (69) 
i 


engineering purposes, although they cannot be con- 
sidered as the exact optimum shapes. 


RESULTS 


Minimum values of Cp/C,?, the corresponding lift 
distributions and wing shapes, under various assump- 
tions have been computed for m = 0.5, 0.8, 1.0. Flat- 
plate data are included for comparison. 


(1) Minimum Drag Coefficients 


TABLE 1 
7 Minimum Flat Plate 
Case Cp/Cr? m(Cp/Cr?) Cpo/Cr? = m(Cp/Cz*) 
= (a) 0.0796* 
=o 0. 1592* 
(a) 0.2476 0.1238 0.2477 0.1239 
m =0.5 (b) 0.2486 243 0.3855 0.1928 
(c) 0.2491 
(a) 0.2214 0.1771 0.2224 0.1779 
m =0.8 (b) 0.2216 0.1773 0.2821 0.2257 


0.2221 


(c) - 
—e (a) and (b) 0.2298 0.2500 
oe ee 


0.2298 0.2500 
0.2310 





* Obtained from reference 2. 
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(2) Lift Distributions 
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(3) Camber Shapes 
; 7 : SUCTION INCLUDED i 
> 2 lined an ~ + 729 = aus RETAUe 
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“ 29 
ees = SUCTION OMITTED BS 
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Se we a0 a $ : = 
z SINGULAR TERMS RETAINED ~*-~~_ a z 
a2 | SINGULAR TERMS REMOVED oS vid s 
¢ aa - eo ae 3 = 1075 6 oa 
Fic. 7. Optimum conical cambers (based on five-term approximation) at JJ = V2 
CONCLUSION AND DISCUSSION lift distribution will contain only a small component 
fa , as Spe ea. oe that would otherwise produce suction. It has even 
The theoretical minimum of Cp/C,? for a_ thin, ; : 
been conjectured that this component would com- 


cambered, conical, delta wing with subsonic leading 
edges, including the full effect of suction forces, is found 
to be very nearly equal to the value for a flat plate.* 
Only when the leading edges are very close to the Mach 
cone is there noticeable improvement in drag by use of 
the optimum shape. The maximum gain in the limit- 
ing case of sonic leading edges is still only 8 per cent. 
The conclusion may be drawn that in this case no ap- 
preciable reduction of the induced drag of a flat delta 
wing will be obtained by using ‘‘conical’’ camber. 

However, the theoretical calculation sometimes falls 
short of the truth for wings moving in real fluid. Un 
less the leading edge is sufficiently rounded, there will 
be flow separation that nullifies the suction forces. As 
a consequence, the theoretical minimum Cp/C,” esti- 
mated by assuming full effect of suction forces is not 
always attainable, and the corresponding optimum 
camber may not be the ‘“‘genuine’”’ one. 

The results for minimum Cp/C,? with suction force 
omitted show an appreciable improvement over the 
corresponding case of the flat plate (with suction force 
omitted). In fact, the minimum values without suc- 
tion are only very slightly higher than the minima of 
the former case with suction. In other words, the use 
of proper camber will almost entirely overcome the 
unfavorable effect of the loss of suction force. 

The lift distribution for minimum drag with suction 
omitted includes only a very small component of the 
type that would produce suction if flow separation were 
absent. 
proportional to (b;' — 6,”)?, and the optimum lift dis- 
0.8) leads to |b,’ — 5,”"| = 
minimum (case of no 


The suction force, according to Eq. (34), is 


tribution (say for m = 
0.0052 C,/r. A 
separation) was also determined with the condition },’ = 


“restricted” 


b,”, and the results show very small differences from the 
foregoing case, where the restriction was not applied 
but the favorable suction force was omitted. 

It appears natural that if the suction force is not 
taken into account, then the corresponding optimum 


* This result agrees with a recent investigation made by Adams 
and Sears® for relatively slender plan forms, as well as the well- 
result of R. T 


known Jones and others for very slender plan 


forms 


pletely vanish for the ‘‘exact’’ minimum if the approxi- 
mation involved by the use of a finite number of terms 
representing the lift could be avoided. This specula- 
tion is practically of little importance, since the dif- 
ference between the ‘restricted’? and the unrestricted 
case is negligible, even in the present approximation. 

The same method of approach in the determination 
of minimum drag for conical wings could be extended 
to nonconical wings by proper choosing of the w and 
w series. For illustration, the simplest case of finding 
the optimum combination of a parabolically-cambered 
wing and a flat plate has been treated in Appendix III. 
The results show that, for cases with suction included, 
there is practically no improvement in drag (maximum 
reduction is only 6 per cent at m = 0.4),7 while for 
cases with suction omitted, the improvement is large 
for small m (maximum is 33.4 per cent at m = 0) and 
is small for large m. These data give some idea about 
the effect of chordwise camber and its comparison with 
conical camber. 

The theoretical optimum conical shape for minimum 
Cp/C.2 by using finite number of terms for approxima- 
tion shows a logarithmic singularity at the centerline 
of the wing. However, the removal of the “singular 
part’ of the shape does not influence the Cp /C,” essen- 
tially; this may be understood, as the singular part 
of the shape acts like a central ‘“‘fin’’ with very small 
effects on the pressure distribution. By calculation, 
one finds that the effect of the singular part on mini- 
mun Cp/C,? is practically negligible. 

Althouch the values of minimum Cp/C,? 
corresponding pressure distribution are insensitive to 
the number of terms taken for approximation, the shape 
of the center portion of the corresponding wing 1s, 
other 


and the 


under certain conditions, sensitive. In 
words, the shape of the center portion of the wing has 
a small effect on Cp/C,? and on pressure distribution, 


and the approximated “‘optimum’’ wing may still show 


very 


very bad convergence near its center, while the number 


the improvement is zero. agrees with 


TAt m = 1, This 
Graham’s conclusion’ that, for delta wings with sonic or super- 
sonic leading edges, a simple camber produced by bending the 
wing along spanwise lines gives no reduction in drag 
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of terms taken is considered sufficient to give close value 
of minimum Cp /C_?. 

The case in which the worst convergence for the shape 
has to be expected is the one with suction force omitted 
and m large. This was concluded from the fact that 
camber is very ineffective in the drag reduction for this 
case, as indicated by Fig. 8. The reason may be 
stated as follows: The ‘‘conical camber”’ is actually a 
specific combination of chordwise camber and _ twist, 
and here it is the twist that is more effective in the 
conical case. Eut, for the same slope of the wing ordi 
nate, the twist is largest near the tip and smallest at 
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the center. Hence, the shape of the center portion may 
exhibit great variations and will not be determined 
sufficiently if only a few terms are taken for the evalu 
ation of minimum Cp ‘C,?”. 

To illustrate the accuracy of the five-term approxi 
mation, values of minimum Cp’ C,”? for the case m 
().8, with suction forces omitted, have been calculated 
They are, in the order of increasing number of terms 
taken for approximation, 0.2821, 0.2420, 0.2257, 
0.2225, and 0.2216. For other cases, the convergence 


is expected to be even better 


APPENDIX (1) 


EVALUATION OF U2, —1Uon! 


Then 
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Fic. 8. Minimum Cp/Cx?2 and minimum Cp/C ,? multiplied by leading-edge slope for wings with parabolic chordwise camber (constant 
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Comparison with flat plate 
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Uae: = (—])"*—!m S cos (2n — 1)@ cos 6dé 
= (—1)"—! (m/4) ((sin 2n6/n) + } [sin 2(7 1)A| (n 1)}) \3 


From Eqs. (3), (12), and (18) 


dWey,1/d¢é = (—1)"—! [(2n — 1)/m?] } [cos (2n — 1)6]/(sin 6 cos 6)f V1 m?* sin? @ \4 


Since U2, 1 1s zero at ¢ = 0, m (or 6 = 0, r/2) and w»,-_; shows only logarithmic singularities at ¢ 0 and is finite 


at ¢ = m, the first term in the right-hand side of Eq. (A2) vanishes. Thus, 


m mm dt . 
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Jo sin ¢ 


So, which are evaluated in Appendix (IT). 


APPENDIX (11) 


EVALUATION OF S», = [ sin*9@ V1 — m? sin? 6 dé 
Since 
d . . . . . l . 
16 (sin*—! @ cos 0 V1 — m? sin? 6) = (2k — 1) sin*®-26 V1 — m? sin? 6 — 2k sin*® 6 V1 — m? sin? 6 — 
¢ 
m? sin**@ m? sin? +2 @ 
(A6) 
V1 — m? sin? 6 V1 — m?* sin? 6 
One has, by integrating both sides from 6 = 0 to 6 = 7/2, 
ae f sin’ @ V1 — m? sin? 6d0@ = (2k — 1) [ sin*-29@ V1 — m? sin? 6 dé — 
v0 ~ Wt 
~“—* sin?’ @ lat sin2*+2 g 7 
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= dé = — = + = = 
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Hence, after rearranging terms 
Soe = 4 1/[(2k + 1)m™*]} (1 — m?) [K(m) — E(m) — m?S, — m'S;,... — m*-4S», s] + 
m*—2(2km? — 1)Sy%-0f (AQ) 
from which one can obtain 
So = [1/(38m?)] [1 — m2)K(m) — (1 — 2m?)E(m)] (A10) 
Sy = [1/(5m4*)] } (1 — m?) [K(m) — E(m)] + m2(4m? — 1)S,! (All) 
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Sw = [1/(11m")] )} (1 — m*) [K(m) — E(m) — m?S, — mS, — m®S. + m>(10m? 1).S5 


APPENDIX (III) 


EFFECT OF CHORDWISE CAMBER ON DRAG 
The same method of approach for determining the minimum C) ‘C,? value of a conical wing can be extended to 
This is illustrated by the following typical case—a delta wing with parabolic camber constant 


nonconical wings. 
2. The results obtained will provide one with a rough idea of the relative 


along span (2 « x? + constant) at WZ = y 
merits of chordwise and conical cambers in drag reduction. 


Similar to the conical case, one may put 


uU = Ayu, + Aolto Al15) 

W = aw, + AsWe (A16) 

where “, w = the resultant nondimensional perturbed velocity components on top surface of the wing; ™, w, = the 
= the components created by a wing with upwash at surface proportional to x; and 


77) 


flat plate components; uo, w» 
= arbitrary coefficients that are linear functions of C,, coefficient of total lift. 


uy, = m/Vm? — & A17) 


a), Qe From references 15 and 9, 


= —(1/m)E(V 1 — m?) 


op 
oO, = 


No, W. can be deduced from reference 18, such that, at J = 7/2, 
u2 = (x/CH) [(2m? — &*°)/Vm? — &? A19) 
We = x/C (A20) 
where C = root chord and 
H = [(1 — 2m?)/(1 — m?)] E(V 1 — m?*) + [m2/(1 — m*)JK(V 1 — mm: 
It can be shown that 
; a" Ps m 
Cy, = | baw gg Mt Ons Sina a | A21) 
m 3H 217° 
for a i _ 2 4 KF =. oe?) 
— | 2E( V1 m°*) V1 m | a2? + = ofp 2E( V1 m*) 
i m 0 hed 
4V1 — m? m ae 
‘vil : m ayaz + - (3 sis as ) at (suction force included) (A22) 
3H 2H H 
and 
IF a: 4 IF —. 3m ) ; ' ; 
Cs = 7 J2E | elt. a a;? + : a 2E (WV 1 = | aw: + ath (suction force omitted (A23) 
{ m 3 HT 2S 
at a lift coefficient 
Cr = 2rfay + (m/Hjae] (A24) 


With Eqs. (A22), (A23), and (A24), one can determine the appropriate set of a; and a» for each case by certain mini- 


mization procedure. 
Results are plotted in Fig. 8. 
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An Investigation of the Bending of Angle 
Beams in the Plastic Range 


H. NORMAN ABRAMSON,7 HARRY A. WILLIAMS,£ AND BRUCE G. WOOLPERT** 


Stanford Unwerstty 


SUMMARY 


An earlier investigation'® revealed the necessity for further 
study of the bending of unsymmetrical beams in the plastic 
range, inasmuch as twist occurs and has a marked effect. 

Angle cross-section beams of 24 S-T aluminum-alloy material 
were subjected to bending in the plastic range with the plane of 
loading at various angles to the minor principal axis of the cross 
section. The moment was applied by means of an eccentric load 
with a relatively long moment arm. 

For a simplified approach, it was found that Cozzone’s method 
of plastic bending analysis® and an analysis based on an exponen- 
tial relationship between stress and strain both gave reasonable 
bending 

Both 


methods are approximate since they assume that the neutral axis 


correlation between theoretical and experimental 


moments in the plastic range beyond the yield strength 


does not rotate or translate and the cross section does not rotate 
The exponential method is not valid in the range below the yield 
strength. 

Other factors concerned with the general behavior of angle 
beams subjected to plastic bending are also considered. 


SYMBOLS 


a = constant dependent upon the material for the 
exponential equation 

a- = distance to point on neutral axis, in 

Aa, = change in dg, in. 

b = length of short leg, in 

( = general symbol for distances from neutral axis, 
in 

Ce, Ct = distance from neutral axis to outer compressive 
and tensile fibers, respectively, in. 

d, g = symbols for particular cross-section dimen 
sions, in. 

e = general symbol for strain, in. /in. 

Oc, ef = outer fiber compressive and tensile strain, re 
spectively, in./in 

e; to ey = strains determined experimentally, in./in. 

Ge = stress determined from compression and tension 
tests, respectively, lbs. per sq.in. 

fos ty = proportional limit stress and yield strength, re 
spectively, from average stress strain curve, 
Ibs. per sq.in. 

fn = outer fiber bending stress, lbs. per sq.in 

‘a, Ie = outer fiber compressive and tensile stress, re 
spectively, ibs. per sq.in 

fo = Cozzone’s intercept stress, lbs. per sq.in 

locs Jes = intercept stress for compressive and _ tensile 
area, respectively, lbs. per sq.in. 
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k ,k = Cozzone’s beam section constants 

K = constant in exponential equation, Ibs. per 
sq.in 

Sn» 9 distance to point on rotated neutral axis, in 

t = Young’s modulus of elasticity, Ibs. per sq.in 

| RE Gn! oe moments and product of inertia of cross section 
with respect to x and y axes, in.* 

I,, Iy = moment of inertia with respect to principal 
axis, in.4 

Ty, I, = moment of inertia of compressive and _ tensik 
areas, respectively, with respect to neutral 
axis, in. 4 

Ky = plastic stress function based on average stress 
strain curve, lbs. per sq.in 

Kye, Ky = plastic stress function based on compression 
and tensile stress-strain curves, respectively, 
Ibs. per sq.in 

iL = moment arm to centroid at mid-span, in 

VW, = experimental bending moment, in.-lbs 

Mu = bending moment computed analytically, in. lbs 


Win, Mer = bending moments computed analytically for 


compressive and tensile areas, respectively, 


in.-lbs 

n = power to which e and p are raised in the ex 
ponential equation 

r = dynamometer load, lbs 

Qmer Ome = statical moment of compressive and_ tensik 
areas, respectively, in.* 

r = distance from neutral axis to an element in the 
beam cross section used in the derivation of 
the exponential equation, in 

t = thickness of angle beam leg, in 

Al = Aa, sin y, in 

a = angle between principal axis and neutral axis, 
deg 

Aa = change in angle a, deg 

0 = angle between plane of loading and y axis, deg 

@ = angle between plane of loading and normal to 
neutral axis, deg 

p = radius of curvature, in 


= angle between x-axis and neutral axis, deg 


INTRODUCTION 


-— BENDING in the plastic range when the applied 
loads are parallel to a principal plane of bending 
has been studied by a number of investigators.’ 

When the loads are not parallel to a principal plane, the 
complexity of the problem increases many fold. The 
complexity is further increased when the beam cross 
section is not symmetrical with respect to a principal 
axis. Goodier’s general solution of the problem for 
open sections! indicates the difficulties encountered 
even in the elastic case. It is seen at once that pure 
bending, as it is usually understood, does not prevail 
except under very restricted conditions.¢7 If the ap 


+t This general subject is discussed further in another paper 
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— Plane of Couple 
eo a =~ y 
d ha yf 
y 
y NA 
a | 
(c) 
SECTION PROPERTIES 
Rectongulor Axes Principal Axes 
I, + 0.2706 Ig = 0.2854 
I, = 0.0456 lj > 0.0308 
ln, = 0.0600 
] x C) 1, l. GG. | Gis oF Ce 
o° 66°33' | 52°35' [0.0330 | 0.0385 |0.0929 | 0.0862 | 0.670 | 0.744 
30° 83°37' | 39°39' |0.0167 | 00173 /00635 | 0.0622 | 0.507 | 0.518 
60° 88°13' | 14°15' | 0.0140 | 0.0177/0058! | 0.0558 | 0.457 | 0.572 
76° 07' | 90°00'|-0°07' | 0.0127 | 20183) 0.0570 | 0.0537 | 0.434 | 0592 
FIG. | — ANGLE BEAM CRCSS SECTION PROPERTIES 


plied torsion and axial thrust are eliminated from Good- 
ier’s solution, it is evident that a torsion couple must 
still act at any cross section except the midpoint of the 
length unless all of the deflection is parallel to a principal 
plane. This couple produces twist, and it is evident 
that the relationship between components causing 
normal stress must vary along the length of the beam. 
Hence, the neutral axis must rotate also. 

If plastic action prevails to any large extent, the vari- 
ous parts of an open section must be relatively heavy 
and the usual assumptions for shear flow might be 
questioned. The major obstacle to the solution of the 
plastic bending problem results, however, from the 
This 


leads to even more complex relationships between 


nonlinear relationship between stress and strain. 


moment and stress and eliminates the use of super 


position. Except in special cases, pronounced curva- 
ture results in considerable rotation. 

It is evident that a rigorous solution of the problem of 
plastic bending of unsymmetrical beams must be de- 


9 


veloped from the concepts of pure plasticity’ e.g., 
Phillips presents the exact theory of combined pure 
bending and axial force for a prismatical bar’ using the 
The 
solution of the differential equations is extremely dif- 


equations of the theory of plastic deformations. 


ficult. Hence, it appears wise, for the benefit of the 


designer, to seek approximate but practical solutions 
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which will give reasonably accurate results. Some at 


tempts in this direction have already been made.'*'® 

In a previous investigation, '® by one of the authors, of 
plastic bending when loads were not parallel to a prin- 
cipal plane, some exploratory tests of angle beams were 
made. A modification of Cozzone’s approach’ gave a 
resisting moment at the mid-span which was less than 
15 per cent larger than the applied end moment until 
plastic instability occurred at which time the per cent 
difference increased rapidly. 

The present investigation is an extension of the pre 
vious work. In planning the program it was decided 
to test beams of rather long span in order to accentuate 
rotation and to apply the end moments by means of a 
small eccentric load with a relatively large moment 
arm.* The primary purpose is to determine to what 
extent a relationship between bending moment at the 
mid-span and outer fiber stress will agree with the 
applied end moment if a simplified theoretical approach 
is used in the analysis. The resisting moment is com 
puted by a modification of Cozzone’s graphical method 
and by an analytical method based on an exponential 
relationship between stress and strain. Both methods 
neglect rotation and other variables in the interests of 
simplicity. Data on rotation of the cross section and 
movement of the neutral axis are also presented to add 
information regarding the general behavior of the tested 


beams. 


APPARATUS AND TEST PROCEDURES 


Test Specimens and Stress-Strain Characteristics 


The test specimens were machined to the dimensions 
shown in Fig. 1 from 24 §-T aluminum-alloy stock, 1 by 
2 in. 1n cross section. Each beam had a clear span of 24 
in. Photographs of some specimens after testing are 
shown in Fig. 2. 
stress-strain 
S-T 


Considerable difference between the 
curves in tension and compression exists for 24 
aluminum alloy, but inasmuch as this material is com- 
monly used in practice, it was decided to use it rather 
than some other material which might have tension and 


All 


analytical computations were based upon an average 


compression curves that are in better agreement. 
stress-strain curve. Tension and compression samples 
from the original stock gave stress-strain data as in 
dicated in Fig. 3. The curves of Fig. 3 were selected as 
being representative of the material and were used for 
is also 


all computations. Cozzone’s intercept stress, fo, 


shown in the figure. 


Loading and Strain Apparatus and Testing Procedure 


The arrangement of the loading apparatus used for 
bending the angle beams is shown in Figs. 4and 5. The 
two lever arms, machined from 3, 4-in. steel plate, were 
supported at one end by cables from an overhead beam 

* As used in this report, rotation refers to twist and not to 


changes in end slope of the longitudinal axis 
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FIG. 5— PLAN VIEW OF PLASTIC BENDING 
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Fic. 6. Location of electric strain gages, 


and at the other end by ball bearing rollers. Tracks for 
the rollers were attached to a heavy table immediately 
below the suspended arms and so arranged as to con- 
strain the movement of the fore-end of the arms in a 
path perpendicular to the center line of the apparatus. 

The arms at the cable end were provided with fittings 
which securely held the specimen in a position perpen- 
dicular to the two arms and parallel to the tracks, thus 
forming a rigid frame. The specimen fittings were 
further designed, by means of a movable dise which 
could be secured with a retaining ring, so that the speci- 
men could be rotated about its lengthwise axis to pro- 
vide any angle of bending desired. The load was ap- 
plied through cables attached to the arms as shown. It 
is obvious from the loading arrangement that pure 
bending was only approximated, due to the presence of 
a small axial load. 

Strain readings were taken with an SR-4 type indi- 
cator which is shown to the left of the bending apparatus 
in Fig. 4 The strains were measured at various points 
around the center cross section by SR-4 type A-12 elec- 
trical strain gages cemented to the specimen, as shown 
typically inl Fig. 6. 

The specimen was placed in the bending apparatus 
and adjusted to the desired angle of rotation. The 
normal testing procedure was to apply an increment of 
load, stop the testing machine, and immediately record 
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the strain indicator readings; particular attention was 
given to the latter step, when in the plastic range, be- 
cause of creep. 

As the plastic range was reached, the decrease in 
moment arm resulting from the curvature of the beam 
was measured, at the same time that the strain readings 
were recorded. Rotation of the beam was also measured 
at that time. 

Specimens were removed when the reduction in 
moment arm was sufficient to tend to prevent any in- 
crease in bending moment when additional load was 
applied. Generally, one or more of the electric strain 
gages had become inoperative at this point, and the 
maximum outer fiber strain was on the order of 2 per 


cent. 
ANALYTICAL PROCEDURES 


Cozzone’s Method 


From reference 16, Cozzone’s method® may be modi- 
fied to apply to an angle beam if the rotation of the 
beam in the plastic range is neglected and the neutral 
axis is assumed to remain stationary. Neglecting these 
factors, an approximate solution may be found by 
dividing the cross section into a compression area and 
a tension area. The resisting moment of the compres- 
sion area is obtained by assuming the symmetrical cross 
section shown in Fig. l(b). The resisting moment of 
this beam is computed and divided by two. This is 
equivalent to computing /y, for one-half the cross sec- 
tion only. Then k, = [Qn./(/y,/¢-)] and the moment 


becomes 


Mi = (Inxe/€e)(1/c08 ) [fme + foclke — 1)] (1) 





Plane of Couple 





For 24 ST Al Specimens: 


a = 91,450 psi 
n= 1.1228 
k = 25 psi 
Cross-section notation as used in the exponential 
solution for Mz, 


Fic. 7 
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Table No. 2 


BEAMS 
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- THEORETICAL BENDING MOMENT 


Angle Beam M- 4 
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In the case of angle beams, an expression for bending 
moment may be found as follows: 


l a l 1 
V, n—1 (d ( -_ 
cos @Lp (mn + 1) (wz + 2) 
Kp l¢ 
g ¢ (d Cy — ges) + x 
2 cos ¥Y 


; tC; 
~ Belt - 
sin ¥ 


a C3 
é A 4 


where all notation is referred to Fig. 7. 


( a Ca 
n—1 
p (n + 1) 


A derivation 


(6) 


ol Eq. (6) may be found in the Appendix. 


Comparison of Bending Moments 


It is easier to compare experimental and computed 
bending moments in the plastic range than to compare 
outer fiber stresses; hence, the general procedure of 
this paper has been to take the experimental strain 
corresponding to a given load, determine the necessary 
stress values from a representative stress-strain curve 
and then insert these in the appropriate formula to 
compute the so-called ‘‘theoretical bending moment.” 

The ratio of My,/.W, has been taken as the criterion 
lor comparison of bending moments. 





Mth = é 1 - 27.05 # 


f 


€ : 7a 
Counterclockwise, referred to Figure l. 


Exponential equation: 


theoretical value, 


RESULTS AND DISCUSSION 


In the interest of brevity, detailed information is pre- 
namely, M-4, 
and only the important results for certain of 


sented for a single typical angle beam 
6= 0 
the other specimens tested are given. In particular, re 
sults are shown for beams tested with the plane of load- 
ing at angles of 0°, 30°, 60°, and 76° 07’ to the minor 


principal axis of the cross section. 


Strain Measurements 


Typical load and strain data are presented in Table 
1, where the particular strain gage locations are shown 
in Fig. 6(b). 

In general, strain readings agreed reasonably well for 
With 


some exceptions, the difference was not more than 5 per 


two specimens tested under the same conditions. 


cent in the elastic range and 10 per cent in the plastic 
range. The latter might be partially explained by the 
nonuniformity of the material itself and by the diffi- 
culty in locating SR-4 gages exactly the same on dupli- 
cate specimens. Such discrepancies, however, did not 
greatly effect the computed bending moment because 
even large differences in outer fiber strain did not affect 
the stress values when working on the higher range of 
the stress-strain curve, where the slope is small. 
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Maximum Strains and Neutral Axis Locations 


Maximum compressive strain (e,), maximum tensile 
strain (e,), location of neutral axis intersection with the 
outside of the short leg (x,,), and location of neutral axis 
intersection with the outside of the long leg (y,) were 
all found by a graphical method which assumes a 
linear variation of strain across each face of the beam. 
for beam M-4, @ = 0°—Allustrating 
The procedure was to 


A sample solution 
the method is given in Fig. 8. 
lay out to scale a leg of the beam with proper gage loca- 
tions and for a given load plot the strains for each gage 
perpendicular to the leg. A straight line connecting all 
strains for a particular bending moment was drawn, 
from which e, (or e,) and x, (or y,) could be read 
directly. Table 2 presents these data for beam M-4. 

One would expect the outer fiber strains to become less 
exact when determined as described above as the rota- 
tion and warping increases. However, such errors 
should have little effect on bending moment in the 
plastic range. 

The x, and y, values are approximations only be- 
cause these distances are quite sensitive to small change 
in strain. When a gage was near the neutral axis and 
the strains were relatively small, the effect of errors in 
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the strain was quite pronounced. However, the results 


when considered collectively indicate definite trends. 


Theoretical Versus Experimental Bending Moments 


Theoretical bending moments were computed from 
Eqs. (3) and (6). The per cent differences between 
theoretical and experimental bending moments were 
computed from the expression (\/,,/.17, — 1) X 100 per 
cent. A computation table for theoretical bending 
moment by both Cozzone’s method and the exponential 
method is presented for the typical angle beam M-4 in 
Table 2. The values for all beams are presented graph- 
ically in Figs. 9 and 10. For comparative purposes, 
results for different angles of loading are shown on the 
same figures. The elastic, partially-elastic, and plastic 
ranges are shown on these graphs by indicating on the 
moment scale the approximate values of the expert 
mental moment at which the proportional limit, f,, and 
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the 0.2 per cent offset yield strength, f,, were reached in 
the outer fibers of the beam in question. 

An inspection of Figs. 9 and 10 shows that theoretical 
bending moment was generally within approximately 10 
per cent of the experimental moment in the plastic range 
beyond the yield strength for both Cozzone’s method and 
the exponential method. The differences in the range 
below the yield strength were approximately the same 
for Cozzone’s method but very much greater for the ex- 
ponential method, which is to be expected in view of the 
fact that the exponential relationship between stress and 
strain should not be used below the yield strength. In 
certain cases where the twist was large, the differences 
were slightly greater. 

Beam M-5 which was tested at 60° gave results con- 
sistent with other specimens. It is seen from Fig. | that 
the neutral axis makes an angle of less than 2° with the 


principal axis. There was little rotation of the cross 
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section. It was then decided to test beam N-7 in such a 
position that the neutral axis coincided with the principal 
axis (9 = 76° 07’). The results showed negligible rota 
tion of the cross section although the rotation and trans 
lation of the neutral axis were consistent with other 
specimens. However, the difference between the com 
puted and experimental bending moment appears out 
of line with the results for the G0° beam, M-5. 


Change in Position of the Neutral Axis 


Figs. 11 and 12 show the rotation and translation of 


the neutral axis for the various beams. The values of 


Aa and AT, may be computed from the expressions 

Aa = (90 a tan—(£,/9%,) + 1s 38 (4 
Aa, = 0.261 + 0.160(%,/¥2) — x, (S 
io ae (9 


AT. = Aa, sin (a 


These expressions are derived from simple geometrical 
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relationships referred to in Fig. 1. The angle Aais con- 
sidered positive when the neutral axis rotates clockwise. 
Aa, is positive if it represents an increase in a,. When 
AT, = Aa, sin @ is positive, it indicates that the neutral 
axis is translating toward the compression side of the 
beam. It should be noted that all of these changes are 
with respect to the axes of the cross section. Obviously 
they are intimately related to the cross-section rotation 
as will be discussed in the next section. 

In general, rotations were all counterclockwise and 
were of the order of magnitude of 3° to 4° for all beams. 
The translations were all positive and of an order of 
magnitude of 0.01 in. except for angle beam N-7 for 
which it was approximately 0.02 in. 

The same information is presented in another manner 
in Figs. 13 and 14 by plotting the values of x, and y, 
against applied moment. Since x, and y, were obtained 
by graphical methods, the values are approximate only. 
Curves have been indicated to show general trends. 
When the curves move in the same direction, translation 
is indicated; when they move in opposite directions and 
cross over, rotation is indicated. 
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Behavior of Resisting Couple at Mid-Span as Cross Section 

Rotates 

The torsion couple which causes the beam to rotate 
must be equal to zero at the mid-span because of sym- 
metry. Hence, the only couple acting at this section 
is the conventional bending moment couple. Also, 
equilibrium conditions demand that it be equal in mag- 
nitude and opposite in direction to the external couple. 
It must also lie in a plane which is parallel to that of the 
latter. 

The location of the resisting couple for beam M-4 is 
shown in Fig. 15 for three stages of loading. The re- 
sultant tensile and compressive forces are F, and F,, 
respectively, and are located approximately to scale. 
The forces are evaluated and located by graphical inte- 
gration. The difference in magnitudes was as follows: 


External Resisting 

Moment F, F Arm Couple 

(in.-kips) (kips)  (kips) (in.)  (in.-kips) 
(a) 6.23 5.0 4.7 1.24 5 ome 
(b) 11.90 10.3 10.6 1.00 10.45 
(c) 13.08 13.7 13.2 0.86 II .30 


It is seen that the percentage difference remains es- 
sentially constant with increasing load. Some error is 
introduced by the graphical solution, and also by the 
fact that the true external moment couple does not re- 
main in a horizontal plane but rotates slightly as the 
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BENDING OF ANGLE BEAMS 


beam deflects vertically. Hence, its magnitude be- 
comes increasingly larger than PL and the rotation is 
counterclockwise in Fig. 15. It will be noted that the 
resisting couple rotates in the same direction. 

The sketches of Fig. 15 also show that the neutral 
axis rotates considerably less with respect to the cross 
section than the latter rotates with respect to its original 
position. This difference depends on the shape of the 
cross section and on the degree of plasticity. In this 
case, when the left-hand tip of the tension leg becomes 
stressed beyond the proportional limit, a ‘‘softening”’ 
results and the resultant, /,, moves toward the neutral 
axis and away from the “‘soft’’ area. When the stress 
at the outer compressive corner moves into the plastic 
range, a similar shift is made by F,. Since the re- 
sisting couple must remain parallel to the plane of the 
external couple, the neutral axis rotation must be ad- 
justed with the cross-section rotation so that the stress 
distribution satisfies equilibrium conditions 


Deviation from Pure Bending 


As mentioned earlier, the bending apparatus used 
only roughly approximated pure bending. Actually, an 
axial stress having a maximum value of approximately 
1000 Ibs. per sq. in. existed. Quite often in practice 
combined loadings of this type are encountered; to this 
extent, the investigation of this paper more nearly 
approaches practical applications. In contrast, a 
rigorous method of analysis for the loading used is 
considerably more difficult than for pure bending. 
However, it has been shown that an analysis which 
assumes pure bending gave fairly consistent results 
nonetheless. Since rotation was neglected, the agree- 
ment resulted, to a certain extent, from the fact that 
errors were partially compensating and that moment 
was not particularly sensitive to strain variations in the 
plastic range. 

As a matter of interest, strain gages were placed on 
beam M-5 at the stations A, B, C, shown in Fig. 6(a 
The detailed results are not shown, however the 
moment arm was corrected for the different degrees of 
curvature existing at stations B and C. Inspection of 
the data shows that strains were essentially constant at 
the three stations. This is partially accounted for by 
the fact that while the resultant moment couple in- 
creased from A to C, the bending couple about the neu- 
tral axis remained approximately constant after the ten 
sion component was subtracted. However, the x, and 
y, values show that the neutral axis location was dif- 
ferent at each station as might be expected. The per 
cent difference in bending moment for the three sta- 
tions are shown in Fig. 9(b); they are in close agree- 


ment. 
CONCLUSIONS 


Since rotation of the beam cross section and of the 
neutral axis were neglected in the analysis, any conclu- 
sion must be confined to the particular beam shape used 
in the experiments. If the analysis is extended to 
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shapes having essentially different proportions from the 
one reported in this investigation, the differences be 
tween a theoretically computed bending moment and 
the actual moment might easily be several times as large 
as these reported.* This statement also applies if 
another type of loading is used. The method should be 
extended with caution to materials having essentially 
different stress-strain characteristics. 

The following conclusions can be drawn from the 
results of this particular investigation: 

(1) A modification of Cozzone’s method of plastic 
bending analysis® gave approximate but reasonable re 
sults. For the particular material and configuration in 
vestigated, the correlation between computed and ex 
perimental bending moments in the semiplastic and 
plastic ranges was within 10 per cent. 

(2) An exponential relationship between stress and 
strain also gave reasonable results. For the particular 
material investigated, the correlation between com- 
puted and experimental bending moments in the plastic 
range beyond the yield strength was within 10 per cent. 
The method is not valid for the range below the yield 
strength. Simplification of the method appears to be 
possible but requires further study. 

(3) Rotation of the cross section was the primary fac- 
tor limiting the ultimate strength of the beams. This 
rotation cannot be included in a simplified analysis such 


as was used in this investigation. 
APPENDIX 


Derivation of the Exponential Equation for Theoretical 

Bending Moment of an Angle Beam+ 

A number of equations relating stress to strain beyond 
the proportional limit have been suggested. Several 
of these were investigated by the authors with the view 
of developing a suitable flexure formula for an angle 
cross section. Many of the expressions lend themselves 
quite readily to a simple formula for a rectangular beam 
with loads in a principal plane but become quite compli 
cated for more general solutions. 

It was noticed that when the product fe was plotted 
against e for a number of materials, the outer portion of 


the curve approximated the form 
fe = ae 


It was noticed further that this expression could be 
modified by conventional methods” so that it also 
roughly approximated the straight line portion of the 
stress-strain curve by writing it in the form 

° » n—1 - 

fe =ae"—K or f = ae — K/e Al) 
It is evident that this curve does not pass through the 
origin when e = 0; however, A can be adjusted so that 


* Some results for rectangular and I-section beams are given in 
another paper. !? 
+ The exponential equation has been derived for beams of 


rectangular and I-section in another paper 
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it passes quite close to the origin. Using this equation as 
a basis, approximate expressions can be derived for an 
angle beam. 

If no warping of cross sections occurs (all notation re- 
ferred to Fig. 7) 
é = (r/p) and f = (ar"—'/p"—') — (Kp/r) (A2) 


Considering the area above the neutral axis—-t.e., the 


compression area: 
Outer triangle: 


M, = S fredr, where 


F actt" ( l =e) 
= / = 
ao (n + 1)(m + 2) 2 


Inner triangle: 1/7, = 
ge,”" ; l K pc, 
g n—1 - (A4) 
p (mn + 1)(2 + 2) 2 


M,=M,-—M;, 


and 


~ > (dey"** — ge,"*) — 
p (n + 1)(m7 + 2) 

Kp 
- (dc, — gc,) (Ad) 


Considering the area below the neutral axis—1.e., the 


tension area: 


h — leg: 
M, = JS fredr (A6) 


n- | 


= t/cos y [(a/p "es" tT Im +1) — Koes] 
b — leg: 

M, = (t/sin y) {(a/p"~')[es"*'/(n + 1)]—Kope,} (A7) 
The component \/,, cos ¢ of the bending moment re- 
sisted by bending about the neutral axis is the sum of 
Eqs. (A5), (A6), and (A7). This leads to the value of 
My given as Eq. (6) in the body of the paper. 


Discussion of the Exponential Method 


Values of m and A for materials similar to the one of 
the present investigation are not difficult to determine, 
by the method of reference 17, when applied to the 
stress-strain curve of the material in question. 

An exploratory investigation aimed at simplifying 
the exponential method would appear to indicate that 
for the portion of the stress-strain curve beyond the 
yield strength, Eq. (6) may be modified by setting 


DECEMBER, 19585 


AK =0. The effect is to increase slightly the value of the 


computed bending moment. Again, such a generality 
applies only to the particular conditions of this report.* 
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Transient Temperature Distributions and 
Thermal Stresses in a Skin-Shear Web 
Configuration at High-Speed Flight for a 
Wide Range of Parameters 


H. SCHUH? 
Royal Institute of Technology, Stockholm, Sweden 


ABSTRACT 


The temperature distribution and thermal stresses are calcu 
lated in a configuration consisting of a skin and identical webs 
equally spaced for which only particular examples have hitherto 
been calculated. The basic differential equations for the tem 
perature distribution are brought into a suitable dimensionless 
form, and it is shown that the variables determining the problem 
can be grouped into two dimensionless parameters provided the 
distance between the webs is not much smaller than about half 
their height, which is the case in most practical applications 
By this reduction of variables a general solution of the problem is 
made practicable. Further, two asymptotic solutions of the 
temperature distributions are shown to exist which are well 
known from certain heat flow problems. For cases with the 
same material in the web and the skin, diagrams are given for the 
thermal stresses as functions of the two parameters mentioned 
above and of an additional parameter which comes in when cal 
culating the stresses. Diagrams are also given of the tempera 
ture distribution, and, with their help, complete stress distribu 
tions can be calculated by a simple formula even if the skin and 


the web are of different materials 


SYMBOLS 


length in Fig. 1(c), see also Eq. (1) 
= constant occurring in Eqs. (25) and (26 


length in Fig. l(c), see also Eq. (2) 


( specific heat, skin-friction coefficient (c, 

d = thickness of skin or web 

g = parameter defined in Eq. (17 

h a/( peCod2) 

H = altitude above sea level 

k = heat conductivity 

length of turbulent boundary layer on aircraft surfacc 
L = half the pitch of the web in Fig. 1(b 

1. = Mach Number 

N = half the height of the web in Fig. 1(b 

p defined in Eq. (30) 

Pr = Prandtl Number 

R parameter defined in Eq. (21) 

Re Reynolds Number of boundary layer based on boundary- 


] 


layer length / 


time 
Tr = temperature 
“u = speed of aircraft 
x = distance from middle of web 
a = coefficient of heat transfer from aircraft surface to air 
8 = coefficient of linear thermal expansion 
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= ratio of specific heat at cénstant pressure and constant 


volume 


6 = thickness of heat insulating layer 

€ = strain; emissivity 

# = dimensionless temperature 

6 = dimensionless mean temperature defined in Eq. (31 
x = diffusivity 

& = dimensionless distance defined in Eqs. (11) and (18 
p = density 

o0 = stress 

7 = dimensionless time 

¢ = constant defined in Eq. (29 


Subscripts 


1 = web 
2 = skin 
a = air 


f = boundary layer, friction 
g = total (at heat-transfer coeflicient 
= initial values 


= heat insulating layer 


m = maximum 

Pp = constant pressure 

r =} occurring in Eq. (34 
s = static temperature 


(1) INTRODUCTION 


ERIOUS PROBLEMS arise in supersonic flight due to 
S the heat produced by friction of the air. One of 
these is thermal stresses that occur through temperature 
differences within the aircraft structure. These stres 
ses may be due to steady temperature differences 
for instance, those produced by the heat of the engine 
or they may be of a transient kind, occurring through 
temperature differences set up by a change in aircraft 
speed or altitude. Because it takes only a short time 
to change the speed of an aircraft in comparison with 
the time required for leveling temperature differences 
in solids by heat conduction, the temperatures of dif 
ferent parts of the aircraft structure may differ appre 
ciably, giving rise to thermal stresses. The present 
paper deals with some aspects of these transient phe 
nomena. Evidently the first step is to determine the 
temperature distribution within the aircraft structure. 
To do this for a whole aircraft is a time-consuming 
task and results will have little general value because 
of the number of parameters involved. General infor 
mation can be obtained if simple models are investi 
gated. As such, an airfoil section is assumed here con- 
sisting of a skin reinforced by shear webs [see Fig. 
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Fic. 1. (a) Wing structure consisting of skin with reinforcing 
webs. (b) Wing structure of section (a) replaced by identical 
webs equally spaced. (c) Simplified model used for analysis 
corresponding to parts of section (b) lying within dotted line 


I(a)]; it is further simplified for the present purpose 
by assuming the webs to be identical and equally spaced 
as shown in Fig. 1(b). This model was _ previously 
treated in references 1, 2, and 3 and a number of solu- 
tions have been obtained for certain particular values 
of relevant parameters. Altogether about nine ex- 
amples have been investigated, however the variation 
of the parameters has been limited, so that general 
conclusions drawn from these examples are rather un- 
certain. For instance, only one value of the heat- 
transfer coefficient has been assumed in all the examples. 
In this paper solutions have been obtained for a con- 
tinuous variation of the parameters within such a range 
as appears of practical importance. In the present 
case the great majority of the work is involved in the 
calculation of the temperature distribution: once this 
is known the thermal stresses can be calculated by a 
simple formula. 


(2) ANALYTICAL TREATMENT OF THE PROBLEM 


It is assumed that the same amount of heat is ex- 
changed at the upper and the lower surface, hence the 
temperature distribution will be symmetrical and only 
half the model need be considered. Further, the model 


is replaced by a bar [see Fig. 1(c)] consisting of two 


parts with different cross sections: one (from x = 0 
to x = a) corresponding to the web and the other 
(from x = atox = a + b) tothe skin. The tempera- 


ture is assumed constant in the direction normal to the 
surface of the skin and the web. In the immediate 
neighborhood of the junction of both this assumption 
does not hold and in order to account for it, if neces- 
sary, the relation of the quantities a and 6 in Fig. 1(c) 
to half the height of the web N and half the pitch of the 
web L is left open. For most purposes it is sufficient 
to assume 


a= N = (dy } ) (] 
and b= L (2 


where d; is the thickness of the web. 


The differential equations for the heat flow in the 


simplified skin-shear web configuration are”: 


oT; ot = K,(O77; Ox") Q < * Sa 3 


and 


OT»,/Ot = k2(O?T2/Ox?) + h(T, — To) 


axsx<act+od (4 


7 and 7», are the temperatures, «x; and x. the thermal 
diffusivities in the shear web and the skin, respectively, 
7, is the boundary-layer friction temperature, » the 
distance from the middle of the web, ¢ the time, and h = 
a@/(pxCeds), where a is the external heat-transfer coef- 
ficient due to convection, p», ¢2, and dy the density, the 
specific heat, and the thickness of the skin respectively. 
The heat transfer to the inside of the wing has been 
neglected since it is usually sufficiently small. The 
heat transfer due to radiation has also been neglected 
and this point will be referred to later. The boundary 
and initial conditions are [see also Fig. 1(c) | 


x = 0; O7;,/O0x = 0 5 


x = a; T; = i fe dik\(OT; Ox) = 


2dokx(OT» Ox) (6 
OT>/Ox = 0 (7 
t = 0; ,=7,=T (8) 


Since the heat-transfer coefficient a changes with 
speed and altitude, @ is in any particular case a fune- 
tion of the time, as is also the boundary-layer friction 
temperature 7. In order to obtain general solutions, 
instantaneous start of movement with constant speed 
and horizontal flight is assumed, so that constant values 
of 7; and a can be used. The same assumption has 
also been made in references 1, 2, and 3. It has been 
shown” ‘ that the temperature distribution and the size 
of the thermal stresses are not much affected by the 
times of acceleration which occur in high-speed flight, 
although the time when the maximum stresses are 
reached may depend on the time history of the aircraft. 
In an acceleration assisted by rockets or afterburners, 
instantaneous start is at any rate a good approximation 
for both maximum stresses and their time history. On 
the other hand top speed can often be maintained only 
for a short time and then the speed decreases; in such 
cases the present solution in its simple form can be 
used only during the time when the top speed is main- 
tained. With guided missiles it should be noted that 
the projectile may hit its target before thermal stresses 
have reached their maximum value. Cases where the 
speed or the altitude of the aircraft changes during the 
major part of the flight cannot be covered by the solu- 
tions of this paper. 

In an aircraft structure, rivets and screws are gener- 
The thermal 


ally used to join the web and the skin. 
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contact in such junctions is not perfect, so that a tem- 
perature jump occurs across the interface between the 
two parts, which may influence the temperature dis- 
tribution in the structure. However, such cases will 
not be treated here and instead perfect thermal contact 
is assumed between different parts of the structure. 

It has already been observed in refernce 2 that in 
any of the particular examples calculated the tempera- 
ture gradient in the skin vanishes with sufficient accu 
racy ata distance from the junction of the skin and the 
web which is much smaller than }. In view of this 
fact, the boundary condition Eq. (7) at x = a + bis 


replaced by 
v7 @; O12/0x = 0 (9 


which is general and does not refer to a particular case. 
After obtaining general solutions, the minimum permis- 
sible size of b = L will be determined for which Eq. (9) 
holds; anticipating a later result, this size of L does not 
limit the applicability of the present solution in most 
cases. One advantage of this modification of the 
boundary condition is that one of the parameters of 
the original problem as determined by Eqs. (3) to (7 
disappears and hence a general solution of the problem 
is made practical. Another is that the same solution 
covers also cases with different materials for the skin 
and the web. 

The following dimensionless temperatures, coordi 


nates, and time are introduced: 


vd; (7, — 11) (7, — 7;)] (10 
O z re. @ 
E vd \ (11) 
v 1 / (iy l,) (12 
aZxZat+b 
& (x — @ a|Nv kK, Ke \ (1:3) 
T {ky a- (14) 


Introducing these quantities in Eqs. (3) and (4) one 


obtains 
Ov; OFT OW’, O(E,)° GZ é.z 7 (15 
OW. Or = [0° /O(é)?] — go OZ Z om (16 
where a parameter 
g? = a*(h ny) = a*(apicy/dok poco) (i7 


has been introduced, which appears to be of general im- 
portance in problems of heat flow in structures whose 
solidity is small. The boundary conditions are re- 


written and read: 


& = O Ov, Of = 0 (18 
s=1 and & = 0: I, = do: 
Ov; Of; = R(OdI. OF) (19) 
f2—> ©; Ovs Of = 0 (20) 
where 
R = 2(d2/di)V Repoce ki pics (21) 


he initial conditions become 
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On inspecting Eqs. (15) to (21) it follows that the 
solution is only determined by two parameters* besides 


the two independent variables and is of the form 


dv = a(t, 7; R, g (23a 
J. = Velho, 7; R, g (23b 


Attention is drawn to two asymptotic solutions, 
which can easily be found. The first asymptotic solu 
tion is obtained when g ~ ©, provided that d. ~ 0. 
If do 0, there would be no external heat transfer as 
inspection of Fig. l(c) shows. For large values of g the 
second term on the right-hand side of Eq. (16) is greater 
than the first term, and the solution of Eq. (16) 1s 


ve exp. g-T) 24 


Consequently J». is independent of & and because of 
the solutions of J 
the thickness ratio 


the boundary conditions Eq. (19 
and J» are independent of R—1.e., 
of web and skin. From Eq. (24) it follows that, for 
g — ©, % becomes zero before J; has appreciably de 
parted from its initial values. Hence this asymptotic 
solution is equal to the temperature distribution in a 
slab of thickness 2.V which is initially (at / Q) ata 
constant temperature different from zero and which is 
kept at a surface temperature equal to zero for all times 
{> 0. The solution of this problem is well known. 
The second asymptotic solution is obtained when 
L—0O; then the skin-shear web configuration becomes a 
solid wing or—-in the simplified model of Fig. 1(b)——a 
slab with thickness 2.V, as in the first asymptotic solu 
tion, but with external heat transfer taking place from 
the surfaces of the slab with a coefficient a to a medium 
whose temperature is 7,. The solution of this problem 
is also known.’ The solid wing can be considered as a 
extreme case of the skin-shear web configuration shown 
in Fig. 1(6); however it does not fit the simplified model 
of Fig. 1(c), since there would be no external heat trans 
fer for b = L > O according to Eq. (4 The parameter 
determining the temperature distribution is g (aa) k 
as follows from Eq. (17) by putting ds a there, but R 
disappears as parameter. This asymptotic solution 
will not be discussed here although it could be a valuable 
supplement to the solutions presented later, since these 
require values of L to be bigger than a certain lower 
limit. 

In references 2 and 3, solutions have been obtained 
in the form of a series expansion with the help of a 
Laplace transformation. However, in any particular 
case the computing work seems to be rather time 
consuming, particularly since the convergence of the 
series is rather poor at small times. This point seems 
to be confirmed when observing that only one numerical 
example has been calculated in reference 3 and 8 ex 
amples in reference 2, although in preparing the latter 

* Parkes? found also by a somewhat different approach that 


two parameters determine the temperature distribution in the 


examples he calculated 





paper a computing staff was available. It has been 
shown in reference 4 that the computational work can 
be considerably reduced by using a suitable numerical 
method. Essentially, that method consists in replacing 
the differential Eqs. (15) and (16) by finite difference 
equations and using a step-by-step procedure. Good 
accuracy can be obtained with this method, even for 
rather large spatial steps. This numerical method be- 
comes the more accurate the smaller the finite differences 
of space and time are, but at the same time the amount 
of numerical work increases; hence a reasonable com- 
promise has to be made. A high degree of accuracy is 
not necessary in the present problem, because certain 
errors are unavoidable in applying to actual aircraft 
structures the results of the simplified model treated 
here. In order to save time a maximum error of about 
+().03 for temperatures and thermal stresses was al- 
lowed for in the present solution. 
The thermal stresses are, according to reference 6, 


o = EiB(7T; — 71 + Aj) 0 42 & 2 I (33) 
o. = bhep(T, — TV. + Ad) 0 2 & < @o (26) 
o1, 02, He, Ey, and 6), B: are thermal stresses, Young’s 


modulus, and the thermal coefficients of linear expansion 
for the web and the skin respectively; A; and A: are 
two constants to be determined by the following condi- 
tions: (a), that the strain «; and ¢ for web and skin is 


the same—1.e., 


€;) = Ap = €& = A>vB> (27) 


and (b), that there are no external forces at the ends 


a i+b 
dy [ 0} dx + 2d» [ oo dx = 0 (28) 
0 /a 


1.¢., 


It is convenient to introduce the following quantities 
g = 2(d.E.b/d, Fa) (29) 


p= (b/a)V k1/ke (30 


1 ca 
| é 

6, = [ 0, di, 0(&) = dv» dé 51 
/70 & J0 


Eqs. (25) and (26) can be rewritten as 


0; = kieiu(T;, = iz 


f A, + o(C2/P1)O0(p) + ¢| 1 — (@2/P1) | = 
v; = (32) 
| l+Q9 
0» = hoB>( (a iil 7’;) 
GA. ( p ) + | Bi Bo )6; i. i 3 2) 
i, — (.3:3 ) 


The stresses can therefore be easily found, if the tem- 
peratures 3, J, 0, and 6 are known. 


(3) RESULTS 


Numerical solutions have been obtained of Eqs. (15) 
and (16) with the boundary and initial conditions Eqs. 
(18) to (22) for values of g = 0.923, 1.85, 3.70, 7.38 and 
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for values of R = 7.5, 2.0, and 1.0. The result of caley- 
lations for R = 2.0 are presented in Figs. 2-S where 
each figure refers to a certain time 7; tiie temperatures 
3; and J» are plotted against g with dimensionless dis 
tances £ (measured from the middle of the web) and 
& (measured from the junction of the web and skin 
as parameters; the figures contain also 6;, the mean 
temperature of the web, and 6,(&,), and the mean tem 
perature of J) between O and &,. The value of &, in 
the diagrams has been chosen so that within the ac 


curacy of present solutions 


rer 


O2(E2) = (ko,/E2)O2(fo,) + [1 — | £) Io « 34 


for any values of & > &2,. 

The temperatures for large values of g were obtained 
from the first asymptotic solution discussed previously. 
From these and similar figures for R = 7.5 and 1.0, 
which owing to lack of space are not reproduced here, 
an estimate can be made of the minimum size of ZL, for 
which the present solutions are valid. The tempera- 
tures d. differ at & = 1 only a little or in some figures not 
at all from those at &— ©. Hence the exact boundary 
condition for any particular case Eq. (7) can be replaced 
by the simplified boundary condition Eq. (9) if & > 1; 
this condition can be relaxed to & > 1/2 for reasons 
given later. Hence it follows that 


pb = (b/a)V k/Ko => 1/2 (35 


and thus the smallest size of b = L is determined for 
which the present solutions are valid. 

With the help of these figures the stress distribution 
can be found for any time 7. If ¢ is not too small, the 
maximum stress occurs in the middle of the web, hence 
it suffices to calculate stresses for that point and plot 
them against time in order to find the maximum. These 
curves are rather flat near their maximums, which can 
thus be determined with sufficient accuracy from those 
times to which Figs. 2-S refer. 

If the web and the skin are of the same material, 
further simplifications are possible. Because the posi 
tion of the maximums of the stress-time curves does not 
vary much with ¢ and because these curves are rather 
flat, the time 7,, at which these maximums occur can be 
assumed approximately the same for all values of the 
parameter ¢ ranging from g¢ = © down to the lowest 
corresponding to the smallest possible size of L. Values 
of dimensionless time 7,, vs. parameter g are shown in 
Fig. 9 for R = 7.5, 2.0, and 1.0. It is important to 
know these times 7,, for judging whether the present 
solution for constant speed after the start can be ap- 
plied to cases where the operational diagrams apprect- 
ably differ, for instance, when the top speed can be 
maintained only for a short time, as occurs with many 
high-speed aircraft. 

Again for cases with the same material of the web and 
the skin, diagrams can be given for the maximum 
stresses occurring in the middle of the web as a function 
of the parameter g for various values of ¢ and three 
values of R. The times taken to reach maximums of 
stresses being independent of the parameter ¢ greatly 
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The dimensionless times in the figures are 0.05, 0.1, 0.2, 
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transfer coefficient due to convection in air for a turbulent bound 
ary layer as a function of altitude above sea level /7, Mach Num 
ber J and Reynolds Number Re of boundary layer based on 
boundary layer length /. fis a function of the Prandtl Number, 
its value for air not vet accurately known. A value of f = 1.10 
is tentatively proposed in this paper 


simplifies preparation of these diagrams. ‘The results 
can be seen in Figs. 10-12, which were obtained for R = 
2.0 from Figs. 2-8 and for R = 7.5 and 1.0 from calcu- 
lated temperature distributions not published in this 
The influence of R on the stresses is small, par- 
For values of g > 6 the 


paper. 
ticularly at bigger values of g. 
maximum stresses are independent of g—1.e., inde- 
pendent of the heat-transfer coefficient a, height of web 
a, and thickness of skin d,. For g < 3 they are about 
proportional to the square root of the heat-transfer co- 
efficient, to the height of the web and inversely pro- 
portional to the square root of the thickness of the skin. 
The 
parameter ¢ can, however, not be extrapolated beyond 


decrease of maximum stresses with decreasing 
those values given in the Figs. 10-12, since for suf- 
ficiently small values of g the maximum stress occurs 
in the skin. This is illustrated by Fig. 13, in which the 
maximum stresses in the middle of the skin and in the 
junction of the skin and the web are plotted vs. ¢ for 
the parameter g > ©; below a certain value of ¢ the 
maximum stresses in the skin become greater than those 
in the web. Besides, the times taken to reach the 
maximum stresses are appreciably different for both the 
skin and the web. 


A 


T 


URE DISTRIBUTIONS 835 
Previously the minimum size of L was estimated, 
which is permissible when using the standard solutions 
derived in this paper. This estimate has been checked 
by calculations for R = 2.0, g = 
values of L/N = b/a = 0.25, 0.50, 
exact boundary condition at x a + b, Eq. (7 
used instead of the simplified boundary condition Eq. 


1.70 and 3.55, and 
1.25 in which the 


, Was 


(9). These solutions were compared with the standard 
solutions and little difference two 
found for values of LN larger than 0.50 and even for 
L/N = 0.25 a value not altogether wrong is obtained. 


between the was 


In view of these results a value of | 2 was chosen for the 
lower limit of p in Eq. (35). 

For predicting temperature distributions and thermal 
stresses, the variation of the heat-transfer coefficient a 
with altitude above sea level /] and Mach Number .\/ 
has to be known, and this information may be obtained 
in a convenient form from Fig. 15(a) which refers to a 
turbulent The 
efficient for air was assumed as follows’* 


boundary | layer. heat-transfer co 


a = f(Pr) palpatt(c,/2 56 


where p, and C,,_ are the density and specific heat of air, 
u the aircraft speed, and c, the coefficient of friction for 
a turbulent boundary layer which has been taken from 
The factor /(Pr) allows for the influence 
1 for 


reference S. 
of the Prandtl Number Pr on heat transfer; 
Pr = 1. According to the well-known Colburn for 
mula’ {(Pr) = (Pr)-*’*: 
ficient experiments available to decide whether 


however there are not suf 
this 
formula gives accurate values for air (Pr 0.72) as a 
seem to indicate 


0 


Some recent experiments! 
1.10, which value is about 10 per cent 


medium. 
that f(0.72) = 


* Except for a region near the leading edge the pressure gradi 
ent on airfoils is sufficiently small to be neglected in determining 
the heat transfer coefficient. Further it is assumed that the 
boundary layer remains attached and is undisturbed by shock 
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Fic. 15(b). Diagram for calculating Reynolds Number as 


used in Fig. 15(a) as a function of altitude 7 above sea level, 
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TABLE | 
Ratio of Maximum Thermal Stresses o;»1/o¢1m with and Without 
a Heat-Insulating Layer of Thickness 6 and of Heat Conductivity 
k,. Characteristic Parameters: g = 3.0 (bare wing), R = 2.0, 
and ¢ = 10 


(a) Flight conditions: HT = 50,000 ft., 17 = 2, and Re = 1 X 105 


F B.t.u 

6 (in.) ky = 0.134 = 0.0336 0.0134 aa 
. ft.hr.°F 

0.002 1.00 0.99 0.96 

0.01 0.98 0.91 0.82 

0.02 0.96 0.85 0.71 

0.1 0.82 0.56 0.39 

(b) Flight conditions: H = 15,000 ft., M7 = 2, and Re = 1 X 105 

6 (in.) k, = 0.1384 0.0336 0.0134 

0.002 0.98 0.93 0.85 

0.01 0.91 0.74 0.55 

0.02 0.85 0.60 0.438 

0.1 0.55 0.31 0.21 


lower than that of Colburn. The heat-transfer co- 
efficient in Fig. 15(a) decreases appreciably with increas- 
ing altitude, but does not change much between J = 3 
and J = 5 at constant altitude. Fig. 15(b) facilitates 
the estimate of Reynolds Number, where / means the 
length of the turbulent boundary layer on an aircraft 
surface. 

For completeness’ sake, the temperature differences 
occurring in the formula for the stresses are given. 


T,-— T, = T31 + [(y — 1)/2])(Pr)'? a7} 
Toil 1 + [Cy — 1)/2](Pr)'7(.M,)?} (37) 


provided the aircraft was initially in equilibrium with 
the surrounding air. y is the ratio of the specific heats 
at constant pressure and constant volume. 


and s refer to boundary-layer friction temperature, 


Indices f, 7 


initial temperature, and static temperature, respectively. 
Again, a turbulent boundary layer is assumed. 

In order to estimate the effect of a heat-insulating 
layer mounted on the outside of the aircraft skin, the 
heat-transfer coefficient a is to be replaced by a, in 
calculating the value of the parameter g, where a, is 


l/ag = (1/a) + (6/k,) (38) 


6 is the thickness of the heat-insulating layer whose 
thermal conductivity is k;. The thermal capacity of 
the heat insulating layer has been neglected, which is 
correct for thin, but not for thick layers. The values of 
k, for three typical heat-insulating materials are: (a 

plastics with asbestos fibers, 0.1534; (b) best inorganic 
heat-insulating material, 0.0336; (c) air at +32°F., 
0.0134 (B.t.u.)/ft. hr. °F.). Air has been given since 
it is the best heat insulator and the effective component 
of porous materials. The influence of the thickness of 
layers of these heat-insulating materials on maximum 
stresses can be estimated from Table 1(a), which refers 
to a typical example of flight at 50,000 ft. altitude and 
Re = 1 X 10%. It can be seen that appreciable 
effects can be obtained at altitudes of 50,000 ft. with 
layers of about 0.05-in. thickness and the best of heat- 
insulating materials. However, at lower altitudes heat- 
transfer coefficients are higher and the thicknesses of the 
insulating layer can be less for the same stress reduction 


[see Table 1(b).] 
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TABLE 2 


Maximum Thermal Stresses for the Examples Calculated in Ref 
erences 2 and 3 Compared with the Corresponding Results of 
Present Solution 


Tin 
FE, Bit T; - 7 
Source other 
(refer-  Ex- Br» auth present 
ence) ample g ¢Y By R ors author 
l l 3.56 20 l 7.5 0.79 0.80 
y 2 3.56 5.33 1 2.0 0.72 0.73 
2 3 6.89 5.3: I 2.0 0.85 O.85 
2 4 6.89 2.785 |] 2.0 0.76 0.76 
2 a) 6.88 5.33 ] 2.0 0.85 0.85 
2 6* 3.97 5.33 l 2.0 0.73 0.76 
2 ‘i 16.00 0.5 1.414 0.76 0.74 
2 8 3.56 5.33 ] 2.0 0.72 0.73 
3 | 3.56 20 l 7. 0.79 0.80 


* In the original paper obtained by interpolation 


Lastly, the maximum stresses of the nine examples 
calculated in references 2 and 3 are compared in Table 
2 with those obtained by the present analysis. Com- 
parison is very good, particularly when noting that the 
biggest deviations occur in cases where the solutions of 
reference 2 have been obtained by extrapolation and it 
is therefore possible that they are less accurate. 

It should be noted that in references 1, 2, and 3 a 
heat-transfer coefficient of 90 (B.t.u.) (ft-"hr.°F.) has 
been assumed in all cases. In references | and 2 this 
value was stated to correspond to an altitude of 50,000 ft 
and a Mach Number of 3.1. 
The heat-transfer coefficient at the above flight 
condition (see Fig. 15) would be 21 (B.t.u. 
(ft.*hr.°F.).{— Alternatively, the heat-transfer coeffici 
ent used in references |, 2, and 3 corresponds to an alti 
tude of 15,000 ft. and the same Mach Number of 3.1. 
Consequently the stresses found in references 1, 2, and 


However, this is errone 


ous. 


3 would be reduced by between 15 to 40 per cent if the 
correct heat-transfer coefficient for an altitude of 50,000 
ft. were assumed. For the same reasons the prospects 
of reducing stresses by insulating layers at altitudes of 
about 50,000 ft. or above are less favorable than ap 
peared in reference 2. 

The influence of radiation can be easily considered for 
the asymptotic solution when g — ©; then the tem 
perature of the skin reaches equilibrium before an ap 
preciable heat flow occurs into the web. The skin can 
therefore be treated separately and its equilibrium tem 
perature 7), due to convection and radiation can be 
calculated from the following equation: 

eCs(T,,' T,. = a(T, — 7\,) (39 
Cs is the constant of radiation for a black body, ¢ is the 
emissivity and 7’, the static air temperature. The radia 
tion from the body is assumed the same as from a “‘gray 


(Continued on page 866) 


7 In reference 2 Hoff's value for the heat-transfer coefficient 
was checked by estimating it from the drag of aircrafts by using 
Reynolds’ analogy. However only the frictional part of the drag 
can be used for that purpose and that part is much less than the 
total drag in modern high-speed aircraft 

t Assuming a boundary-layer length of about 10 ft. 


4 
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Computation of Vibration Modes and 
Frequencies on SEAC 


WILLIAM F. CAHILL* ann SAMUEL LEVY? 
National Bureau of Standards 


SUMMARY 


A description is given of several codes developed for use in 
he Standards Eastern Automatic Computer for computing the 
free vibration modes and frequencies of aircraft structures. The 
codes are general in nature and apply without modification to any 
aircraft structure whose influence coefficients, mass distribution, 
ind geometry are specified. The codes are applicable when no 
more than 63 mass points are needed to define the displ icemenit 
half-wing. In an example with seven mass points along 
ixis and 24 mass points for each half-wing, the com 


of the 
the fusel ige 
puting time for the first six symmetric and the first six antisym 
modes was 6 hours 


metric 


INTRODUCTION 


A KNOWLEDGE of the vibration modes and fre 
quencies of aircraft is of primary importance in 
determining their flutter speeds and in making dynamic 
analyses of the response of the structure to landing im 
pact and other shock loads. A great deal of effort has 
therefore been expended in developing convenient 
methods for carrying out the tedious computations in 
volved. 

The method of computing vibration mode shapes and 
natural frequencies by matrix iteration is a systematic 
procedure which has been gradually developed from 
the initial work by W. J. R. Collar. 


rhis procedure converges to the lowest frequency mode. 


Dunean and A. 


To obtain the higher frequency modes, use has been 
made of the orthogonality relationship to (a) reduc 
the size of the matrix successively by one and thus ob 
tain convergence to the higher modes, b) “‘sweep 
lower modes from the trial modes and thus obtain con 
vergence to higher modes,* * or (c) modify the original 
matrix without changing its size so that the lower modes 
are successively removed.° 

rhe codes to be described in this paper are based on 
the matrix iteration scheme together with a variation 
of the matrix modification scheme given by Kincaid in 
The rigid body modes are eliminated as 
This over 


reference 5. 
described in Appendix III of reference 6. 
all choice is particularly well suited to a high-speed 
computer such as SEAC because it results in highly 
repetitive computations where speed and accuracy are 
of primary importance. 

The Standards Eastern Automatic Computer, SEAC,’ 
is an electronic digital computer containing a 512-word 


Received October 7, 1954 

" Mathematician 

t Formerly, Chief of Aircraft Structures Group, National Bu 
Now, Manager, Stress Analysis Engineering, 
Special Defense Projects Division, General Electric Co 


reau of Standards 


mercury delay line memory and a 512-word electro 


static memory (not used in this problem). It also 
has internal auxiliary tape storage units which make 
it possible to store the many numbers required in the 
present problem. It computes at an average rate of 


about 500 operations per second. 


METHOD 


The method used in setting up the code follows that 
given in reference 6, Appendix III, but treats the sym 
metric and antisymmetric modes separately so that 
This 


limits the applicability of the code to symmetrical 


only a half-wing need be considered at a time. 


structures. 


Consider a frame of reference coordinates x, y, 2 
fixed in space and oriented with respect to the airplane 
as shown in Fig. 1. Let the mass of the wing be repr« 
sented by a moderate number of lumped masses. Let 
Let the origin 


Vn», Vy, Z, be the coordinates of mass n. 


be at mass 1. The influence coefficient 6, between 
two lumped mass points m7 and » on the airplane is her« 
defined as the vertical displacement at m due to a ver 
tical unit load at m when point | is clamped and when, 


for determining symmetrical modes, the slope dy dx, 








Fiqure |. - Coordinate 


system 
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normal to the fuselage axis, is taken as zero all along 
the fuselage axis; while, for determining antisymmet- 
rical modes, the deflection y is taken as zero all along 
the fuselage axis. (As the relative fuselage rigidity in- 
creases, the symmetric and antisymmetric influence 
coefficients tend to equality.) Let the motion of the 
clamped point 1 be described by X, 6, a; where X is the 
upward displacement of point 1, @ the roll from x to 
y about the z axis, and a@ the pitch from z to y about 


the xv axis. 


Antisymmetric Modes 


Consider the fundamental antisymmetric mode with 


More accurate values, y;’, yo’, .. . Vn’, 


a frequency w. f 
of the displacements parallel to the y axis may be com- 
puted by iteration from assumed values y;, yo, .. . Y, 
in terms of the antisymmetric influence coefficients 
and the inertia forces due to vibratory motion of 
masses M,, Mo,...M, by: 


Nn 


a 
-(y' — 0m) = D Myd, = DS B,,y, 


w r=1 I 

| n n 

(ye’ — Ox2) = yi M,v,62, = p ® Bo,“ y 
Ww 7 1 ; l 


= > B,,(y 
’ l 


had 

(¥n’ — Ox.) = DY My, 
Ww r=! 
The B coefficients are frequently referred to as terms 
The superscript (1) indicates 
To determine 


in the dynamic matrix. 
they apply for the fundamental mode. 
# we form the sum 


to “a 
- } x,M,(y,’ — 0x.) = > x,M, >> By, (2 


Ww l 


To obtain convergence to the next higher mode, it is necessary to modify the right-hand side of Eqs. (1 


ing reference 5, this is done as follows: 


(ny! — 6x1) = 2. B,,“y 


Py 


W 


Eqs. (6) make use of the orthogonality principle to eliminate from any trial mode y;, ys, . . . 


AERONAUTIC 


| “ 
—(vn’ — Ox,) = 2. B,,“Yy, — 
? | 


> 
1 


Ww] 


a 
=i 


WI 
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and make use of the fact that if the airplane is not 


being accelerated as a rigid body 


n 


> x,M,y,’ = 0 3 


1 
») 


giving on substitution in Eq. ( 


n 


> x.M, >, Bs 
| j l 


6 
fone aoe n f 
Ww ) 
> x2 MV 
1 
a hs , . a ? / 4 / 
The more accurate values y;'/w", yo’ /w”, <9 Va /@ ate 


obtained from Eq. (1) using the value of 6/w? given by 
Eq. (4 


Vi 6 a 

=—x+ > By,“y 
Ww Ww } l 
Vo’ 0 = 

= \ > Bo, 
Ww Ww } l 

5) 

Tag 0 a 

= Xn t Dd Ba 
uw we r=1 


The ratio y, (y,’/w”) gives an approximate value of the 
square of the fundamental frequency «,° while the ratios 
(y1' w*) (y,-' w*), ete., give deflections in the funda 
mental antisymmetric mode normalized with respect 
to a preselected station r. These displacements de 
scribe the mode shape more accurately than the origi 
nally assumed values. By repeating this cycle of com- 
putations, in each case starting the new cycle with the 
more accurate values just obtained, convergence to 
the fundamental mode is obtained. 


Follow- 


My Py yy, 


> ILy 


= > B,,y, 
> Myr 77! 


, ¥, the portion which is 


fundamental mode, and thus bring convergence to the next higher mode of the iteration process represented by Eqs. 


(1) to (5) with B™ replacing B“. 
To obtain the third and higher modes, the successive dynamic matrix terms are obtained from 











To ¢ 


and 1 


Subs 








not 


ire 
by 
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My y_" 
B,,{'*+) = B,, - am . 
w? > M,[y,“]? j 
l 
and the iteration process represented by Eqs. (1) to (5) 1s repeated. 


Symmetric Modes 
Consider the fundamental symmetric mode with a frequency w. More accurate values y,", yo’ 
placements parallel to the y axis may be computed by iteration from assumed values yi, vo, ... , y, 


symmetric influence coefficients and the inertia forces by 


l 
(Vy _— \ a » M_y,6; B, \ 
w 1 
l “a a 
(yo’ — A az.) = >, M,yb. > Bz, 
Ww” ] , 1 
XQ 
l a oe 
(y,’-A—a > My > Bay 
a ] r 1 
To determine \ and a@ we form the sums 
os a Sl 
yo ba M LF, —_ \ — az = » M Z 2 B ede 
Ww 1 l y=] 
{) 
Le a 2 
YS My,’ — dX — at > Md By 
Ww l l l 


and make use of the fact that if the airplane is not being accelerated as a rigid body, 












































>, sAf9.' 0; YS May,’ 0 10 
Substituting Eq. (10) in Eq. (9) and solving for X and @ w” gives 
TABLE | 
Mass, x» and s Station Coordinates, and Antisymmetric Mode Shapes and Frequencies 
— a 
Station Mass x z Mode 1 Mode 2 Mode 3 Mode Mode 5 Mode 6 
lb sec*/in, in. in. 
1 | 0.01548656 49.497474 ) 0.1053281 | -0.2994552 0.2696731 1.233706 1.8963363 0.7571406 
2 | 0.02353352 91.923881 ) 0.2487737 | -0.6269135 0.6068175 1.110532) 3.0521377 1.5830190 
3 |0.03005745 | 134.350288 ) 0.3784546 | -0.7872769 0.7848211 | -0.4930628 1.510221 0.8791356 
4 | 0.06446276 | 176.776695 0 0.508352 | -0.6894540 0.8096368 | -2.2153764 | -1.251740k | -0.8962385 
5 |0.05782825 280. 862813 0 1.0000000 1.0000000 1.0000000 1.0000000 0.0002), 80 0.213383 
6 |0.03844221 49497474 52.325902 0.030368 | -0.1327213 0.1280213 0.701223 0.807622z 0.2196500 
7 |.0.05315167 91.923861 0.064833 | -0.2)19023 0.135307 1.0086363 0.8656522 | 0 | 
8 0.06155684 134350288 0.1047606 -0. 3083052 -0.0186641 0.5459477 0.1979509 0. 
9 |0.07419415 | 176.776695 0.1469578 | -0.2855751 | -0.4161736 0.2403583 | -0.6053620 | 0. 
10 |0.11244671 | 250.457221 52.325902 0.2629954 0.0525995 | -2.1532393 | -1.2127619 1.0000000 | 
11 | 0.30318519 9.497474 | 108.89444y | -0.0048429 | -0.0521661 0.0789027 0.434089) | * 0.361473 | 
12 | 0.26769497 91.9238681 -0.0119559 | -0.0864503 0.0581000 0.7980655 0.2799138 
13 | 0.27562399 | 134.350288 -0.0182794 | -0.102)706 | -0.0840490 1.021047 | -0.2527450 
ly |0.15552180 | 176.776695 -0.022572) | -0.0847514 | -0.3719861 0.9786777 | -1.0727818 
15 |0.06074808 | 217.930309 | 108.894444 | -0.0261044 | -0.0171589 | -0.7701349 0.730786 | -1.8185453 
16 |0.90096093 49497474 171.119841 -0.0258761 | -0.006673 0.0651073 0.1182361 0.151290 
17 | 0.68739492 91.923881 -0.0535989 | -0.0030001 0.1065604 0.2307441 0.1756180 
18 |0.5546400) | 134.350288 -0.0838845 0.0107528 0.099668), 0.3266008 | -0.0380381 
19 |0.13027820 | 176.776695 | 171.119641 | -0.1160221 0.0388201 0.0539281 0.2545818 | -0.3468054 0. } 
20 |1.44840038 9.497474 | 233.345236 | -0.0373266 0.0173477 0.0689152 | -0.1383700 0.0754152 0. 
21 |0,.92140026 91.923881 { -0.0737608 0.0407854 0.1307565 | -0.3362522 0.1180312 ; 
22 |0,31614511 | 134.350288 | 233.345238 | -0.1130376 0.071066) 0.1777964 | -0.5746732 0.0767416 0.1038077 
23 |1.18136199 9.497474 | 295.570634 | -0.0437203 0.0305001 0.0693467 | -0.2954392 0.023680 0.1845587 
24 |0.60981791 91.923881 | 295.570634 | -0.0846868 0.064032 0.130668 | -0.6651949 0.0080449 0.4989885 
w (rad/sec) 41.971 110.637 169.813 238.643 351.901 399.850 
— 
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; 0) (E we) (E05) 


It should be noted that in the special case of a thin straight line beam, where ¢ is zero at all stations, Eqs. (11) and 


12) are replaced by 


> MY By 


i ~ 
; il S: z 0 
a = 
MV 
] 
a a if 0 Io! the dynamic matrix and stores it on tape for use in 
, if 2 So edie (1g Ke : : : 
a computing the mode of next higher frequency. 
[he more accurate values yy" w*, yo w*, ..., Yn /@ Formation of Dynamic Matrix and Constants for Removing 
are obtained from Eq. (S) using the values of \ w* and Rigid Body Motions 
aw given by Eqs. (11) and (12). Rewriting Eq. (8), rhe first code reads in, converts into binary form, 
vy! \ - and stores the mass constants, .1/;, J/o,...,. 17,,, and the 
ae 5 - 6 ; . 
=— + —& + 2 Beh geometric constants, 
Ww Ww eo) l 
Vo" r a al hy, the, «ss 5 0 ee Rio ae 14 
— = t oo + DY Boy s 
a w a ' (for antisymmetric modes 











¢ 
Vv,’ nr a a 
=— +-— os. 2, Bar ¥ 
aw Ww WwW , ] 
The ratio y, (y,’ w*) gives an approximate value of . 
the square of the fundamental symmetric frequency eee | 
w;” while the ratios (¥;'/w") /(y,’/w"), ete. give deflec e 
tions in the fundamental symmetric mode normalized bs 
with respect to a preselected station » which are more : aa | 
accurate than the originally assumed values. Iter 3 
ation of this cycle of computations converges to the 2 
fundamental mode. oe | 

To obtain convergence to the second mode and . 
higher modes, B,,‘") in Eq. (S) 1s replaced successively S i 
by B,,"°’, B,,, ete., given by Eq. (7), and the iteration 7 
process represented by Eqs. (S) to (15) is repeated. S 

MACHINE COMPUTATION ° 
» 

The problem is broken down into three parts for pur 3 
poses of programming due to the limitation of memory 2 [ | 
space in SEAC. The function of the first code is to ~ 
form the dynamic matrix [B,," coefficients] and store \ ls — 





it on tape and also to determine the constants used in ? 90 ssc 50 200 


eliminating the rigid body motions and print them for 
use with the second code. The second code computes 


the frequency and mode by the iteration scheme pre- 
Figure 2- Locations of stations on one wing 


sented in the previous section. The third code modifies 
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M1, Me, ~ 21, 22, 2 14’ : — ; 
F 1/D (> vate), one or more u's + 0 S 
for symmetric modes r=1 
1 ." 
It sums the elements of the above vectors as a check — (x VW ), all u's U 
that they have been read in correctly and stores them ; 
internally rhe code then computes the products G 1D (x / ), one or more a’s * 0 19 
; : Vor [1 13 
Z Uy, U2, cope V/ 141, .biolto, o* M it 1+) (), all u's () 
nd the co ) E, F ’ where. for antisy _ 
ind the constants D, E, F, and G, where, for antisym Che numbers obtained in Eqs. (17) to (19) are printed 
nd metre modes, out for use by the second code The first code then 
ilif . 
\- reads in the influence coefficient matrix a row at a time, 
) 7 1G ; P ; ; 
l ae converts 1t into binary form, computes the rows of 
corresponding terms B in the dynamic matrix, and 
I iz stores them together with their row sums on magnetic 
F 0 IX tape for use by the second code. A running sum of 
the terms in the influence coefficient matrix is also 
G 1/D 19 


computed to be printed at the end as a check on “‘read 


ind, for symmetric modes, in. 


a D (x V/ ) (x . “) [s , ) 16’ Iteration 
, The vi 


second code chooses as a starting vector, 

) > - by.4 Tere. 1. It then reads the dynamic 

ng E L/D) Et 2 Dh one or more u's #0 (17 oe io ' 
matrix a row at a time from magnetic tape, checking 


0, all u's 0 the row sum, and performs the summation indicated 
ml, 


he 








Mode I 


w, rad/sec 41.97 169.81 


ae a 














Mode ~ 5 
w, rod/sec 238.64 35/.90 


























= os 


Figure 3- Noda! lines and frequencies of antisymmetric vibration modes of delta wing 
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TABLE 2 
Antisymmetric Influence Coefficient Matrix (10% in./Ib.) 





5 


te 
lw 
= 


Station 2 


6 7 8 9 10 11 





1 76 1.6164 
2 84, 3.2803 
3 95 4.7735 
4 408 6.280) 
5 2804 11.9258 
6 247, 4588 .6071 #&.7234 1.0510 
7 «4351 -8356 1.1291 1.3552 1.9709 
8 -5975 1.21722 1.6190 1.9718 2.8950 
9 ©7368 1572 2.0448 2.539) 3.8208 
10 9796 1.9431 2.7567 3.5005 5.7745 
11 -1766 =. «3328S =F€ 73 5392 -7962 
12 63112 »=.5903—_ «7984 ~Ss«i«w 9662 0-1. 293 
13 04331 -8281 1.1285 1.3737 2.0463 
14 25436 1.0463 1.4364 1.7610 2.6618 
15 -6397 1.2358 1.7058 2.1066 3.2597 
16 1368 2585 - 3489 24222 6279 
17 e232 ©6596 )~=— «6620 )~—S iw 7510 =21.1168 
18 -3L09 6449 -8719 1.0570 1.5764 
19 4309 = 8161 «21.1049 1.3420 2.0123 


20 -1158 =©.2188 = .2953 ©3576 ~=— 5325 
21 -2072 3913 5281 - 6393 -9520 
22 22920 3=.5517.—S «we TG) =. 9017) 1.3439 
23 -1039 = «.1962_—S's—(i«w 2B = 0 3206~—— «776 
2h 1865 = .3523.—S(«w WY 755 —Ss«w™dS 758 ~—¥ B78 








274 4351 5975 -7368 .9796 -1766 
4588 =.8366 1.1722 1.4572 1.9431 - 3328 
-6071 1.1291 1.6190 2.0448 2.7567 .4473 
7234 1.3552 1.9718 2.5394 3.5005 ~5392 
1.0510 1.9709 2.8950 3.8208 5.775 -7962 
-1618 -28)3 -3911 -483h 6451 1193 
2643 .5212 .7256 #£.8990 1.2013 -2133 
-3911 -7256 1.0374 1.3010 1.7510 2965 
4834 -8990 1.3010 1.6721 2.2960 - 3687 
-6451 1.2013 1.7510 2.2960 3.4662 4.953 
-1193 -2133 2965 - 3687 4953 0940 
2093 - 3789 5304 .6612 - 8886 ‘ 
.2903 5303 .7500 .9415 1.2716 .2268 
3636 =. 667 -9520 1.2081 1.6525 - 2840 
4277 = «7864 1.1266 1.4415 2.0163 o 3342 
-0929 ~=.1670 .2332 = .2910—Ss(—=# 392 -0743 
1650 -2970 4153 5184 - 6986 1318 
Pec | 4169 5841 - 7306 -9867 1843 
-2921 8.5272 =.7399=Ss «9274 s1.2589 ~~ «23:29 
-0786 =. 114 -1976 =.2468 =. 3330 .0629 
-1407 2530 =. 33535 «4413 5956 ~=.1127 
198k - 3567 4985 - 6226 - 8408 -1588 
-0705 .1268 + #.1773 .2213 -2988 .0565 
-1267 -2278 - 3183 -3975 5366 =.1014 





by the right-hand side of Eq. (1) for antisymmetric 
modes [or Eq. (8) for symmetric modes]. 

Using Eqs. (14), (15), and (17) to (19), it is seen that 
both Eq. (5) for removing rigid body modes from the /th 
antisymmetric mode and Eq. (13) for removing rigid 
body modes from the /th symmetric mode can be 
written for the deflection at station g as 


/ 


” LED», YD By, - FX MY By, + 
Ww” s=] ) ] 1 ? l 


n n n 


Eu, >> M, > Byy, — Gu, d 2, 2; &.. PY. = 
s 1 ’ 1 


: get 
> By, (20) 
? l 


The second code carries out the computation indicated 
by Eq. (20), normalizes the resulting mode shape to 
give unit deflection at a preselected station and com- 
pares the improved mode with the starting mode. It 
repeats this iteration cycle until the initial and final 
modes agree to a preset small value indicating that 
convergence is complete. The frequency w, is then 
computed, and the frequency and mode shape are con- 
verted from binary to decimal notation and printed. 


Modification of Dynamic Matrix 


The third code computes the terms in the new 
dynamic matrix for higher modes as indicated by Eq. 
(7), using the numbers already stored on magnetic 
tape. It stores this new dynamic matrix on magnetic 
tape together with row sums for checking purposes. 
The machine is then ready to take the second code 
and compute the next higher mode. This process is 
continued until all the desired modes have been ob- 


tained. 


EXAMPLE 


An example of the computation carried out by 
SEAC using these codes is given in Table | for the 
antisymmetric modes of a delta wing airplane. In 
column (1) of Table 1 is given the station number 
(only the half-wing stations | to 24 are given) column 
(2) gives the mass at the station, columns (3) and (4 
give the x- and s-coordinates respectively and _ col- 
umns (5)-(10) give the mode shapes and natural fre- 
quencies for the first six modes. Fig. 2 shows the 
location of stations on one wing. The nodal lines for 
antisymmetric wing vibration are shown in Fig. 3. 
Similar results have been obtained for symmetric modes. 
The influence coefficient matrix used in the computa 
tion of the antisymmetric modes is given in Table 2. 
The computing time for six symmetric and six anti- 
symmetric modes for this delta wing airplane having 
seven stations on the fuselage and 24 stations on each 
half-wing was 6 hours. This time is sufficiently short 
to make such computations practicable, where they 


could not be undertaken at all with desk computers. 


CONCLUSIONS 


It is concluded that the computation of vibration 
modes by high-speed electronic digital computers 
such as SEAC can be done for structures with many 
degrees of freedom. The computation time is rela- 
tively short and the data preparation time is moderate 
since it merely involves tabulation of the mass values 
and locations in addition to the influence coefficient 


matrix. 
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The Flat Plate with an Angle of Attack in a 
Choked Wind Tunnel 


K. GOTTFRIED GUDERLEY* 
Wright Air Development Center 


INTRODUCTION 


ECENTLY, B. Marschner' has shown, that in a 
R choked wind tunnel, the wall influence on the 
pressure distribution over a wedge is rather small. 
This is of interest for the question of transonic wind- 
tunnel testing. It appears desirable to extend the in- 
vestigations to a body with an angle of attack. 

This is done in the present paper for an inclined plate. 
By repeated application of an idea used previously* to 
determine the flow over this body at Mach Number one, 
the numerical work can be kept within reasonable 
bounds 

Aside from the numerical result which illustrates the 
wall influence on the pressure distribution over the 
plate, the general approach also deserves some interest 
as an example for the treatment of unsymmetric flow 
patterns in the vicinity of Mach Number one; this 
theory was developed in reference 3. 

Some details of the investigations which might have 
interrupted the general line of thought have been rela- 
The first part contains the 
Although 


only loosely related to the present subject, some re- 


gated to a second part. 
leading ideas and a discussion of the results. 


marks have been added which show in principle how 

the flow in a choked closed tunnel can be interpreted 

as the flow with a slightly changed angle of attack at a 

Mach Number different from one. 

Part I-——THE PRINCIPAL IDEAS OF THE INVESTIGATIONS 
AND THE RESULTS 


1) The Basic Equations of the Hodograph 


Let w and 6 be the absolute valute value of the veloc- 
ity vector and its angle with a fixed axis. They are the 
independent variables of the hodograph. In transonic 
work the variable w is frequently replaced by another 


variable n defined by the equation 
1 36m] (I, 1) 


here w* is the value of w at the sonic velocity and y is 
the ratio of the specific beats. By a limiting process 
in which the deviations from the sonic velocity are as- 
sumed to be small-—the details of which are carried 
out in reference 4 
y in the 7, 6-plane is simplified to 


the equation for the stream function 
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tory 


Yon — n'Voo = O I, 2 
Consistent with this approximation the pressure coei 
ficient for sonic condition is expressed by 


Cp = (p — p*)/(p*tw*?/2) = —2(y + 1)7"" +n (I, 3a 


Here p is the pressure, p* the sonic pressure, and p* the 


density at sonic speed. Similarly, one finds for the 


Mach Number .J/ 
M-1 (1/2)*(y + 197% =n (I, 3b 


After y is found, the coordinates x and y in the physical 
plane are determined from the relations 


\ (y+ 1 n° wWe I, ta 
ve = (y + 1)'7-¥, I, 4b 
y=y I, 4c 


Here the stream density at the sonic speed is put equal 


to one. 


2) The Boundary Value Problem 


The flow over an inclined plate in a choked wind 


tunnel is shown in Fig. 1. Sonic lines start at the trail 
ing edge and at the leading edge and extend to the walls. 
A stagnation point is obtained at the pressure side (the 
lower side in Fig. 1) of the plate; it les close to the 
leading edge. Streamlines that lie slightly above the 
stagnation streamline experience first a reversal of the 
flow direction and then spill over the leading edge 
The Meyer expansion that arises at the leading edge 
propagates into these streamlines. This explains the 
shape of the sonic line, which starts at the leading edge 
on the lower side of the plate in a direction normal to its 
surface and then bends upward to reach the upper wall 
Theoretically, one obtains at the suction side of the 
plate, right behind the leading edge, a complete vac 
uum. Naturally, other physical effects will interfere 
here, but their influence on the subsonic part of the 
field is negligible. 

Since the pressure on the suction side is lower than 
the sonic pressure, a Meyer expansion will also start at 
the trailing edge. The sonic line extends into the 
lower part of the flow field and finally reaches the lower 
wind-tunnel wall. 

The hodograph is shown in Fig. 2. Corresponding 
Point 


D is the leading edge, point C the stagnation point, 


points in Figs. | and 2 have the same notation. 


point B the trailing edge. The free-stream velocity is 


marked by point P; at this point all streamlines origi- 
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Fic. 1 Flow over a flat plate in a closed tunnel at choking 


nate. For thestreamline that runs from infinity (point 


P) to the stagnation point C and then follows the sur 


face of the body, we set y = 0. 


that lie, in the hodograph, as well as in the physical 


Part of the streamlines 


plane, above the stagnation streamline, will spill over 
the leading edge. 

The map of the lower wind-tunnel wall is given by 
PQO’, the map of the upper wall by PO”. It is im 
portant that the points P and Q do not coincide (as 
they do for symmetric bodies for reasons of symmetry 
At the present we do not know whether point P lies 
at the upper or lower side of OQ. Our description of 
the flow field refers to a location of the point P as shown 
The relative location of points P and 
It is de 


in Figs. 2 and 3. 
(J is an additional parameter of the problem. 
termined during the computations by the requirement 
that at zero velocity, there is a branch point of the 
streamlines. 
pears so natural that a branch point of the streamlines 


From the physical point of view it ap- 


must lie at zero velocity, that one might expect this 
condition to be fulfilled automatically. A counter- 
example is given by the incompressible hodograph for 
which, from the mathematical point of view, the stag- 
nation point is equivalent to all other points of the hodo 
flow with a hodograph 


possess a physical significance 


graph. Moreover, patterns 
according to Fig. 
although only little practical importance 


case where the branch point of the streamlines does not 


even in the 


coincide with the stagnation point. 

Fig. 3 shows the boundary value problem in the 7, 
§-plane. The origin of the coordinate system is placed 
into point B. By the simplifications made in the 
derivation of differential Eq. (2), the boundary CDE 
moves to infinity. Along the line ABC we have y = 0, 
and the boundary condition along CDE expresses itself 
now as the requirement that y vanishes at infinity. 


845 


' 


For the Mach wave EO”, one obtains the condition 
that no singularity propagates along it, and that it 
does not map to infinity. The condition that point 
C is a branch point of the streamlines is now expressed 
by the requirement that infinity is a branch point of the 
streamlines. As a consequence of the approximation, 
the stagnation point coincides in the physical plane 
with the leading edge of the airfoil—i.e., the region of 
reversed flow becomes smaller and smaller as the angle 
of attack is reduced. 

The width of the wind tunnel may be unity, and we 
allow for the possible displacement of the model in the 
Then 
we have as boundary conditions along the wall: y 
1/2 — c for O’QP and y¥ 


direction toward the upper wall by an amount « 
const const 
along PO”. 


3) Method of Solution 


In the following computations, the assumption will 
be made that the choking Mach Number is very close 
to one—i.e., that for a given angle of attack, QO’ is 


very short. This requires that the wind tunnel is 
sufficiently wide. 
that the ratio ng (38 
ng 1s the value of 7 at point Q and @ is the angle of 
attack. 

Under this assumption one can expect the pressure 


The future computations will show 


2)" must be small, where 


distribution over the plate to be expressible as a de 
velopment with respect to ng—1.e., the solution will 
appear as the solution at Mach one with certain correc 
tion terms added. 

The calculation proceeds in the following steps. 

Step 1. One determines a family of particular solu 
tions that fulfill the boundary conditions at the map of 
the plate and are free of singularities within the region 
bounded by EOABC. Naturally, they will then dis 
regard the boundary condition along O’'QPO”. Since 


the boundary condition along ABC and at infinity is 
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Hodograph representation of Fig. 
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homogeneous, there must be at least one singularity 
in the flow field. This singularity will be placed into 
the point O = O’ =O”. Actually in each particular 
solution, three singularities will be superimposed on 
each other at point O in order to fulfill the additional 
conditions that infinity is branch point of the stream- 
lines and that the plate has a given length. 

The most general solution that fulfills the boundary 
conditions at the plate can be obtained by a superposi- 
tion of these particular solutions—within a region in 
which this superposition converges. The map of the 
wall will lie outside of the region of convergence, but 
the map of the plate is within this region provided that 
OQ is sufficiently small. The coefficients by which 
these solutions are multiplied are naturally undeter- 
mined at this point of the calculations. The shape of 
the body in consideration enters only this part of the 
computations. 

Step 2. We determine the most general solution that 
fulfills the boundary condition along O’QPO", but is 
otherwise free of singularities within the subsonic 
regime. Here the boundary condition at the plate is 
disregarded. For this purpose, one finds first an ex- 
pression which fulfills the inhomogeneous boundary 
condition prescribed along O’QPO”, and then super- 
imposes expressions which have ¥ = 0 along this line. 
These latter expressions must have at least one singular- 
ity in the flow field which will be placed to infinity. 
The solution that fulfills the inhomogeneous boundary 
condition will have the location of point P and of point 
Q as parameters. The most general expression that 
fulfills the boundary conditions at the wall is obtained 
by superimposing the particular solutions which have 
Y = 0 along OQ, on the solution which fulfills the in- 
homogeneous boundary conditions along OQ. 

Step 3. If OQ is sufficiently small, the solutions 


found in Steps | and 2 have a common region of validity. 
There they must coincide. This yields an infinite 
system of equations for the coefficients that were left 
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undetermined in the previous steps. This system can 
be solved by a development with respect to the length 
OQ—1.e., with respect to the deviation of the choking 
Mach Number from 1. 

Details of Steps 1 and 2 are found in the second part 
In the following paragraphs we derive the form that 
the results of these investigations will assume, from a 
reasoning of plausibility, and then proceed to calculate 


the final results. 


4) Certain Particular Solutions of Eq. (I, 2 


The particular solutions which are obtained by the 
steps described previously, appear as the superposition 
of other particular solutions which have been investi 
gated in references 3 and 4. Some of their properties 
will be described subsequently. We introduce new 


independent variables 


E= 7/([(3 2)6|" : (I, 5b 
If one makes the product hypothesis 
y = p 1/12) 2 V/3 G(§, dN) (I, 6 


where X is an arbitrary parameter, one is led to the fol 


lowing differential equation for G. 


re) | ola = | ( ° | ) 
a io ck 
de L -! 2 , 


For a given value of \, the differential equation has 
two linearly independent solutions. They can be 
chosen such that the expressions y which arise from 
Eq. (I, 6) are either symmetric or antisymmetric with 


respect to the line & = —o. This will be indicated 
by the notation G’’ or G respectively. We impose 
the condition that, for & = —o#, G”’ = 1 and G"“ 


behaves as (—£)~**. Obviously, if an expression 
¥(n, 9) fulfills Eq. (1, 2), then also ¥(n, 6), where 0 = 6 
const., is a solution. We shall use 6 = 0 — @—1i. e., 
the origin of the 7, #-system will be placed at point O. 
The variables p and & [Eq. (I, 5)] will then be denoted 
by p and £. 

As can be found from the character of the functions 
G for large values of \ (reference 4, Part 3j), the bound- 
ary for the region in which a superposition of the par- 
ticular solutions found in the first step will converge, is 
given in the subsonic region by the line j = const. 
through point Q (line SQT) and in the supersonic 
region by the characteristics TR and MS. The region 
of convergence for the superposition of the particular 
solutions found in the second step is bounded in the 
subsonic region by the line p = const. through point B 
(line UVB), and in the supersonic region by the charac- 
teristics BA and UN. The shaded region in Fig. 3 is 
common to both representations. 

Now we ask which solution of the form (I, 6) taken 
for the p, —-system will occur in the representation of 
the solution in the shaded region. The conditions 
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FLAT FLATE Wite 
imposed are that along AO and EO no singularities 
propagate and that these Mach waves do not map to 
infinity. The particular solutions can be symmetric 
as well as antisymmetric with respect to the line = = 
—o, This discussion is carried out in the second part. 


One finds the following values of V X: 


For the antisymmetric solutions: 


VA = (3/4) +h I 
p= G4 
Viv = -—(9/4) — 3h\ 
For the symmetric solutions: 
ViA= (1/4) +h 
| 
Vd = —(3/4) — 3h) 


The value of p for point O may be denoted by py. 
Expressed by the angle of attack @, we have py = 
9 4) -A°. 
by the last discussion will be multiplied with suitable 


The particular solutions (I, 6) obtained 


powers of pp and then denoted in the following manner. 


Antisymmetric solutions: 


Y— (5/6) —» = +(B/ po)“ ~" -G™ 4 8,[(9/4) + 3h}?} 
(I, Sa 
1/6 } 
y ‘ F = +(p pi x 


+ h)?} (1, 8b 


Symmetric solutions: 


¥— 1/3) —h = (B/po) GF &,[(3/4) + 3h]?} 
(I, Se 
Vi = (p Po) . -G &, | l } a h\?; (I, Sd 


The upper and lower signs apply for 6 > 0 and @ < O 
respectively. The subscripts of y give the power to 


which p/ pp is raised. 


5) Solutions Fulfilling the Boundary Conditions at the 

Plate 

Among the particular solutions just listed, those with 
a negative subscript are singular at point O. The 
particular solutions which fulfill the boundary condi- 
tions at the plate will be denoted by capital VY with a 
superscript which corresponds to the strongest singular- 
ity which occurs at point O. 

Let us first consider the expression which describes 
the flow with a free stream Mach Number one. Here 
the singularity which occurs at point O is given by W_ 5. 
To fulfill the additional condition that infinity is a 
branch point of the streamlines, one needs another 
solution which has y = O along the map of the plate. 
The only singularity which can be used in this case is 
¥_1/3, (reference 3, p. 50, 760). 
Furthermore, the length of the plate in the physical 
plane is prescribed; this condition will be fulfilled by 


or reference 5, p. 


multiplying the hodograph expression with a suitable 
factor. Thus, the flow with a free-stream Mach Num- 


ber one has in principle the form 


y = (L/)I)-Vv (I, 9a) 


AN 


ANGLE OF ATTACK S47 
where 
We = ise tat; °° ¥ ; 
2. (1/6) 4 ; sles y 6) + 3 oT 
> ¥ a3 “y I, 9b 


t 


Here L is the length of the plate, and / and the a's are 
constants—/ 1s the length obtained from the solution 
y ””. The a will be determined in 


rheir superscripts are those of the functions y in whose 


the second part. 
representation they occur; their subscripts are those of 
the function ¥Y whose coefficients they are. The terms 
in the sum give the regular terms which will naturally 
occur in the representation of the solution in the vicin 
ity of point O. 

The other particular solutions which fulfill the bound 
ary conditions at the plate are constructed in a similar 
manner. One first finds a solution which has y = 0 
along the map of the plate surface and is singular at 
The singularity may be given by ~_:5/6 
or by y-1/3)-a (A = 1, 2, 
to this particular solution of the same type with w_ , 
as singularity at point O, in order to make infinity a 


point O. 
One superimposes 


branch point of the streamlines. The expression thus 
obtained will in general change the length of the plate 
Mach 


multi 


when it is superimposed to the basic flow at 
Number one. Therefore, the expression (I, 9b 
plied with a suitable constant is superimposed on it, 
such that the resulting expression does not give a con 
tribution to the length of the plate. Thus these par 


ticular solutions have the form 


5/6 } 
— y = /e . & 5/6 y T 
t — t 
a 3 y 1 T >» aii % 4 
ft) 
Wiis ‘ie B > Ap /3 y I, % 
0 
i " ‘ 
V = y { T Q@ ( . 
13 / = j 
a 1/3 -W ee >: A116 " , 4 
o 
“= 13 
Vays +3) + DL anys ¥i3 (1, 9d 


6) Remarks Concerning a Body of a More General Shape 


For a more general body, the solutions which fulfill 
the boundary conditions at its surface have, in principle, 
the same form. The quantity p) which occurs in the 


definition Lis a charac- 


(I, 8S) can be fixed arbitrarily. 
/ 1s the value of this length 


starts the 


teristic length of the body; 
as obtained from the solution W One 
calculation of such flow patterns with the flow over the 
body at Mach Number one. At point O the same 
singularities will occur as for the flat plate.’ Thus 
the solution will have the form (I, 9a) and (I, 9b 

naturally the constants in Eq. (I, 9) will be different. 
Next those particular solutions are determined, which, 
if superimposed to the basic flow, leave the surface of 
linearization of the 


the body unchanged. Here a 


boundary condition in the hodograph is required.* 
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The solutions thus obtained have the form of the ex- 
pressions (I, 9c) and (I, 9d). 

These considerations yield a simple interpretation 
for the solution ¥~*” in the case of a body of an arbi 
trary shape. Let us assume that the flow over an 
arbitrary body at Mach Number one is known, and 
also that the angle of attack is changed by a small 
amount. The change of the angle of attack 1s convent 
ently expressed in the hodograph by the shift of point O 
along the sonic line. Denoting the change of the angle 
of attack by Aa, and considering for a moment the 
function W~”° 
the strongest singularity of Eq. (1, 9b) in the form 


as function of 7 and 6, we can write 


[7. l(Aa)]-y 5/6(n, 0 — Aa) 


Consider now the first order terms in a development 
with respect to Aa. The predominating term is found 
from the last equation as 


—[L/l(Aa)]-(OW_ 5/6/08) 


Thus the change of the angle of attack is equivalent 
to the superposition of a solution which leaves the sur- 
face of the body unchanged and has at point O a singu- 
larity given by (O~_; 4/06) together with other singu- 
larities Y_5/6 and Y_,,/3 which appear automatically 
in the solution of the boundary value problem. 

Now (OW_;6/ 08) is proportional to Y_4/3 as can be 
seen from reference + where the same particular solu- 
tions are discussed in a slightly different form. The 
factor of proportionality is best obtained from the be- 


havior of the functions along the line — = — ©. There, 
according to Eq. (I, Sa) and (I, Sb) and the definitions 
for G ee 

= 5/6 = 3/2 = 5/6 
y 5/6 = \P/ Po) a eal = (p/ Po) x 


[(3/2)6 (—n)*’7] 
Because of the definition of p, [Eq. (1, 5a) | 


(0p/06) = O for@=0 
1/2 


and thus, (OW_5/6/00) = (3/2)-(p/ po) *+ pp 


According to Eq. (I, Sc), this is the behavior of 
3/2) > po sy 13 at — = —o 


Thus it follows that the effect of a small change of the 
angle of attack Aa is expressed in the vicinity of the 


body by a superposition of 


— (3/2)-(L/I)-e~**- a */*-Aa (1, 9e 


One is tempted to interpret also the other particular 
solutions W~''’°, w~7*, Eqs. (I, 9c), and (I, 9d) in a 
similar manner—1.e., as higher order terms in a de- 
velopment with respect to the angle of attack. For a 
general body this leads to difficulties because of the 
linearization required in fulfilling the boundary condi- 
tions at the surface of the body. In the case of the 


plate where the boundary condition is /inear, and where 
the solution is known for any angle of attack, the par- 
ticular solutions ¥~ °°’ ~" and ¥~"’® ~* are in princi- 
ple contained in the expression for the flow at Mach 
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Number one. Obviously, expressions of the form 
(II, 16) which are used for the determination of the 


. . 56 . 
particular solutions V etc., arise by repeated 
differentiation of the expression (II, 16a), which gives 


the flow over a flat plate at Mach Number one. 


7) Solutions Fulfilling the Boundary Conditions at the 

Wind Tunnel Walls 

The solutions which fulfill the homogeneous boundary 
condition Y = 0 along O’QO", have one singularity at 
infinity, given by a particular solution—Eq. (I, 8 
with a positive subscript. Besides, they contain a 
superposition of functions with a negative subscript 
which are regular at infinity. Particular solutions 
which are antisymmetric with respect to — = 
fulfill the condition Y = 0 along O’QO” immediately 
Denoting particular solutions which have y = O along 
O’QO" by capital Y with a superscript corresponding to 
the subscript of the particular solution which is singular 


at infinity, we thus have 
gt" ¥i1/G) 4 ; I, 10a 


For a singularity which is symmetric with respect to 


Seri 


= —o, the entire expression yY will be symmetric 
with respect to this line. In such particular solutions 
no antisymmetric terms (I, Sa) will appear. 


3 in (pa po)” = si Xx 
Wi-1/3 (I, 10b 


Here fg is the value of jf at point Q. The constants 3 
play a role similar to that of the a. It is not surprising 
that in these particular solutions the ratio (je py) ap 
pears as a parameter. 

The determination of the inhomogeneous particular 
solutions is simplified by the fact that for a sufficiently 
small model, the ratio of the length PQ to the length 
OQ is small. In the vicinity of point Q, the differential 
Eq. (2) for Y can be simplified to 


Van v NQ Woe = U 


nq is the value of y at point Q. This is Laplace's equa 
tion. Ina region where this equation is applicable, an 
expression which fulfills the inhomogeneous boundary 
conditions along OQ can be represented by 


y = (1/2) — c — (1/n) X 


Im log (*\ n + 16 \no ae na| — V |no| — InP 


Here np is the value of 7 at point P. The second root 
in the argument of the logarithm is positive; for the 
path O”’PQ, the first root is positive, too. The sign 
of the difference between the two roots is positive be- 
tween P and O”, negative between Q and P. As one 
travels around point P on the path indicated in Fig. 
3 by the dashed line, the argument of 


| 2 1/2 | 
V 7T 16 Ne TT ie! = \ Soi —~ IGP 


increases from 0 to z. Thus we have y = (1/2) — ¢ 
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rm along PO”, and y -(1/2) — c along PQ. As one derived in the second part. It can be characterized 
the travels around point Q along the path indicated, the by its behavior for 6 0, 7 < ne 
ted first root changes its sign and the expression 2 

. l/r)-( ine np ‘(in No 
ves f 
*% sill 10 V |" at Phis particular solution is symmetric with respect to 
remains negative, thus y 1/2 c along QO the line & = It will not contain terms which are 
the singular at infinity. Thus one will expect it to have 
If V ne} — |ne|| < |V at 28\ng| " + |ne the form 
iry 
: then the expression I, 11) can be developed with re — ,; 
u Vi 2, 0 Pa Pp ¥ 
spect to N Ne ae np 
2 I, 14 
a y 12)-— <¢ 1/27) Im log (n + 26 \n@ t 
pt ne - (1/mr)-V ine np| X rhe expression yy derived in the second part behaves 
ns 1/2 1/2 - = 
” Im (n + |n@| + 16- \ne - t, 42 lor 0 0 and p > pg as 
ly Thus the desired inhomogeneous expression is repre 1/V 3)- Inei . n Ne for 
ng sented for small values of the distance PQ by the super 
to position of two singular expression whose local behavior 8) Solution of the Boundary Value Problem 
ar is described by the analytic functions occurring in The most general expression which fulfills the condi 
Eq. (I, 12). tions at the plate surface is given by a superposition of 
The logarithmic singularity in Eq. (1, 12) combined the solution (1, 9), which may be written in the form 
- with the constant 1 2 fulfills the boundary conditions 
am y = 1 2 along O"Q and y = —1,2 along O’Q. This y I E ws FT A 6) — kh yy — 6/6 
; solution (except for a factor 1/2) has been considered ; 
‘IC : : : 2 ; — 
. in reference | in conjunction with the flow over a sym : ar ed _ V 4 ‘iy 
metrical body in a choked tunnel. For jp > fg, this 0 
expression is represented by Here the A’s are constants which will be determined 
subsequently. The most general expression which 
yy = we Y -(5/6 hn’ (Pe/ p : V_656 j fulfills the conditions at the wind-tunnel wall is given 
D — 13 by a superposition of ¥, (Eq. I, 13), of wy (Eq. I, 14 
’ ” ~ . - ° 
3; of the expressions (I, 10), and of a constant —c which 
1. The upper sign is valid for 6 > 0, the lower one for takes into account a possible displacement of the model 
D 4< 0. The constants y can be found from reference | from the centerline of the wind tunnel. 
iS 
! ae © a DS ae wilt 
Baas, at ; ' , : ¥ YI C Wir Lo : 
al | I (0 3 i h| | | 1/3) + h| 0 
\ PE ' at agg = x , 
h w-V 5 r(2/3 r(1/3 > 2 WV , I. 16 
, r(5/6 ri /6 0 
ul ° 
r{(11/6) +h] P[(7/6) + kh] where C’ V3 4) In@ “(ine NI : 
| 2 / aa) eo 23) th I, 16b 
y 6-11 6 D6 ™ 
xv 3 5/6- ( 5/6) +h = Ss is ol : 
eS 2s a ' The coefficient C’ of the second term and the B’s are 
; 1/3-4/3...C1/3) +a unknown. In the common region of validity the ex 
. 1 6-76 1/6) +A pressions (I, 15) and (I, 16) must coincide. A sufficient 
V » -° ° - ° . e 
I, la condition is that the coefficients of the individual par 
: ‘ sag E , ticular solutions (I, S) are the same his vields the 
rhe expression (I, 13) does not contain terms which 2 ; 
; ee following equations 
are singular at infinity. ; a 
: : , terms with W_ 5, hy (A iz 
The other singular term in (I, 12) yields y = O 
along O’'O and O"O. A representation for it will be é -A ‘he ? 7 ( ‘(pa Pp I, l7a 
t 
Cc ° 
terms with y_,; 
n ; 
Lplit+ = AP? -* ean ™ pF BO tg gO ne y ltalod”” Gh 
> 0 | 
terms with W_ (4/3 Lh 0, 1 
, 3 j ye yy +i = ; Pi ae 3) + (4/3) +3 » 
(L/l)-A7~“° = C'-§_ (4/3 n° (Bo po)” + > B’’”’-B_«43 n° (po/p [, 1%7e 
0 
terms with Y_1/3 
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5/6 1 5/6 5/6 ! 1/3 k 1/3 } 
(L/l) E 1/3 as z, A *Qu31/78 a. } A *@..4/3 | = 
} A 0 


i « = 1/3 k/3 k/3 (a k/3) 4 3 — 
C’-6_1/3 (fe/ po) ss » Br" - B13 {PQ po ; I, 17d 
k=0 


terms with y, 3, h = 1,2 
(Li) E = + Zane ae ge | B*/* I, 17e 
a h ] h 0 
Since yy = | [see the definition (I, Sd) and the first paragraph of Part II], the coefficients of Y give 
(L/l) E alll ay i +s. " -F 2 ates * sa se | =-c+B I, 17f 
k= 0 
Finally the terms with y'° * “* yield 


A l 


Ch |< 1/6) + (h/3 sili + > A 5/6 ‘iis tie dai 5/6 4 


a ”C”C—ClC ha |-30"" , Cie 
0 
The first of these equations gives directly the A~'” ~" (h = 1,2...). Furthermore, the quantities B'' ° ap- 


pear only in Eq. (1, 17g), thus they can disregarded until the other unknowns have been determined. 
Important for the pressure distribution over the plate are the quantities A—they yield the correction terms to the 
. . 5/6 - . - . 7 . ? 
basic flow given by WV . Among the value of A obtained from Eq. (I, 17a), the one with the lowest power of 
] 


po/pvois A~''’®; it is proportional to | jy po|''°, the other terms are presently unimportant. 


For further computations we make the assumption, justified by later results, that the coefficients A » ~" are 
multiplied with powers of jg/ po larger than 5/6, Then Eq. (I, 17b) yields 


L/l = y~5/6- (Be ay” 4- higher order terms I, 18a 
or Bo/po = (y-36) -L”?-1-©” + higher order terms I, ISb 


Accordingly L// can be considered in the following computations as O(jg/ po) One obtains from Eq. (I, 17e 


B® = Ol@a/p)*) 
and from Eq. (I, 17f) 
B® = c + O[(fe/po)”’*] 


With these results and Eq. (I, 18a) one finds as the lowest order terms in Eq. (I, 17d 


°o l o 


we 5/6 5/6 Yo = l ’ 0/— l 
Y —5/6(PQ’ Po °a-1/3 = C6_, 3(Pa’ po) + CB_1/3 (pa po) 
- 1/2 5/6 : ( 
Hen Cc’ [(Be/ po) “~* ¥-5/6°@-1/3 C°B-1/3 
ence = _ 
0-1/3 6 | 


Inserting this into Eq. (I, 17c) one obtains 


= (11/6 + h 5/6 « = 1/3 + h ° 
L (4/3) —h PQ Y -5/6° @-1/3 "66 1/3 h Pe ( 9 9-(4/3 j 
A = = = °C B_-1/3 2 = B~(4/3 
l po 0-1/3 Po \ 0-1/3 


I, 18 


Presently this expression is of importance for h = 0. The solution for the flow over the plate with the lowest order 
terms for the wall influence included, is obtained by inserting these results into Eq. (1, 15). The quantity pg pp» is 


given in formula (I, 1Sb) 


L 5/6 Di. j 6 6 5 /¢ 6 3 3 
y = E aad that + Pa ( ul -W 11 = a_1/3 * ioe ; -W “ x 


l Po \V-5/6 0-1/3 
~ 1/2 © 
p € 0 3 : 
( ‘) = (s 1 3’ a — G4 :) -W | I, 19 
Po 76/6 0-1/3 


If higher order terms for the solution were determined, then one would expect fractional powers of jg py» to occur, 
even if the displacement of the model from the wind-tunnel axis is disregarded. The values of the y are given in 
Eq. (I, 13a), and the other constants are derived in the second part—one has 


S -11/6 80 a 19 
7 J —o/6 . 
Y-5/6 = ; aoe =—; a1, =- (II, 48) 


5e-V 3’ Y-5/6 ‘4 16 
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6 13 ] I] 
ian. 2 — 

3 3°r(5/6 
f 34> - I] Ola 
r(1/3)rd 2 
3 3 (5/6 
i — (II, 6lb 
8 r(1/3)rd 2 
{9 13 ae r(5,6)-T(5/3)) 
= 34(7+ “Bo ’ 
Ga 7 ; 2 r(5/2 y sae 


and thus finds 


and 


For simplicity numerical values may be inserted as far as practical. 


So/p = 0.988-6,° °° L*" (I, 20b 
Since Po = (9/4)87 (I, 20c 
we have Ze = 2.2234," .L°"* 
and ne| = Vibe 1.305 @°/* 7°’ (I, 20d 


Furthermore, because of Eq. (1, 3b), we obtain for the 
choking Mach Number J/¢h01 


hok 


1 — Ment (sae. 7 (I, 


20e 


9, The Pressure Distribution Over the Plate 


According to Eq. (I, 4), the coordinates in the physi 
cal plane arise from a given y by a linear operation. 
Therefore, if one superimposes different expressions y 
on each other, then the corresponding values of the co 
ordinates which belong to a given value of n and @ are 
also added. We apply this to the coordinates x which 
are found along the side of the plate BC for a given value 
of7. The contribution to x belonging to a particular 
’ The 
0 is assumed to be the leading edge of the 


solution V " may be denoted by x‘~” ° 
point v 
plate which is, in the present approximation, identical 
with the stagnation point (ny = — ¢ With these nota 


tions one obtains 


(1,23 


The solutions x * have been constructed in such a 


manner that only x ”’” gives a contribution to the total 


length of the plate. The total length of the plate is ob- 
i.e., determines the value of 


The 


tained if one sets 7 = 0 
x for the rearward starting point of the sonic line. 


Then 


us| (1. 
J 


We choose y = 1.4. 


al Bo 1/2 
Vv a -¢- 0.9094 W 
p 


length /, 


20a) 


introduced previously, is given by the value 
of x " for n = 0. 


The last equation can be written as 
. sg. (2 zs" 16 
8 l p vi 32 


= 1/2 i 
PQ x 

r ( ) -O.9094 
Pp. 


The quantities , " depend only upon 9, (36), 2)” 


as can be seen from the second part. Introducing 
for brevity 
9 = 9/(30/2)°" (I, 22) 
i 
and (9 \ 


1/2 13 
Q \ 
0.9094 [, 23 


where € is a small quantity, introduced to show that the 


expression €-/s is a correction to fj), Eq. (I, 21) is writ- 
ten as 
x/L = filn) + e-foln (I, 24 
Since 7 is related to the pressure coefficient by Eq. 
I, 3a), we are interested in a result of the form 
n = o(x/L) + €-go(x/L) (I, 25 


For this purpose the last equation is inserted into Eq. 
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This equation is developed with respect to e. The 
zero order terms yield x/L = fi[gi(x/L)]|—i.e., g; is the 
inverse of fy. 

The first order terms yield 


[ol Jew) +«[a)]=o 
| 


x fo(gi[x/L]) 
s2 7 € , 
a & Ih (gi [x L]) 


According to the results of the second part we have 


x oF t5/2) 


- a re | uw’ ~(1 + u3) "du 
'y r(5/6)-F(5/3) « 


The solid line in the upper part of Fig. 4 shows the in 


verse function g thus obtained. Furthermore, we find 


€-fo(#) Se 30 (—7)3 35 
a = jie ~ = —+ i+ 
fi’(n) po (71+ (—7) LS 
41/2 _9 
( (**) 0.9094 ( *) (—¥ 
pv ss 


In this expression one must, according to Eq. (I, 26), 
replace » by g; and then insert the result into Eq. 
(2. Bo. 


sure distribution 


n (*) [ Bo (2 g,3 3) 
y= 7; ioe i a 
(30) /2)° ~ Vee po \il 1 — gy 1S 


= 40/8 - 
p 

c( ‘) 0.006 
po 


One thus obtains an expression for the pres 
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Fic. 5. Functions needed for the evaluation of the pressure 
change caused by the wind tunnel walls 
The value of pg py» is obtained from Eq. (I, 20b). The 
expression g; } (30/77) [gi°, (1 — g:*)] + (35 48)}, which 


occurs here, is shown in Fig. 5. 

These formulas give the pressure distribution on the 
subsonic side of the plate. It has been found previ 
ously that the free stream pressure distribution on the 
suction side is obtained with sufficient accuracy by 
analytic continuation of the solution obtained to super 
sonic speeds; we assume that the same is true in the 
present case. 


Accordingly, we have for the suction side: 


F(z x 31'(5/2 i aie 

(n) = = _— = = “ ee 

; L r(5/6)-T(5/3 

= du n> 1 


The function g(x, L) is then the inverse of this expres- 
sion (solid line in the lower part of Fig. 4). The pres- 
sure distribution on the suction side is expressed 


by 
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tat : Frese 
[he expression 2, (30. 77) [g:° 
is included in Fig. 5. 


Discussion 


The results obtained here can be discussed from two 
different points of view. First, the computations give 
some data which illustrate the effect of the wall on a 
plate in a choked wind tunnel. Secondly, the present 
investigations illustrate for an unsymmetric body the 
treatment of a flow field which deviates slightly from 
a flow field with a free stream Mach Number one. It 
thus can serve to illustrate the essential features of 
such flows. 

Fig. + gives a comparison of the pressure distribution 
at Mach one with the pressure distribution in a choked 
tunnel for a plate whose length is equal to 0.1 of the 
wind-tunnel width and with an angle of attack of 0.1 
radian (5.7°). Actually the quantity /(3)/2)°” is 
plotted; according to Eq. (I, 3a) it is proportional to 
the pressure coefficient. The solid lines give the pres- 
sures in a free stream at Mach one, the dashed lines 
give the results from the choked tunnel. The differ- 
ence in the pressure distributions is rather small, which 
shows that even unsymmetric bodies can be measured 
in a choked tunnel without too strong an influence of 
the wall. 

According to Eq. (I, 19), the angle of attack and the 
plate length enters the correction by means of the value 
of fg po, Which is determined by Eq. (I, 20b). It is 
natural that the increase of the plate length at a fixed 
wind-tunnel size will increase the wall influence. That 
a reduction of the angle of attack will increase the rela- 
tive error may be somewhat surprising, but it is ac- 
tually clear from the transonic similarity law. 

The expression for the choking Mach Number has 
the form found in reference | for a symmetric wedge 
in a choked tunnel. If one denotes by LZ the length of 
the wedge from the nose to the shoulder, then the ex- 


pression found in reference | assumes the form 


1.488-@"""- LE 


l — M enor 


Even the coefficients at the right side of this equation 
are quite close to the value given in Eq. (I, 20e).. That 
the symmetrical wedge chokes at a slightly lower Mach 
Number than the flat plate can be ascribed to the in- 
fluence of the spill-over which occurs at the leading 
edge of the plate. 

Within the approximation, a displacement of the 
body from the wind-tunnel centerline is exactly equiva- 
lent to a change of the angle of attack. A general 
formula will be given in connection with the flow over 
a body of a more general form. 

The flat plate set at a fixed angle of attack may be 
considered as an example for an unsymmetric body. 


Eq. (I, 19) can immediately be considered as the solu- 


AN 


ANGLE OF ATTACK O35 
tion for a general body if the proper expressions V 

v~**, and W~''’® are inserted. The properties of the 
body enter this solution in the form of these particular 
solutions, by the quantity / which is a characteristic 
length of the body and in the form of the coefficient 
a_,3 °° whichisrelated to the unsymmetry of the body. 
The latter quantities are immediately connected with 
»® ete. (For the plate, the unsym 
Thus the 


the solutions V 
metry is caused by the angle of attack.) 
extension of the present approach to other bodies would 
require only the determination of the function ¥~°”° 
vw? and wo'®, 

One general conclusion can be drawn. The expres 
sion (I, 9e) together with Eq. (I, 20a) show that the 
displacement from the wind-tunnel axis is exactly equiv 
alent to a change of the angle attack, Aa, which 1s 
given by 

Aa = —(2/3)-0.9094-f9'*-« 


Here jg can be expressed by the choking Mach Number 
and one obtains for a body of an arbitrary shape 


> 


Aa = —(2/3)-0.9094[277 /(¥ + 1)](1 — Menor)?’ * + = 
—0.716-c(1 — Meno.) 
This effect is quite small—e.g., for a model with a 


choking Mach Number 0.9 and a displacement from 
the centerline of 0.1 times the tunnel width, one finds 
— 0.00227 

-0.15 


Aa — 
or Aa = 


In reference 3, the flow in a choked tunnel has been 
used as an illustration for flow patterns which are close 
to a flow with a free stream Mach Number one. It 
appeared plausible that no matter which boundary 
conditions are imposed at a distance from the body, the 
solution would always appear in the form (I, 19 
naturally with different values of the constants which 
depend upon the boundary conditions imposed at a 
distance from the model. The flow pattern which is 
most interesting in this regard is the flow at a high sub 
sonic Mach Number. Naturally, a displacement from 
a fictitious center line will be without influence in a flow 
field which extends to infinity, thus the solution for a 
flow with a high subsonic Mach Number can be ex 


pected to have the form 


y = 


an : T Co*Q-1/3 om 


WV MS)] ‘t. 2a 


(L/1) [Wo + (Bar po 


The constants c; and c. are independent of the form of 
the body considered. jy is the value of p which be 
longs to the free stream Mach Number. py can be 


chosen arbitrarily——e.g., it can be defined by a charac 
teristic angle of the body, then it would be possible to 
apply afterwards the transonic similarity law. It can 
also be set to equal unity. A direct evaluation of this 
formula is impossible in most cases, even if the con 
stants c; and c. were known, since for each body con 
figuration one needs the functions V 5/6 ew "and 
wv *°*. But it is useful if we compare the influence of 


different boundary conditions imposed at a distance 
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from the given body, for then these functions WV are 
the same. In particular the question can be answered 
as to which free flow corresponds to the flow in a choked 
wind tunnel. 

For a symmetrical contour, the second terms in 
Eqs. (I, 19) and (I, 28) are zero, and the expressions 
will agree if 


A 


C:°pu = Pa’ Y—-11/6/ ¥—5/6 
Naturally, this formula is applicable only if the con- 
stants c; and y_11/4, y—-5/5 have the same sign; then we 


have the relation 
| ee (8 Big yi a/erescidt” 


To obtain agreement for an unsymmetrical contour, 
a change of the angle of attack must be admitted. If 
both the deviation of the free stream Mach Number 
from one and the angle of attack change are small, 
we expect a simple superposition of the effect, and one 
obtains them by a combination of Eqs. (I, 19) and 


(I, 9e) 


L ee DB / : ie 
y= E ie ‘ge € er ey Oe sia” & 


Y-5/6 


6 1/3 ‘ 3 9 2 
pos) — po 2. daw? (I, 29) 
6_1/ Z 
1/3 - 


The expressions (I, 28) and (I, 29) will agree if 


~ 
l 


In these relations the constants c; and c» are unknown 
but they can be determined once and for all by experi- 
ments. But furthermore, the quantity a_, 3 °/® enters 
the solution, and this quantity is determined by the 
body contour. On the other hand, this quantity is the 
lowest order term which causes the pressure distribu- 
tion along the lower and upper wall of the wind tunnel 
to be different. Thus if the pressure distributions 
along the walls agree for two different bodies, then these 
. .. and, therefore, 


bodies have the same value of a_/3 
the same Mach Number and angle of attack correction. 
Thus the corrections to be applied can be found from 
the wall pressure distributions. The Mach Number 
correction is found from the average of the pressure dis- 
tribution along the lower and upper wall, and the angle 
of attack correction from the difference of the pressure 
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wall pressures. Furthermore, one determines for the 
same body the pressure distributions in a free flight 
(or in a much larger wind tunnel) at different Mach 
Numbers and for slightly changed angles of attack, 
The Mach Number and that angle of attack for which 
the pressure distribution over the body in free flight 
agrees best with the pressure distribution in the tunnel 
gives first, the coordination between the wall pressures 
and second, the effective Mach Number and the angle 
of attack correction. 

If one accepts the linear dependence of Avon a_,,3 ‘ 
one does not need further experiments. An important 


5/€ 


extension is obtained by the application of the transonic 
similarity law. 

Naturally, the procedure of determining Mach Num- 
ber and angle of attack corrections described here is not 
restricted to a choked flow but can be applied also for 
other Mach Numbers.® 

The corrections considered here are the lowest order 
terms in a development with respect to the model size. 
Even if the first order effects are equal to zero, the higher 
order effects will differ. In the present case one must 
keep in mind that the lowest order terms are propor- 
tional to pg, while higher order terms would —— according 
to the form of the system (I, 17) contain powers of 
Bo, ete. Since the powers of pg are thus rather 
closely spaced, it nay happen that higher order terms 
have an influence comparable with the lowest order ex 
pressions. Here it is rather reassuring that even the 
lowest order influence is within reasonable bounds for a 
body of technically acceptable dimensions; otherwise 
there might be doubt if the results of a wind-tunnel 
test, even when corrected in the manner described 
above, would have a practical significance. 


Part II—-DETAILS OF THE INVESTIGATIONS 


1) Formulas for the Particular Solutions of Eq. (5 


(a) In Eq. (I, 6) of the main text, particular solutions 
for y were given by the hypothesis 


y= p-/l2 V/3) G(E, d) 

Depending upon whether the function yY which arises 
from this hypothesis is symmetric or antisymmetric 
with respect to the line £ = — o, the functions G will 
be denoted by G’”’ or G’“’. As mentioned in the first 
part, we impose at & = — © the condition that G 


1 and that G'” behaves as (—¢)~”~. The values of 


VX which are used in the present analysis are deter- 


distribution. 

Practically, one will measure the pressure distribu- mined by requirements imposed to the function wy at 
tion over a sample body in a choked tunnel and also the € = — o and along the Mach wave é = 1. 

According to Eq. (19d) of reference 3, the function G™ is given along — = 1 by 


G!) = r( ) ( r'[(2/3)VYA]-sin} r[(4/d/3) — (1/4) ]f 
a 2) \P[(W/d/3) + (5/12) ]-PIA/X/3) + (1/12)] 
r-(2 3)VYA]-sin} r[—(/A 3) — (1 Hy} GQ — vim tare 1400 —¢ ) 
P[—(+/d/3) + (5/12) ]-P[—(V/0/3) + (1/12) ] 


(1 — €3)-(V/3) + 1/12) [1 + O(1 — &)] + 
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Using the definitions for p and &, Eqs. (I, 5) one finds 


1 r((2/3)-VaA]-sin} r[(Wd/3) — (1/4) ] 5 
eae y+ var ) (7 -W/SI- ia(VWA/3 I; +00 —e)]4 
“ 2) \T[(VA/3) + (5/12)]-T[(VA/3) + (1/12)] 
r{—(2.3)-VA]-sin} r[— (V/A,/3) — (1/4) ] 
= ee me (1 — 8)2/)- va [1 + O11 — &)]) (II, 1 
C[— (VA/3) + (5/12)]-Tl- (VA/3) + (1/12)] 
Similarly one finds for the functions G’” the following representation in the vicinity of & 
l r[(2/3)VA]-sin} r[(4/d/3) + (1/4) 5 
y= Ot VAL “ ~ l+ Oo —z 
2 r[(V/A/3) + (11/12)]-Tl(VA/3) + (7/12 
r(— (2/3)-VA]-sin { x[- \/3) + (1/4) ]5 
V \I t a. Di 1 — £3)(2/3)-VA 1] + O(1 £3) ] II, 2 
r[— (WA/3) + (11/12)]-Tl— (VWA/3) + (7/12) | 
b) For the particular solutions which represent the solution in the vicinity of point O, we have the requirement 
that along the characteristic £ = 1, neither a singularity propagates nor this line is mapped to infinity. Obviously, 


the second term in the parentheses of Eqs. (II, 1) and (II, 2) describes such a behavior because of the presence of the 
factor (1 £\2 3), so we obtain the condition that the coefficient of this term must vanish. This happens 1 
either the sine function in the numerator is zero or the T’-functions in the denominator are infinite. One thus ob 
tains the following values of 7/A 

3/4) +h and -(9/4) — 3h h oO 1.3 for the antisymmetric solutions li, 3 


1/4) +h and —(3/4 3h h 9, 1,2 for the symmetric solutions II, 4 


[These are the same values as given in the main text. 
For these values of 1/A we need the behavior of the functions G along the sonic line. Formulas (1c) and (19d 


of reference 3 give 


G | = - ¥* A 
2TL(V/A/3) + (11/12) ]-Tl[—(VWA/3) + (IL 12 
( V/Xr l V/r cee) } (2/3)-T(1/2 
} + - : ge) 4 x 
12 ; 12° 3 IT (VrA3 7 12))-Tl—(VA/3) + (7/12)] 
\ 5 N 11 4 
1 — g*)c/t2 v4 er(-% + —> . :) (II, 5 
3 12 12 3 
and 
3 1(4 3)-1 1/2 
G - - 1 — ¢3)a/1 v»/3) XK 
PT l(V/dA/3) + (5/12)]-T I V/d/3) + (5/12)] 
N h 72 3 F(2/3)-1 1/2 
3 12 3 12’ 3 IP ((V/A/3) + (1 /12)]-1 V/r/3) + (1/12)] 
() — gare vwa.g-r(— V4 >, WO uve) — 
3 12 12 3 


Inserting the above values of ¥/ A one obtains, with some simplifications in the ['-functions, as representations in the 


vicinity of the sonic line 


(a) J | (") > v3 P[(5/6) + hy /T(5/6 a 
ee 42,1 = - dh = (—)"-(3/8)-2 - (—)*-(3/8)-2°°"-& X 
JI r((5/3) + Al/T(5/3 
r((7/6) + h]/T (7/6), 
. ) II, 7a 
> ((4/3) + h]/T(43), 
j } Sy v3 F[3) + Al/Td/3 ; 
G” <¢ + 3h = — (—)"-(3/2)-2 - 3 /4)-2°"-E &X 
() lt ri(7/6) + h]/1(7/6 
r{(2/3) + A] 1T(2/3)} 
— " ——. (II, 7b 
r{[(5/6) + A] (5/6); 


r[(43) + A) T4383 
[b= - mee a II, 7d 
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j | (*) 27 Fr rl(5 6 h\/T(5/6) 1 on r((7/6) + h]/T(7/6 . 
oe Seb = — 3h 2 =i (2) ne§ S 
\ 1.6 r[(2’3) +h] T(2/3) 4 rf(/3) + kA) rda/3 
ry (2 } h| r(2/3 


ri(i/2) + ayPrdy2 
9 . , {(4/3) + / } re r{(a/3 / ae 
G 181 (3) on | | » oa TI eee ¢ b.92/8 pg I . II, Sd 
| Be | IS r((7 0) — h] (7 6 T} a/o — h| Pr(5 ) 


Eqs. (II, 7a), (II, 7b), (II, Sa), and (II, Sd) originate from each other by a suitable choice of hk. The values of V/A 


given in the first parts of Eqs. (II, 3) and (II, 4) lead to the expressions (II, 7b), (II, 7e), (II, 7d), and to (IT, Sb), 


(I1, Sc), and (II, Sd). 


In (II, 7c), (11, 7d), (II, Sb), and (II, Sc), either the value of G or its derivative with respect to & vanishes at (), 


Actually the corresponding solutions y are polynomial in é or also in 7 and @, as can be seen by inserting p, £ and the 


proper values of +A into Eqs. (II, 5) and (II, 6). 


) 


One obtains for G™ and YWA/3 = (1/12) +h 

y , F[(2/3) + Al/T(2/ 
= (— . 

r(q /2 

for G’”’ and ~/A/3 = (5/12) +h 
Tf(4/3) + AY TU 
r{(l 

for G’ and VA/3 = (7/12) +h 


h 
¥ =(-) Sa 


r((3/2) + hj/T(3/2) 


r[(2/3) + h]/T(2/3) ( 
Ls 


and for G” and YA/3 = (11/12) +h 


, ri 3) + hy 4 3 : 

a ae. rf(3/2) + Ay 13/2)" ( 
(c) Another family of particular solutions is deter- 
mined by the condition y = 0 along = — © and along 
£= 1. The first condition requires that the particular 
solution be antisymmetric with respect to & = —o. 
The second condition requires that in Eq. (II, 2) WA 


be positive and that the first term in the parentheses 
should vanish. This leads to YA = (9 4) +. 3h. 
The behavior of the corresponding functions G along 
the sonic line is given in Eq. (II, 7a). 

(d) We need furthermore asymptotic expressions 
for the function G at high negative values of \. They 
are obtained from the differential equation (I, 7) for G 
by the method of R. E. Langer.? Some details of the 
derivation of the following formulas are given in refer- 
ence 4. For large negative values of X, the differential 
equation (I, 7) is integrated asymptotically by 


G~ (1 = ¢-9)""" ofp)" -e7 8" -T, 5f(V [al /3)-08)) 
(TI, 10a 

and 

ee | oe 3)! 24g)! 2 eit/6 7 v/3[t(V | A} /3)0(E)] 
(II, 10b) 


with v(t) = 2-tg-'[(—#)°? (8d, 11) 


2) +h) TO 2 


> | 
” (6) P| =, -—h,. =i | II, 9b 
a _ > . J 


l : 
g)* tip] — a eae II, 9d 
2 3 (9 '4)6" 


These expressions represent those exact solutions which 
agree with them asymptotically with respect to the 
value of the functions and the first derivative at & 0 


There the above expressions give 


and [v 7. os")? (2/8) ) II, 12b 


For — < 0 and large values of || the Bessel functions 
can be represented by their asymptotic series and thus 
one obtains other asymptotic expressions for the same 
functions G. One finds in either case 

1/12 a 3 V2e-/!! ‘Ye " 


2) Particular Solutions Which Fulfill the Boundary 
Conditions Along the Map of the Plate 


We seek particular solutions of Eq. (1, 2) that ful- 
fill the following conditions: 

(1) wy = O along the line ABC of Fig. 3. 

(2) y = 0 at infinity. 

(3) infinity is a branch point of the streamlines. 

(4 
(5) The characteristics EO does not map to infinity 


) point O isa singular point. 


Except for the character of the singularity at point 
O the problem is the same as that in reference 2. The 
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FLAT PLATE WITH AN 
treatment is also quite analogous except for some sim 
plifications in the representation. 

First one determines particular solutions of the form 
I, 6) which are suitable for the representation of the 
solution in the vicinity of the point B and at infinity. 
In the vicinity of point B, the boundary conditions 
y 0 along BC and AB are prescribed. 


¥ 


problem treated in Part II, le, and one obtains as a 


This 1s the 


possible set of particular solutions 


y= p’? t*.G{s [(9/4) + 3h]?} (II, 14 


For the representation of the solutions at infinity, the 
condition yY = O along BC restricts us to solutions 
which are antisymmetric with respect to the line & = 

condition 2 restricts us to negative values of 
V/\, and the condition that no singularity in y or in 
the higher derivatives propagates along EO leads us to 
the problem treated in Section II, lb. The values of 
VX given in the second part of Eq. (II, 3) are then 
applicable, and one obtains as a possible set of particu- 
lar solutions 


y= p 9 t4. Gls (9/4) + 3n]?} (11,15 


Condition 3 imposes a further restriction. At a large 
value of p, the one particular solution with the ab- 
solutely largest power of p will prevail. Now one can 
see that among the particular solutions (II, 15) only 
: “lé, (9/4)?] has no branch 


the expression p ”'-G'"’[é, 
If this particular 


point of the streamlines at infinity. 
solution is excluded, infinity will always be a branch 
point of the streamlines. 

The solution investigated in reference 2 is repre- 
sented along the line BC by the expression 


y = (-£ ve (2) [ ~ (2) II, 16a 
Po, Po 


Developing Eq. (II, 


r[(5 
V — - = i : 
g Los 73 
} 0 I ») 7 
on rl[(5/2 
or p < po and y —£ - ¥ (-)! os 
0 I (5/2 


for p < po. 


From these expressions one obtains as analytic continuation 


r((5/2) + » + h] 


ANGLE OF ATTACK SO 
here py is the value of p at point O. For the present 
purpose this expression may be generalized in the fol 


lowing manner 
— SS 
y = (—£)~*/2 (! \ 
po/ 


m and n being integral numbers. Since 


solutions obtained by chatiging either m or n» are not 
independent from each other. 

We must show that the analytic continuation of the 
expression fulfills the conditions given at the begin 
ning of the paragraph, and we must study the character 
of the singularities which arise at point O. The final 
aim is to determine the coefficients a which occur in the 
expressions (I, 9), to compute the change of the x 
coordinate caused by these particular solutions, and 
finally to determine their influence on the pressure dis 
tribution. 

That the analytic continuation of the expression 
II, 16) is free of singularities in the entire subsonic 
region can be shown with the method of reference 2 
i.e., by means of the particular solutions (I, 6) with 
negative values of \ and asymptotic expressions for 
large values of | A This part of the work will not be 
repeated here. The remaining computation, however, 
will be carried out since the numerical results are 


needed. 


16) with respect to p or p~', one obtains 


mn _ [= +." G 
¥ i= © r((5/2) + n]-Tih + 1 
r[(5/2) ++ h] j | (°) 7) 
for p < po and ye" = > (—)A ae + 3(1 +n m +-h x 
Pare . =e r((5/2) + n]-Tih + 1 ' | if 


for p> po. 


Then Eq. (II, 7a) yields the following representations of y”"’ 


along the sonic line 
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© 2/3) +mt+h ee ; peo Goat! Vi se Yi 
ae » ® §o1/3 p r((5/2) th + nj-T[(5/6) + m + h]-T5/3 
v = r » ( ) 
} 0 


r((5/2) + n]-T(h + 1)-T(5/6)-T[(5/3) + m + hh] 


2/3) +m+h (2 /o i c -Ti(7/eR) + dk BI ) 
22D (*) r((5/2) + n+ hj-Tl(7/6 m h\-T(4/3) 
/ 0 


r((5/2) + n|T(h + 1)-T(7/6)-T[(4/3) + m + v}h 


po 


for p < po, and 


‘ t6) bn Se SF ene ; , ae 

y””’ nt¢mt1 2 Joao (? , r((5/2) + n+ hj]-T[A1/6) + 2 — m+ Al-T(5/3 

at sso len C((5/2) + nJl(h + 1)-T(5/6)-T[(8/3) + n — m + | 

11/6 + » / ’ - ’ ‘ . ’ ‘ 
2/3 4 r((5/2) + 2 + h]-T[13/6) + n — m +h]- T(4/3) | 

a os 
= S Lad . - ' ’ (= ’ . -—/9 . 
h=0 \po r((5/2) + nj-Tih + 1)-1(7/6)-T[(7/3) + n — m + hj 
for p > po. 


This can be written as 


Par 5) dias lll foi3 C((5/6) + m]-T(5/3) _ fd D a p 
V (=) s—s “42 Sea ee ae +f + n,- + m,- +m, 
S \po | P(5/6)-T[(5/3) + m] 2 6 3 p 


for p < py, and 


a 11/6) +» ) . ht he: fee ‘5 
se 292i fi T((11/6) + 2 — mj-T(5/3) [5 11 8 
y “<2 Sr Ge -F +n,—+n—m,-+ 
S \p | r'(5/6)-T[(8/3) + n — m] 2 6 3 
eo rec tia a ae 5 a e 17 
p oj, T[(13/6) + n — m]-T(4/3) _] 5 1S. ( / p 
n Mm, — 2 a ee Ff oh, oe =e, rs wd | r (1h. 2g 
pO r(7/6)-T[(7/3) + n — m| 2 6 3 pr ab 
for p > pv. 


From these equations we shall derive representations for y in the vicinity of point O and thus study the singularity 
arising there. By means of the following relation for hypergeometric functions (see e.g., reference 8, p. 257). 


5° 


C(a)-T(b , l(a)-T(b — a) ee : 2 
; ‘F(a, b, c, 3) = : (—z)—“"- F(a, 1 +a—c¢,1+a—5),2-') + 
l(¢ r(c — a) 


r(b) l(a — bd) Pe 
(—2Z ’.F(b,1+6-— 6,1 +. f ome ar ie (TE. 2) 


~ ‘ > vail to. 5) 5) p 
one finds that the expressions -F + 2, — + mM, — + m, 
» . » 
po ys 18) >) po 


p 11/6 + n , ie 1] Q p l 
and ff + n, +n—m,-+n—™, 
po yA 6 3 Po 


which occur in Eqs. (II, 19) and (II, 20) fulfill the same hypergeometric differential equation. The same holds for 


2/3) +m - oy 
p | 2 i f p 
( ) F| ~ 8, - + Mm, - + m,( ) from (II, 19 
po 2 6 3 po 
p 11/6) + m n 5 13 - p l 
, bs 2 ‘ - 
and ( ) :F + t+ 4S = OM, ) from (II, 20) 
po 2 6 3 py 


As a consequence the expressions (II, 19) and (II, 20), if developed at p/p) = 1, will be represented by linear com- 


the expression 


binations of the same hypergeometric series. To establish that the inner and outer representation ‘“‘fit’’ at point O, 
one must merely show that the coefficients of corresponding hypergeometric series are the same. 
By means of the formula (reference 8, p. 252) 
: (ic) l(c —a—b) _ 
Pa, 6,.c, 2) = : 2 -F(a,b,1+a+6-—c¢c,1-—2) + 
lic — a)-Tle — 0) 


(ic) (a + 6b — c) eT 
A. ag ~*"* i «geo i tec ~e -4h1- Oo 1S 
'(a)-T(d) 


one obtains, using known relations for the T'-function, 
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. , 2 p 2/3) + m - . (5 31 ((5 6) + m)- r {i te ie n| r 5 5 Ss p 
y - ems ai agree =a +n, + m,- +n, 1 - 
S r(5/6)?-T[—(5/6) + m — n] 2 6 3 p 


r(7/6)-T1/6)-T[—(7/6) + m — n] 2 6 } p 
2/3 P(A 3)-T((7/3) + n] - p ie 7 l | p\t om 
&-2 l F{ —- + m n n, | it. 
r(7/6)P[(5/2) + n] pr \ 6 6 ) po) I 
and p < py and 
, li /t >) oT 7 , - 
} p {_. T(5/3)-T[(5/6) + m]-T[—(5/3) — n] ) 1 | 
y (— 2 = } + + nH ” 
S p { r(5/6)?-T[—(5/6) + m — n] 4 6 


S ny 1 +1 arg 2 (5/3)-T[(5/3) + n] a aS 5 2 p 
+ n,1l— + (- -2 re i — ‘i t= n, | 
3 p r(5/6)-T[(5/2) + n] p 5) 6 3 p 


ye 1'(4/3)-T (7/6) + m|-T[—(7/3) —n] _f5 13 lO p 
£9" wi PR — +2 — m,— + 82,1 


7 (7 /6)-T(1/6)-T[—(7/6) + m — n] 2 6 3 p 
r(4/3)-T((7/3) + n] 0 beat »\) 
n+l , 2/3 . . , 
— -£-2 re en 1 — -Fi —-— m,-, n, | I], 24) 
r(7/6)-T[(5/2) + n] p 6 ( 3 p/' 
for p > po. 
rhe coefficients of the hypergeometric expressions which correspond to each other in the representation for p < po 
and for p > py are the same. The factor (— which occurs in the second and fourth line of Eq. (11, 24) is caused 


by the fact that the particular solutions Y which can be used to represent the solutions in the vicinity of point O are 
either symmetric or antisymmetric with respect to the line 6 = 6, depending upon the power of [1 pp) | which 
occurs in front of the hypergeometric functions. 

For the determination of the a occurring in Eqs. (1, 9) we need the development of y”"" in the vicinity of point O 
in terms of the particular solutions (I, 6). Here the following difficulty arises: The particular solutions (I, 6) are 
defined in the p, £-system—.e., in a system which has point O as its origin. According to the definition of p given in 


Eq. (1, 5) we have along the sonic line 


pi po = (6 4)* and (p/ po) = 1— 6 @)- = (1 NV pp 
rhus in order to find the desired development in particular solutions (I, 6) we need a development of y'"" along the 
sonic line in terms of (1 — V p/po). 
We proceed in the following manner. Along the sonic line the function y”:" fulfills a hypergeometric differential 


equation. We make in this diTerential equation a transformation of the independent variable 


1 — V p/po = u il, Zo) 
Then u = Ois the point O. The expression y”’” represented as a power series development with respect to u has the 
desired form. The coefficients can then be determined by a recurrence relation. Along the sonic line, the expres 


sion for y”’” is given (except for a constant) in Eq. (II, 19) as 


2 


v= (pp 2) P(e /9) 4 n, (9 6) + m, (5/3) + m, p pol 
It fulfills the differential equation 
v"p(l — p) t+ W’[(1/3) — m — p38 — m+ n)] — ¥-(1/6)-[(111/6) + m 0 
By the transformation (II, 25), one obtains the differential equation 


d*y/du?)(—u® + 3u2 — 2u) + (dW /du)} —u? + 2u — 2(3 — m+ n)(u? — 2u + 1) + [(2/3) — 2m]} - 


With a hypothesis v= Dd a’u’* II, 26 


where v is a constant which will be determined subsequently, and the a,” are suitable constants, one obtains the fol- 


lowing recurrence relation for the a,”: 
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a,’ (k + v)[—2(k + v) — (10/3) — 2n] + 
k+v-+ (1/6) + (2/3)(n — m)/? — (103/36) — 4(n — m) — (4/3)(n — m)?! ~ 
a, ~ 2° [k + v — (5/3) ][k + v + (5/3) + 2(n — m)] = O (II, 27 





vie 
a; 1 1) 


The series starts for such values of v, for which the coefficient of the term a, vanishes for k = 0. This happens for 


vy = Oand vy = —(5'3) — n. The lowest powers for u which occur then in Eq. (II, 26) are naturally the same as 
the powers of (1 — p py) in Eq. (II, 23). For v = 0, the first line of Eq. (II, 23) yields directly 
o — ¢ym 3 oi3 P(5/3)-T[(5/6) + m]-Tl—(5/3) — 1] . 
Go =~) *=*24 ee ae BT (II, 28a 
S r(5/6)?-T[—(5/6) + m — n] 
For y = —(5 3) — n, we transform the expression which occurs in the second line of Eq. (II, 23), in the following 


manner 
(5/3) — n j 5/3 n 
(1 — p/ po) = (1 + V p/ po) “(1 — VY p/ po) 
and thus obtain 


3 13 T(5/3)-T[(5/3) + an] 


~/e — (5/3) n ) 
do 5/3 a ») 5/3 —_ ym a —— oe (II, ILb 
8 r(5/6)-T[(5/2) + n] 
The second term in Eq. (II, 19) is essentially given by 
. 2/3) 4+ aye pa si ‘ 
£-(p/ py) ”".F[(5/2) + n, (7/6) + m, (4/3) + m, p/po] (II, 29 


The corresponding expression in a development around point O will contain a factor —. Along the sonic line we have 


fms a i be 3 3) (5 po)! a o)! E.(p po)! 3( 05 p)' 3 (II, 30 
Thus we consider now the expression y = (p/p)? *™-F[(5/2) + n, (7/6) + m, (4/3) + m, p/po] (II, 31 
It fulfills the differential equation 
y"(1 — p)-p + y’'[(2/3) — m — (4 +0 — m)-p] + v(—5/6)[(13/6) + n — m] = 0 
which, by the transformation Eq. (II, 25), assumes the form 
(d?y ‘du*)(— u? + 3u? — 2u) + (dy/du)} — (2u — u*®) + (4/3) — 2m + (—u? + 2u — 1)[S8 + 2(n — m)]! + 
y(—10/3)[(13/6) + (n — m)(u — 1)] = 0 


With a hypothesis 
y= D> du" tT" (IT, 32 
one obtains the recurrence relation 


b“(k + w)[—2(k + mw) — 2n — (14/3)] + 
) “1 3[(k + uw) + (5/6) + (2/3)(n — m)]? — (103/36) — 4(n — m) — (4/3)(n — m)*} — 








b, — o*(kR+ mu — (1/3) ]}k + w + (7/3) + 2m — m)[k + w — (1/3)]f = 0 (II, 33 
The values of u are found as uy = 0 and uw = —n — (7/3). If one replaces in Eq. (II, 23) & by — by means of Eq 


(II, 50), and also inserts the expression for Eq. (II, 31) just found together with Eq. (II, 30), into Eq. (II, 29)—then 
the results have the same form. The initial coefficients 6 are obtained from Eq. (II, 23) as 


ma. (3\ 93 T(4/3)-P((7/6) + m]-T[-(7/3) — n] 
by = (-) 2°". —— — (II, 34a 
S r(7/6)-T(1/6)-T[—(7/6) + m — n] 


éa oe + aa 2/3) — » '(4/3)-T[(7/3) + n] 
Ra = (—) oP ise Ali = (II, 34b 
S (7/6) -T((5/2) + n] 


Since the representations for y’”"’” fulfills for p > pp the same differential equation as for p < pp, an investigation of the 
. . . ~ . “* a 
expression occurring in Eq. (II, 20) would have given the same results. Since 1 = V p/p) we have as representa- 


tion for y””” in the vicinity of the sonic line, if 0 < 4 


yy = >> ay’ (p po) 4 (p/ po) wie dha = Paha ~ (p/ po) “4 
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PLAT PLATE With 


The particular solutions of the form (I, 6) that are 
suitable to represent the solution in the vicinity of 
One must 
only replace the quantities p and & The 
powers of occurring in the last equation corveaponil to 
the values given in Eqs. (II, 3) and (II, The be- 
havior of the corresponding functions G along the sonic 
Thus it is 


point O have been discussed in Part II, 1b. 
by p and é 


line is expressed in Eqs. (II, 7) and (II, 8). 
possible to find from Eq. (II, 35) the superposition of 
particular solutions (I, 7) by 


which the solutions 


y”” are represented in the vicinity of point O. 
Here the following fact is remarkable. The terms 
: 5/3) — : ~(7 , . 

with a; "and with b, "in Eq. (II, 


have the same power of f and belong therefore to the 
same particular solutions (I, 6). The coefficient of 
these particular solutions are determined independently 
by the two different equations. It appears desirable 


to show that the results are the same. We proceed in 


the following manner: Using the definitions (I, 6) 
and the expression for G, Eqs. (II, 7) and (II, Sa), one 


finds that for symmetric as well as for antisymmetric 


particular solutions 


ve »/6 k/2 = — YX 2/3 P(2/3)-T[(7 3) + k| 
baa <an r(43)-P[(5,3) + &] 
(II, 36) 
Hence it follows that 
h 7/3 ’ - —a; 3 ’ 9 2/3 x 
r(2/3)-T((7/3) + 2 — k] se 
\ iT. 34 
r(4/3)-T[(5/3) +n — k] 
This relation is fulfilled for ao "and by 
Eqs. (II, 28b) and (II, 33b). Furthermore, one can 
show without difficulty, that the recurrence relation 
II, 33) for b,~ originates from that for the 


/d 


a, °* ss Eq. (IT, 2 
Eq. (II, 37) 


w= —(7/3) — 1. 


, by inserting into these relations 
n and 


and the values »y = —(5/3) —- 


6 


For the future determination of the functions ¥~” °, 
and W~'' © we list here the first terms of the 


development of some of the expressions y” 
sonic 


reer lOO” 


higher order aie’ (II, 38a) 
4 


= 1/3 
westeel OCs 
9-V3 L 3 Po) 

higher order one} (ea: 

-_ 11 6 9 
eel -C> 

ays Aes po 378 
() 

po 


Furthermore, we 
valid along the sonic line which originate from the defi- 


VY a/s 
m 


along the 


line. 


3Sb) 


higher order terms | (II, 3Sc) 


mention the following relations, 


nition (I, 6) and the formulas (II, 7) and (II, 8) 


AN 
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YV-se = + (?).2 (? yo II, 39a) 
py 
¥-11/6 = We (PR) (II, 39b) 
¥-13 = g-10(2) 11, 39e) 
p 


tae e (JC) cram 


To determine the final expression for the functions 
v°** w-*?) etc., we need the 
the plate. 
According to this equation, a certain linear combina- 


ation of functions y leads to the same linear combina- 


x-coordinates along 
Their evaluation is based upon Eq. (I, 4a). 


tion for the corresponding function x and y. 


Thus we determine first the contributions to x 
caused by the functions y°°, y"' etc. They may be 
nin" 


0.0 el, i 


denoted by x*°, x 
to Eq. (II, 16), we have along the side of the plate 


vy nin 


or in general by x According 


. 3/2 ) + > , 
(—€) (p/p) [1 + (p po)] 
II, 40 


or in terms of 7 and 6 


- 36 — 
= > 3 v( 2 
2(-—n [(3,/2)0.]> 


Hence with Eq. (1, 4a 


Ox" 13 [3\ (36 
O 7 r | , ») ») x 
3/2) + 3) 7 r 3 » 2 
—" at dé gh, 
(j 3, 2)60)° ) | ( 3 2)60)° ) | 


It is practical to introduce a new variable by 


(30,/2)°'* II, 41 


7 * 


The last equation is then written as 


and thus, denoting a variable of integration by u 


/ 


_ 
] 


eines | x se " du II, 42 


Here the constant of integration is chosen such that 
the stagnation point and, in the present approximation, 
the leading edge lies at x = 0. The contribution to 
the plate length of these expressions is found, 


(equivalent to 


if one 
evaluates the integral from 7 = — 
©) to 7 = 0(n = O gives the point on the pressure 
the trailing edge) 


1 = 
side where the sonic line starts——1.e., 


Let us introduce 
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0 
[%" = fi us? rowdy +)" ~- "de Gt, 43) 


By the transformation u* = (v, 1 — wv) this expression 
assumes the form of a B-integral and one finds 


as -(1/3)}T[(5/6) + nj]-T(5/3)/T[(5/2) + n}} 
(II, 44 


In particular we have 


9 = —(1/3)-(P(5/6)-T(5/3)/T(5/2)] (II, 45a) 
J! = (1/3). (II, 45b 
[22 = (11/63)-/%° (II, 45c) 


After these preparations, one can find the expressions 
for WV w** ete. w-”® fulfills the boundary 
condition at the surface of the plate and possesses as 
strongest singularity at point O the solution ~_; ¢% 
Thus, from Eqs. (II, 38a) and 


with a coefficient one. 
(II, 39a) 


y~*® = —(27/8)-V3-y" (II, 46) 
The constant / is the plate length obtained for this solu- 
tion, it is found by means of Eq. (II, 45a) as 


9 /~ 1(5/6) -T(5/3) 
| = -V3 z= 4 
S rt5/ 2) 


Furthermore, one needs for the evaluation of Eq. 
5/6 


(II, 19) the quantity a_j/3 which is found from 


Eqs. (II, 46), (II, 38a), and (II, 39c) as 


a1/3 = —19/16 (II, 48) 
1/3 . ° ‘s 

WV has as strongest singularity at point O the ex- 

pression ¥_4/3 with a coefficient one, and, furthermore, 

its contribution to the plate length is zero. Hence 

from Eqs. (II, 3Sb), (II, 39d), and (IT, 45) 


v4 — 1(135-V'3) S][y'! — (1/3)W"] CII, 49) 


Finally, the solution W~''’° is characterized by its 
strongest singularity given by ~_1,,/4 with a coefficient 
one, by the absence of the singularity Y_4/3; and by 
the fact that its contribution to the over-all length of 


the plate is zero. 
ue  27°63-V3.... 143 1 
yp 11/€ eS |v _ yr! J - | 


(II, 50) 


Then the corresponding values of x (along the plate) 
are given by replacing the functions y by the expres- 
The results are listed subsequently 
“1/8” actually all the 


sions (II, 42). 
3 
and x 


. . —4 
(the expressions for x 
higher x’s, can be evaluated by integrations by parts): 


Because of their simplicity we mention furthermore the 
following expressions which are needed in evaluating 


Eq (I, 26) 


x 13 » ‘ 11/6 
5/6 — (— n = = 
(d dn)(x ) 5 (d dn)(x 
63 2 (—n) 1] 
(—9% — . + = (Ti, 52 
16 21 1+ (-—7)?’ 135 


The form of the first of these expressions is not surpris- 
ing after the interpretation which has been given to 


vw? in part I, 6. There we found that for a change of 
the angle of attack by Aa, the solution for the flow 


field is given by 
56 | /¢ 13 1/2 
(L/1)-[W " — (3/2)- pr -Aal 


On the basis of the considerations of part I, 9, the change 


of the pressure distribution is then given by 


An = (3/2)p0- *. Aa [x 9 (da dn)(x si 
or with Eq. (II, 5la) and py (9 4)007 
An = (2/3)-n- (Aq 6, (II, 53 


On the other hand the pressure distribution over a flat 
plate can, in principle, be written in the form 


) 


n = (xX) + (3500/2) 


If one changes the angle of attack # by Aa then one 
obtains 


[An /(30)/2)7*] = a(x)-(2/3)(Aa/O 


This is Eq. (II, 53). 


3) Particular Solutions Which Have y = 0 Along OQ and 
Are Singular at O or Q 
We seek particular solutions which are symmetric 
with respect to the line § = —o and have y 0) 
Naturally one will have along OQ a change 
Along the line § = — © from Q to 
The solutions will be singu- 


along OQ. 
of the sign of Wo. 
n = — ©, We will be zero. 
lar either at point QO or at point O. The second possi- 
bility is a generalization which is not needed in the 
present applications. 
We start with an assumed form of y along = = 

for p > pe 

ihe * (IL, 4 


v = (B/Be) [1 — (p Be) "| 


here k is a positive or negative integral number. 
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This expression will be continued analytically for p > pe. To obtain this solution we choose in Eq. 
through the subsonic region and it will be shown that (II, 54) k = 0. Then one finds by a comparison of 
the continuation fulfills the desired boundary condi- Eq. (1, 14) with Eq. (II, 55), using Eq. (I, Se 
tions. 

1 . _ "7 =-/= > - 6—4/3/6—-1/3 = 1/2 (II, 57) 

Along the line § = — - and for p/fg > 1, Eq. (II, 54) : , 
is developed in the form To show the convergence of the expression through the 

ayn r((1/2) +h+e] “ subsonic region, we use the particular solutions (I, 6) 
dt owe . é e ee ae a The . ‘ ae | at 
Zo no [(1/2) + k]-P(h + 1 = for negative values of . lhe two terms on the right 
(II. 55) side of Eq. (II, 56) will be treated separately. For 
negative values of \ the particular solutions (I, 6) have 
By means of the particular solutions given in Eq. along the sonic line, the form 


I, Sc) and with the aid of Eqs. (II, 9a) one finds as the 
1 


: a dibs _—1/12 1a ‘i = 2 

representation along the sonic line p cos [(V || /3)-log p] or 6 - 
: 3 op 31 (V 1) 3)-log 3 
y = (1/2)(B/Ba) 34 O'/3. FI(1/2) + sin | A} /3) -log 3B] 

= » a po ye —1] = 9-1/3 . ° , . > 
k, 5/6, 2/3, (—B/Be) | — €-2 > F[(1/2) + Accordingly, we develop the expression p'/'*-y taken 
, w/e a ate 172 Ke) : ; ay : 

k, 7/6, 1/3, (—B/Bq)"]§_ (II, 56 from Eq. (II, 56) with respect to cos [(V |) /3) log 5] 


As a series representation the expression converges for and sin [(V |A| /3) log p] and then replace each term 


3 > pa only, but as a hypergeometric function the by the corresponding particular solution (I, 6). As- 


representation along the entire sonic line is obvious. sume that p ey is represented along the sonic line in 
From here on we shall proceed in the following man- the following form 

ner: By a development of the expression found along 

the sonic line with respect to particular solutions (I, 6) p ey = [ gi(u)-cos(uw)du + [ go(u)-sin (uw) -du 
for negative values of A, we form a representation for si vali II. 58 
the subsonic region which shows that no singularities os ie 
will occur except at the line — = — ©. Since this pro- where w = log p, where u replaces log f and where the g 
cedure differs somewhat from the one applied in refer- are functions which arise by the development. If 
ence 2, some of its details will be shown. In a further p''*-y does not fulfill the conditions necessary for the 


step we develop the solutions for p < fg in terms of the convergence of the Fourier integral then a finite num- 


particular solutions (1, Sb) and thus show that along 
OQ the boundary condition y = 0 is fulfilled. beforehand. In the complex j-plane, p' '"y has the 


ber of suitable particular solutions must be subtracted 


But first, before we carry out these investigations singular points —o, 0, and —s since it is repre 
whose sole purpose is to justify the assumption, Eq. sented by hypergeometric series Eq. (II, 56). There 
II, 54), we derive the expression 6_4/3/6—1,/3 needed for fore the representation (II, 58) converges in the complex 
the evaluations of the first part. The function to be w-plane in a strip —7 </mw < 7. This gives an esti- 
represented is shown in Eq. (I, 14); it is characterized mation for the g. For |/mw) < mr — e the expression 

. . > . a = 1/2 - . 
by its behavior along § = — ~—+.e., Y ~ (pf — pag (II, 5S) can be written as 
)2i(u) [cos (u Rew)-cosh (u-Jmw) + 7 sin (u Rew)- sinh (u Imw)] + 

go(u) {sin (u-Rew)-cosh (u Imw) — i cos (u Rew)-sinh (u Imw)) {du 


This expression would diverge if the coefficients of cos (wu Rew) and of sin (« Rew) were unbounded for large values of 
u. Using the approximation for hyperbolic function of a large argument, one finds 


V g7(u) + g27(u)-e <M gi(u) < M-e g(u) < JM-e 


and (II, 13), that a function G(£, \) which assumes the value of G(0, A) along the 


Now it follows from Eqs. (II, 12) 


sonic line, and for which dG dé is zero, is represented for large negative values of \ and for — < 0 by an expression 


"/ ‘ 1/12 d 5 
G(é, A) B G(O, A)-const.-!Ar “¢” 
or, using AX = —9ux" 
"/e ‘ ° 1/12 ( 
[G(é, —9u*) ‘G(O, —9u?)] & const.- |r -e 


rhe analytic continuation of the expression y from the sonic line through the subsonie region will be free of singular 
ities, if 
daft 3 G(é, — 9u? Git, — 9u? 


ie 1 12 . s 
Y= > o,/N) Cos (uw) = — + go(u)-sin (uw) = —Idu < 
0 G(O0, — Vu 7(Q, — Yu- 


p _ [ M-|\|~'''* const. -e oe ea 
0 
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converges. This is the case for a — e€ where e > O but can be taken arbitrarily small. According to Eq 


Vv 
(II, 11), v(&) = wecorresponds to & = —«. This completes the proof of the convergence. For the contribution dye 
to y; along the sonic line one proceeds analogously. 

To show that the analytic continuation of expression (II, 56) fulfills the condition Y = 0 along OQ, we use the de. 


velopment for 


— °+* pl/3) — k]-1(2/3) {(:) - (°) - (?) ou ( a) £ 
"= \(2) r[(1/6) — k]-1T(5/6) 7 2} NG is i Be 
p\'* Pik — (1/3)]-T(2/3) (?) (‘) (:) " ( ai\= 
2 PL = G/@)/-FIG/2) + HY rf Nel \3a) ~ A pod If > 
ye hime) ++ O)-a()+m(- 2) 
hae eee lag d 
eae et 
Da r(—5/6)T[(1/2) = - s) 


The powers of p/ jg occurring here are those of the particular solutions which have y = 0 along OQ and are regular 
along the Mach wave EO, Eq. (I, Sb). The continuation of the second and fourth term of the last equation to ( IQ is 
obvious. The first and third terms contain the same powers of p jy, it can be shown by a straightforward computa- 








tion that these terms are represented by the same particular solutions (I, Sb). 


(4) Particular Solutions Which Have y = 0 Along OQ and Are Singular at Infinity 


The particular solutions in question are obtained by a generalization of the assumption, Eq. (II, 54). Denoting 
momentarily the expression (II, 54) evaluated for k = 0, by y, they are given along § = — © by 


= C) J @ ’ be & — 
7 oO S@) i _ (7) | od . 7 


The choice of the constant of integration is an important part of this analysis; it will be discussed later. 

If one writes the last expressions as a series, then the powers of (j/fg)~! which occur in it will be the same as in 
Eq. (II, , but the coefficients will be changed by the integration, depending upon the powers of j/jg occurring in 
the various terms. Thus, the analytic continuation of this function (II, 59) is accomplished by tee same particular 
solutions (I, Sc) as previously and the representation along the sonic line is obtained by carrying out in Eq. (II, 5 
the operation corresponding to Eq. (II, 59). Naturally, the operation can also be carried out on the individual 
terms of the powers series representations of the hypergeometric functions occurring in Eq. (II, 56). Here the con- 
stant of integration which will occur, will cause a term with a power (j/jg)" to be present which was not contained 
in the original series (II, 56). This term will represent the desired singularity at infinity; it is the main task to de- 
termine its coefficients. 

In the previous analysis, it was important that for the hypergeometric expressions of Eq. (II, 56) that describe 
the behavior along the sonic line, there exists a relation between the representation at p jg = 0 and at p/fe = 
This relation can be derived by means of Barnes’ contour integral for the hypergeometric function. Here Barnes 
representation will be useful for the determination of the constant of integration. 

According to Barnes’ formula (see, for instance, reference 8, p. 256) we have 


(a) -T( b) ‘ l { l(a + s)T(b + s)-T(—s) 7, 
(—3) ds 


a, b, C¢,2) = 
(c + s) 


I'(c) 2m . 


where the path of integration runs essentially along the imaginary s-axis with such loops, if necessary, that the poles 
of I'(—s) lie to the right and the poles of (a + s) and I'(b + s) lie to the left of the path. The behavior of the in- 
tegrand at infinity of the complex s-plane is such that the path of integration can be closed at will through either the 
right or left half plane. 

We obtain for the hypergeometric functions occurring in Eq. (II, 56) 


~ 9 = —1 9/9 4-9 x . 9) I 5 (x > 7 —¢) _ s 
Foo (- °) |- (2/3) J [ r[(1/2) + s]P[(5/6) + s]T\ (7) ds (II, 60a 
26 3 Ba T(1/2)-1T(5/6) 2rt J -;i. r[(2/3) + s] By 


and 


~ -\-—! 1/1/92 Lio mn 9 -1.T1(7/e@ elca —<) - 
F[ie »(- * |- r(1/3) ff I'((1/2) + s]-T{(7/6) + s]-T(-: (7) ds (11.60b 
2°6 3 Be I'(1/2)-1(7/6) 2mi J -i. r[(1/3) + s] Be 


If one closes the path of integration in the right half plane, then the calculus of residues gives the series representation 
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in powers of (6/pg) 
(p/ pq): 
s] and r{(5/6) + s] or T[(7/6) + s]. 


we obtain for y along the sonic line 


AN ANGLE OF 


In the first case, the residues are formed at the poles of ['(—s); in the second case, at the poles of T'|(1/2 
If we carry out in Eq. (II, 60) the operation indicated in Eq. (II, 5g) then 
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| if one closes it in the left half plane then one obtains the series representation in powers of 


») 


, 1 F(2/3)-(B/pe)~ f r((1/2) + s]-P[(5/6) + s]-1 
/ ee z/s 
é 2. P py ad 
. 2at V(1/2)-T(5/6 : r((2.3) + s]-[s + 1/3) +h 
_ oa 1 TU/3)-(@/BQ) { r[(1/2) + s]-P[(7/6) + s]-M(-s) _ , 
."< ; , " = e DP pe d 
° 2ni -1(1/2)-T J =i r[(1/3) + s] [s + (1/3) + A] 


The integrands have besides the poles which occurred in the original expressions (II, 60), a pole at s 
The residues at the original poles are changed because of the factor [s + (1/3 
The contribution of these poles give the integrated terms of the hypergeometric series 


because of the integration. 


pre ceeding in powers of either p Pe OT (p/pa 


We choose the path of integration in such a manner that point s = (—1/5 
lies outside of the contour which leads to the representation in powers of p jg 
tion contains solely the integrated terms of the hypergeometric series. 


9 


+ h] occurring in the denominator 


h lies to its right. Then this pole 
i.e., around point O, the representa 


This shows immediately that the analytic 


continuation of this representation will give y = 0 along OQ. 
The representation for large values of p jg will contain an additional term caused by the residue at point s 


—1/3 —h. Itis found to be 
2/3 r'(2/3) r((1/6) — h]-T{/2) — A]-P[/3) + A). (?) 
y= -2 eo : : 
r(1/2)-T(5/6 r{(1/3) — hj Bo 
«et (1/3) Td 6) — hj-T[(5/6) — A]-T[1/3) + h] (*) 
slic Re »)-T(7 6) r —h Pe 


The next question concerns the analytic continua- 
tion of this term through the subsonic region. Ob- 
viously, the analytical continuation is given by a par- 
ticular solution of the form (I, 6) with VA/3 = 1/12 + 
h; in general it will be represented by a linear combina- 
tion of the symmetrical and antisymmetrical expres- 
sion y. 
symmetric particular solutions can be found by means 
One finds for the coef- 


The coefficient of the symmetrical and anti- 


of Eqs. (II, 5) and (II, 6). 
ficient of the symmetrical portion 


ir [(1/3) + h1-T(1/2)/P[(5/6) + A]; 


and for the coefficient of the antisymmetrical portion 
9 Jo 1/2fp 9 9 ’ 9 { 
—(2/V 3)-r’"jT[(1/2) — 3h]/T(— 3h);5 


The antisymmetric coefficient is zero, if 3h is a non- 
negative integral number. Then the solutions are en- 
These are 
the particular solutions used in Eq. (1, 10b). If (1/2 — 


3h) is a non-positive integral number, the ratio of the 


tirely symmetric with respect to § = —o@. 


antisymmetric coefficient and the symmetric coefficient 
is infinite—this gives the particular solutions (I, 10a). 

The present approach defined first the function y 
along the line § = — © by Eq. (II, 59) and then used 
the idea of analytic continuation. One can also carry 
out the operation indicated in Eq. (II, 59) in the entire 
D, E-plane with an additional condition regarding the 
behavior of the constant of integration in its depend- 
ence upon £; one thus finds a solution of the original 
differential equation for y. Since the function y de- 
fined by Eq. (II, 59) is regular along the Mach wave 
OE, the function originating by the operation Eq. 
(II, 59) will also be regular along this Mach wave. This 
condition applies also to the term with (j/jg)” and thus 


determines independently the ratio of the coefficients 
of the antisymmetric and symmetric portions. 

The idea just presented leads also to an independent 
determination of the coefficient of the symmetrical par 
ticular solution, at least in a heuristic manner. The 
solution gives Y = 0 along OQ. This determines the 
constant of integration in Eq. (II, 59); one obtains 


for large values of p/je 


= ee ee 
r* tf , i} (7) ‘ 


The integral which occurs here can be transformed into 
a B-integral. 
Thus 

Y = (p/pg) T(1/2)-Tl/3 


+h) T((56) + Al} 


Remembering that along — = — on the contribution 
of the antisymmetric particular solution to y is zero, 
we find immediately the symmetrical coefficient given 
above. 

The procedure just described gives the necessary re 
sults quite rapidly, but the author feels that the pre 
vious method is more suitable to derive these results 
rigorously. 

For the evaluation of the results of the first part we 
need the function (II, 59) for h = 0. 

Here one finds for the representation at infinity 


r(1/3)-T(1/2) 5 f/p\'* — /3\/5\-*” 
-* asa a . - 
(5/6 Pa ‘ PQ 


or with the notation of Eq. (I, 10b) 
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B 1 3° = —3-T(5 6) (rl 3)-TQ1 2)] 
B_4/3° = —3-T(5/6)/(8F(1/3)-T(1/2)]. 
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Transient Temperature Distributions... 


(Continued from page 836) 


body and the energy received from the surroundings by 
radiation is assumed equal to ¢eCgs7‘, which is in 
most cases sufficiently accurate for the present purpose. 
In cases where the thermal stresses are nearly equal to 
their asymptotic value for g > @, it is only necessary 
to replace in the previous formulas the fluid friction 
temperature 7’, by the wall temperature 7), as calcu- 
lated from Eq. (39). In other cases—i.e., if g < 6—the 
maximum thermal stresses occur before the skin has 
reached its equilibrium temperature and consequently 
at a time when the temperature level in the whole struc- 
ture is lower than for the asymptotic case. Since the 
heat radiated from the surface is proportional to the 
fourth power of the absolute temperature, the radiation 
effects will be less at lower temperature levels and con- 
sequently they will also be less for values of g < 6 than 


for the asymptotic case g—> ©. Thus an upper limit of 


the radiation effects on thermal stresses can be found as 
regards the effective temperature difference 7 — 7°; in 
the formulas for the stresses. 


The influence of radiation on the parameter g can be 
determined approximately by using an effective heat- 
transfer coefficient 


Heft, = Aeon + Grad (40) 


where the first and second term refer to convection and 
radiation respectively. The second term in Eq. (40) is 
equal to the left-hand side of Eq. (39) divided by 7, — 
T». 

In the case treated in reference 1 the material of the 
wing was steel and the flight data were JJ = 3.1 and H 


= 50,000 ft. Assuming the emissivity of steel as « = 


0.80, the upper limit of the error due to neglecting radia- 
tion would be 10 per cent as regards the effective tem- 
perature difference 7, — 7; and ag. would be about 14 
per cent higher than a,,,, Hence in this case the ther- 
mal stresses would not be materially effected by taking 
radiation into account. 
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The Load Carrying Capacity of Initially Flat 
Circular Steel Plates Under Reversed Loading* 


A 


R. M. Haythornthwaite and E. T. Onatt 

Assistant Professor of Engineering, and Research Associate, 
Respectively, Division of Applied Mathematics, Browr 
University, Providence, R.! 

May 17, 1955 


SUMMARY 


The testing of a simply supported, centrally loaded mild steel plate, 20 
in. diameter, is described The plate was loaded well into the plastic range 
ind then reloaded in the reverse sense Under reversed loading, temporary 
plastic instability was observed, as predicted by applying a new large de 
flection theory for plates composed of rigid plastic material Earlier tests 
have shown that the limit load obtained by neglecting changes in geometry 
is often of little physical significance for a monotonically increasing load 
These tests show that, in reversed loading, it becomes a measure of the mini 


mum load carrying capacity available during a phase of plastic instabifity 


INTRODUCTION 


I OADING TESTS Of structural steel frameworks!» ? 


that good estimates of their load carrying capacity can be 


have shown 
obtained by the methods of limit analysis In the theory, a 
rigid-plastic material is usually assumed and the load carrying 
capacity is estimated as the load at which a structure composed of 
this material would begin to deform It can be shown that, if 
the material is perfectly plastic—i.e., nonstrain-hardening, and if 
the accompanying change in geometry is disregarded, plastic flow 
continues under constant load.6 If, however, changes in ge 
ometry and strain-hardening are considered, then in general plas 
tic flow can continue only under either increasing or decreasing 
load. Hence the methods of limit analysis may require modifi 
cation before they can be applied successfully to other types of 
structure 

In the case of circular plates, recent tests? ® % have shown that 
limit analysis based on the initial geometry and neglecting strain 
hardening leads to estimates of load carrying capacity that are 
far too low. It has been suggested that change in geometry is 
probably the significant factor To investigate this point theo 
retically, estimates of the load carrying capacity of the plate in 
the deformed state must be obtained. For simply supported 
plates loaded centrally the shape remains close to the conical form 
suggested by the incipient (zero deflection) velocity field 
Consequently, an approximate load-deflection curve was found 
issuming the shape remains exactly conical. Each point on this 
curve represents an upper bound to the load carrying capacity of 
i conical shell composed of rigid ideally plastic material. This 
inalysis has enabled the writers to predict with reasonable ac 
curacy the load deflection curves observed for steel plates over a 
wide range of the thickness/diameter ratio.!! 

The theory leads also to important conclusions concerning be 
havior under reversed loading. These will be discussed in the 

* The results presented in this paper were obtained in the course of re 
search sponsored by the Office of Ordnance Research under Contract No 
DA-19-020-ORD-3172, Project No. TB2-0001 (1086 
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light of new data obtained from tests of a plate supported freely 
near the outer edge and carrying a central load 


THE R1iGip-PLASTIC PLATE IN REVERSE LOADING 


Fig. (1) shows the theoretical load-deflection r« lationship for 
various stages in the loading and unloading of a rigid-plastic 
plate 

When the load is smaller than a certain yield load P,, the plat 
remains rigid. Above P,, the theory for large deformations 
leads to the locus GAB, the right-hand branch of which will be 
followed The right-hand path is chosen to satisfy the condition 
that the plastic work done in the plate should be positive BC 
is a typical unloading path 

If loading is continued in the reverse sense after C is reach 
the plate will remain rigid until Pp, at which the plate will again 
deform plastically, but in the opposite sense Due to assumed 
symmetry of the yield condition, Pe = —Pp If the conical 
shape is again assumed, the locus obtained is DEF, which is 
symmetrical with BAG. When strain is controlled during a test, 
the locus DEF will be followed, implying a decrease in the applied 
load. On the other hand, when the load is maintained, as under 
dead weight loading, the locus jumps directly from D to F, as in 
dicated by the dotted line in Fig l Thus the decrease in load 
along DE indicates the possibility of a temporary plastic instabil 


ity. 
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Fig. 2. Mild steel plate under 
First cycle (load range 





reversing - load 
t 30,000 Ib) 








15 fio -3 | 








Fig. 3. Mild steel plate under 
reversing load - 
second cycle 
(load range * 50,000 Ib) 


Fig. 1 is drawn for a material that can be considered to have an 
indefinitely large modulus of elasticity. If the modulus is finite, 
as in the case of all metals, the load-deflection cycle will be 
modified by the elastic deflection; but it may be expected that 
the general features of Fig. 1 will be preserved, since it is drawn 
for what can be regarded as a limiting case. Thus the reversed 
bending test provides what may prove to be a unique opportunity 
for observing temporary plastic instability under easily controlled 
conditions. In addition, the exact symmetry of the loci GAB 
and FED enables the role of strain-hardening, the Bauschinger 


effect, etc., to be assessed. 


THE TESTING OF A MILD STEEL PLATE 


The theoretical considerations described above have prompted 
a new program of experiments. One test that serves to illustrate 
several important points will be described briefly. A complete 
report of the program will be made elsewhere when all the data 
are available. 

The plate was cut from rolled steel '/. in. thick. For testing, 
it was freely supported on a steel ring, inner radius 10 in., and 
loaded at the center through a 2-in. diameter hard steel punch. 
The test was carried out in a 120,000-lb. capacity universal test- 
ing machine fitted with continuously variable control of the dis- 
placement rate. The central deflection of the underside of the 
plate relative to the lower pattern of the machine was measured 
by means of a 0-1-in. dial gauge reading to 0.001 in 

The load-deflection curves for the first two cycles of loading are 
shown in Figs. 2 and 8. For the first cycle, Fig. 2, the plate was 
loaded monotonically to 30,000 Ibs. This load was then re- 
moved, the plate turned over, and the same load applied in the 


* During both loadings, the deflection after 


opposite direction 
each increment of load was recorded as soon as the time rate of 
deflection had fallen to 0.001 in. per min. or less. For the 
second cycle, Fig. 3, the plate was loaded to 50,000 Ibs. in each 
direction, a steady rate of displacement of approximately 0.04 in 
per min. being maintained in each case. The zero deflection 
datum for Fig. 3 was fixed during the second half-cycle by ob- 
serving the deflection reading at which the plate appeared to be 
most nearly flat, as judged with a straight edge. 

Tensile control tests were performed on standard flat test pieces 
and the average tensile yield stress was found to be 32,000 Ibs 
per sq. in. On the basis of this yield stress, the theoretical 
load-deflection loci described in the previous section were com- 
puted for the known dimensions of the test plate. They are in- 


dicated by the dotted lines in Figs. 2 and 3. 


DISCUSSION 


As has been anticipated, the finite elastic modulus of steel 
leads to substantial modification of the theoretical load-deflec 
tion cycle; however, the general features indicated by the analy 
sis for rigid-plastic material can still be traced. Although the 
corners are considerably rounded off in a manner that recalls 
the buckling of elastic shells, quite definite instability can be ob- 
served in reverse bending. The trend of the two cycles shown 
suggests that the instability will be inhibited completely for 
cycles with lower maximum loads, but its presence at high loads 
is beyond question. 

For the quadrants in which the plate was stabie, the general 
shape of the experimental curve is well represented by the theory 
This provides further confirmation that the dominant physical 
factors have been properly introduced in the theory 

In the first cycle, Fig. 2, the load at which the ordinate i 
crossed is of particular interest because it is less than expected 
This is probably yet another example of the Bauschinger effect 


* The plate was supported on a ring 20 in. in diameter internally The 
plate itself was 21 in. in diameter, so it overlapped the edge of the ring by 1/2 
in. When the plate was turned over after the first half cycle, it would be 
supported at its outer diameter This would slightly increase the rate of de- 
flection in the elastic range, but the ultimate carrying capacity would be un 
affected In any case, the plate would be supported again at 20 in. diameter 


as soon as the slope at the edges changed sign 
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READERS’ 


The discrepancy disappears after the higher load has been ap- 
plied, Fig. 3. The Bauschinger effect may have been cancelled 
by further strain-hardening; on the other hand, it is likely that 
the plate was never truly flat at any stage following the applica- 
tion of the higher load. This would result in a slight increase of 
strength without the aid of changes in the properties of the steel 

It is remarkable that the large deflection imposed during the 
second cycle appears to produce an inconsequential increase in 
strength as the ordinate is crossed. The load carrying capacity 
remains very nearly the limit load obtained by neglecting changes 
in geometry (the apex of each locus Thus this particular limit 
load, which is of little significance for monotonically increasing 
load, remains of continuing significance for reversed loading. It 
becomes a measure of the minimum load carrying capacity avail 


able during a phase of plastic instability 
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On the Decay of a Flow Nonuniformity 
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y ips DECAY MECHANISM of a flow nonuniformity is a topic of 
considerable complexity, but a simple calculation is capable 
of yielding interesting information on the flow changes resulting 
from such decay 

Consider the case of a steady one-dimensional nonuniform 
compressible flow in a constant area duct with negligible wall 
friction. Let it be assumed that at each section the flow is non 
uniform in its velocity but that the pressure and temperature are 
uniform 

Application of the three conservation equations to this case 


produces the relations 


or y—-1 a, 
M*\ 1 + i) 


(1+yiW 


= const (1 


FORUM S69 


where by definition 
V/ u/V ygRI n = (u u? 
and with 
t= ] ” ¥ udni s- = ] mf ul dm 
In these, u is the local instantaneous velocity, 7 and p the stati 
temperature and pressure respectively, m the mass flow rate 
with dm = pudA 
Inspection of relations (1), (2), and (3) shows that, under the 
stipulated conditions, the decay of a nonuniformity is accom 
panied in the subsonic domain by an increase of Mach Number 
and a decrease of pressure and, in the supersonic domain by a 
decrease of Mach Number and an increase of pressure and tem 
perature. The static temperature in the subsonic domain may 


increase or decrease depending on the initial conditions. It 


can be seen that the parameters used vary in a similar manner 
as in the case of a Raleigh process with heat addition rhe 
entropy rise reaches a maximum at a Mach Number of one 
For example, if at a section 1 
V7, = 0.80 n, = 0.90 


then at a section 2, if a complete decay occurred (» 10) one 
I 


obtains 


VW = O.835 7./T, = 1.002 poi} 0.959 


An Electrochemical Method of Visualizing 
the Boundary Layer* 


Edward W. Geller 

Assistant Professor of Aeronautical Engineering, University of 
Alabama, University, Ala 

July 16, 1955 


A DIRECT METHOD of obtaining the velocity profile of the 
boundary layer is accomplished by marking a line element 
of fluid perpendicular to the surface over which the boundary 
layer is to be studied. As it moves downstream such a fluid 
element which has been made visible acquires the shape of the 
velocity profile, which can then be recorded directly on film by 
taking a photograph at the appropriate time 

The problem in applying this idea is to find a method to mark 
in element of the fluid. At the Institute for Technology at 
Stuttgart, Germany, Wortmann! accomplished this by using an 


d by Edward W. Geller as research for a 
lance of Dr. August Raspet 


* This technique was develope 
Master’s thesis and under the helpful 


Aerophysics Department, M State College 
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Fic. 1. Velocity profile is formed by the row of bubbles 
released at the wire. Wire diameter is 0.001 in rhe free 
stream velocity is approximate ly 1 ft. per sec 
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electrochemical reaction along a small wire placed in water flow. 
This work prompted the development of a similar technique at 
the Aerophysics Department, Mississippi State College. A 
short description of the method follows. 

The working fluid for this technique is necessarily water be 
cause of the chemical reaction used. A fine platinum wire is 
mounted perpendicular to the surface over which the velocity 
profile is to be obtained. The boundary layer along one wall of 
a water tunnel was used to demonstrate the technique. When 
a voltage pulse is applied between the wire and another electrode 
in contact with the water, a row of small bubbles is formed along 
the wire by the electrolysis which readily takes place if a small 
amount of salt has been added to the water. In this manner the 
necessary fluid element is marked. The row of bubbles is dis 
torted into the shape of the velocity profile as it is swept down 
stream. By taking a microphotograph at the appropriate time 


i velocity profile such as shown in Fig. 1 is recorded. 
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Comment on ‘‘Calculation of Span Loading 
for Arbitrary Plan Forms’’ 


Johannes Weissinger 

Ordentlicher Professor, Institut fur Angewandte Mathematik, 
Technische Hochschule, Karlsruhe, Germany 

June 28,1955 


RECENT NOTE of de Young! begins with the statement, ‘‘The 

Weissinger method (for example, see NACA Rep. 921) is 
limited to use on wings with a quarter chord line of constant 
sweep” which is correct only for my preliminary report.22 My 
full paper on the subject® contains on pp. 93-98 nearly the whole 
note of de Young (for instance integral equation, numerical 
method of solution, choice of m and A/) and also some numerical 


tables for the aid of the computer. 
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Aerodynamic Derivatives for Both Steady and 
Nonsteady Motion of Slender Bodies 


R. M. Wood and C. H. Murphy 
Ballistic Research Laboratories, Aberdeen Proving Ground, Md. 
July 1, 1955 


ive APPLICATION Of Munk’s slender body theory to the cal 
culation of aerodynamic coefficients describing steady motion 
(a, q constant) has been made by a number of authors.!~* Dor 
rance* has obtained slender body values of those coefficients 
associated with nonsteady angle of attack as a by-product of a 
linear analysis of the potential equation for nonsteady super- 
sonic flow. The stability derivatives for the most general motion 
of slender bodies with arbitrary cross section have been calcu 
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lated by Sacks.6 This treatment is, however, quite complex 


This note will present a simple derivation of all non-zero stability 
derivatives of a slender body of revolution 

We shall make the usual assumption that the cross flow, uw, at 
each cross section is a two-dimensional potential flow, and inde 
pendent of both the axial flow and the cross flows at neighboring 
cross sections. For a circular cross section, the potential of the 
flow with respect to stationary air is that of a doublet rhe 
kinetic energy per unit length can be found by the surface integral 
of a function of this potential,® and is of the form of u?/2 multi- 
plied by an ‘‘apparent mass,’’ which in this case is precisely the 


mass of the displaced fluid. Hence, 
momentum of the fluid per unit length = pAw (1 


where p = air density, 4 = cross-sectional area 

If we consider the momentum of the fluid at adjacent cross 
sections we see that it will depend on the local area and fluid 
velocity. Now the axial velocity, U, will move the fluid from 
one cross section to the next and hence from a region of higher 
momentum to one of lower momentum or vice versa. This local 
change in momentum gives rise to a differential force which can 
then be integrated to give the total force on the body Since we 
are then actually interested in the location of the cross section of 
the fluid and its motion, we will interpret x to be its location with 
respect to the base. (See Fig. 1.) 

We shall make use of a nonrolling right-handed coordinate 
system YYZ with X pointing forward along the missile’s axis 
and the Y axis initially in the horizontal plane, where a nonrolling 
coordinate system will mean one whose x-component of angular 
velocity is zero. (We make use of the tilde superscript to dis- 
tinguish these axes from the standard aerodynamic missile- 
fixed axes. In the case of a spin-stabilized body of revolution, 
those axes would be spinning rapidly with respect to our non- 
rolling axes and, hence, with respect to air at infinity 

As is shown in reference 7, the treatment of time derivatives 
of the dynamic variables—i.e., @, gG—is much more natural in 
nonrolling coordinates. Sacks® uses a missile-fixed system but 
specifically measures & and 6 in a nonrolling system. However, 
he does measure g and r in the rolling system and therefore ob 
tains a non-zero value of Ci, and C,,,,, for a body of revolution. 

If we consider projections of the velocity vector and a fixed- 
direction vector on the XZ plane, and with a and @ as defined in 


Fig. 1, then 


u = —Uat+ (x — Xeg)q (2 
where U = —dx/dt = axial component of the velocity of the 
missile, x., = the distance of the center of gravity from the 
base, x = the distance of fluid cross section from the base, g = 6 


By Newton’s Second Law and Eqs. (1) and (2) the force/unit 
length in the 2 direction, dZ/dx, is the time rate of change of 
momentum per unit length. Hence, dZ/dx = (d/dt) (pAu) and 
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Since we assume that the total velocity V is constant, its axial 
component is not necessarily constant lU’ is of second order 
however, ind is neglected 


Integrating over the length of the missile, one obtains 


Z po -aS; — (a/U)vd* — 
s230e q [2 \ \ d } 
whert 
v] 
dA 
>) = — dx base area 
J0 dx 
ey : 
vd = J 1 dx is dimensionless volume 
0 
nl | ; 
(a Y-g)vd? = { (x — 2 1 dx (x, is the x-coordinate of 
e 0 
the centroid 
d = diameter of the reference cross section 
If we multiply Eq. (3) by —(x — x.y) as is required by the indi 


cated positive direction of 7 (see Fig. 1), and integrate, we can 
obtain the moment about the center of gravity: 


VU »L? }al(vd® — xSv) + (a/l vax. — x _ 
(qg/U) [Xeg2So + (Xe — Xeg)vd*] — (G/U?)k*vd5} (5 
PZ 
where vk?d> = I, (x — Xeg)?A dx; (k is the dimensionless trans- 
e 


verse radius of gyration of the homogeneous model about its 
center of gravity ) 

An examination of Eqs. (4) and (5) shows that, under our as 
sumptions, Z and A depend linearly on a, &, g, and g. In order 
to obtain the relations for the various stability derivatives we 


first write the equations defining them. 


Z = (pV?S/2) [CZqa a + Cia (ad/2V) + Cizq (qd/2V 

C74 (gd*/4 V? (6 
M = (pV?Sd/2) [Cina a + Cina (ad/2V) + Ching (qd/2V 

Cing (Gd?/4 V? (7 
where § the reference area, V = the total velocity 


If Eqs. (4) and (5) are compared with Eqs. (6) and (7), the 
nondimensional { = x/d is introduced, and it is assumed that for 
small angles of attack and yaw U can be replaced by V, then the 


following relations are valid 


Cia 251,/S Cia = (—25/S)Xeg + 20/( S/d? 
Cie Wv/( S/d?) Cia = —40/(S/d?) (Xeg — X 
Cz¢ 45,/S)Xeq Cig = (—4S0/S)x 
$v /(.S/d? Veg — X s 
( Cing —Rrk2/(S/d 
_ S S/d? if 
Reg Be 


The corresponding relations for ) and V can be found from sym- 


metry considerations. It can easily be shown that 


Cys = CZa Cig = —Cihe 
Cyr = —Ciq Cie = Cyr | 
Cyp = Cia Cag = —Cra ? 
Cy, = —CZ; Cixy = Cirg 


Since the derivation of Eqs. (9) was based on the assumption 
of nonviscous flow, this theory makes no prediction of Magnus 
force and moment derivatives. We see, however, that Munk’s 
elementary theory based on the concept of apparent mass of each 
cross section will provide slender-body values of all 16 first-order 
non Magnus aerodynamic coefficients 

Reference 8 includes a conversion of the above results to ballis 
tic nomenclature for two types of motion and also gives an 
order-of-magnitude discussion showing that, for projectiles or 
missiles in air, the angular acceleration terms (arising from the 


4 coefficients) are negligibly small and are usually neglected 
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Comment on ‘‘Hypersonic Similarity and the 
Tangent-Cone Approximation for 
Unyawed Bodies of Revolution’’' 


V. Halbmillion 
Naval Ordnance Laboratory, White Oak, Silver Spring, Md 
June 29,1955 


I REFERENCE I, the authors have successfully attempted to 

combine the tangent-cone approximation with analytical 
solutions for the surface pressure on a slender cone in order to 
obtain an analytic representation of the pressure distribution on 
ogives. A justification for the use of the tangent-cone approx 
mation was made by referring to a report by Halbmillion and 
Kulishek,?, and the method was described as follows each 
transverse section of the body is replaced by a section of the 
coaxial cone tangent to the body surface the pressure 1s then 
assumed to be that which would exist on the cone at the same 
free-stream Mach Number, and the surface pressure on the 
cone is obtained by a numerical solution of the Taylor-Maccoll 
equations 

The implication that the method as described above was used 
and justified by Halbmillion and Kulishek is not quite correct 
The curves of reference 2 originally presented by Halbmillion 
and McGraw in reference 3 are based on the following ipproxl 
mate law 

At a given Mach Number, the pressure ratio on the surface ol 
an ogive at a given point is dependent only on the flow inclina 
tion—i.e., the tangent-cone angle—at that point To a given 
local flow inclination fer any ogive corresponds a single value of 
the pressure ratio (The pressure ratio is defined as the ratio of 
the local static pressure to the static pressure in the undisturbed 
stream 

Fowler and French‘ had made the additional assumption that 
this value of the pressure ratio is equal to the Taylor-Maccoll 
value for the tangent cone. But a check against experimental 
and analytical data which became available later indicated that 
this approximation is not fully satisfactory. Fortunately, the 
major discrepancies in pressure ratios occur at low angles of 
flow inclination and therefore do not affect substantially the 
wave drag calculations. Nevertheless, the writer felt that, for the 
purpose of predicting wave drag and pressure distributions, the 
Taylor-Maccoll values should be corrected For ex imple , at the 
junction between ogive and cylindrical body, where the tangent 
cone angle is zero, the Taylor Maccoll value of the pressure ratio 
is unity. Experimental data and computations by the method 
of characteristics indicate an average value of 0.83, and this was 
one of the key values selected in establishing the pressure ratio 
curves of reference 2 Other values of the pressure ratio for 
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various values of the Mach Number and the tangent-cone angle 
were similarly obtained statistically and the curves were inter- 
polated and extrapolated by assuming functional variations 
similar to those derived from the Taylor-Maccoll theory. 

Therefore, the curves presented in reference 2 should generally 
yield pressure distribution and drag values closer to the true 
values than those derived by Taylor-Maccoll theory without 
introducing experimental corrections. Since these curves ap 
peared, a large number of experimental checks were made at 
Mach Numbers ranging from 1.5 to 5.0, and the results were 
always very satisfactory for engineering purposes. In order to 
secure greater accuracy, it is necessary to resort to the second 
order theory of Van Dyke or to the method of characteristics 
used with a sufficiently fine mesh. The amount and complexity 
of the work involved are then very considerably greater than 
for the case of the tangent-cone method used with the curves 
of reference 2. 
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Variation of Maximum Thermal Stress in 
Free Plates 


P. J. Schneider 

Research Associate, Department of Mechanical Engineering, 
University of Minnesota 

April 11, 1955 


URRENT DEVELOPMENT of high-temperature structures em 
C ploying ductile and brittle materials has not only renewed 
past interest in the general theory of thermal stress and shock, 
but has stimulated, as well, a systematic study of several classes 
of structures with regard to susceptibility and full-range vari- 
ation of maximum thermal stresses for a variety of typical oper- 
ating conditions. 

A technically important structure within the class of one- 
dimensional systems is the plane plate whose thickness 26 is small 
relative to its face dimensions. In this case the plate temper- 
ature, ¢, varies only with «—the distance from and normal to the 
mid-plane of the plate. If, additionally, the plate is free of ex- 
ternal constraints, then the nonzero normal components of ther- 
mal stress, ¢, and o,, away from the plate edges, are identical 
and equal to! 
~] = *] 

KN )dX — (XN) +" X } «(X)XdX 

: 1 


J —1 ~ 


(l — pv) : 
: o(X) = 
J He} 2 
(1) 


Here v, E, and 8 are the Poisson’s ratio, modulus of elasticity, 
and coeflicient of linear thermal expansion for the plate material 
(each considered uniform), and XY = x/6 is the dimensionless 
distance from the plate mid-plane Y = 0. If the temperature 
distribution, ¢(x), or history, ¢(x, 6), is symmetric about Y = 0, 
then the third term in Eq. (1) vanishes, and the local steady or 
instantaneous thermal stress is simply proportional to the dif- 
ference between mean and local plate temperature, (7 — ¢). 

Six free plates are considered for a comparison of full-range 
Each plate is of 


resistance to maximum thermal stress alone 
semithickness 6, and, in addition to assuming uniform mechanical 
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Fic. 1. Maximum thermal stress in free plates 


properties v, E, and 8, the material for each is considered to be 
of uniform thermal conductivity k: 

Case I. Solid plate; uniform unit-volume heat generation, 
q’’’; faces X = +1lat ty 

Case II. Solid plate; positive linear heat generation, g,'’’|1 + 


a(t — t,)|,a >0;faces YY = +latt, 
Case III. Solid plate; negative linear heat generation, 
q’’ [1 + a(t —t,)],a <0; faces X = +l ath. 


Case IV. Porous plate; nongenerating; plate porosity P; 
fluid coolant of uniform thermal capacity C;, thermal conduc- 
tivity k;, and specific mass velocity G(— © < X¥ < — 1); coolant 
temperature ty,at Y = — ©; face Y = lath. 

Case V. Porous plate; uniform heat generation; porosity and 
coolant asin Case IV. 

Case VI. Solid plate; nuclear heating by nonuniform gamma- 
ray absorption, qi’’A exp[—A6(1 + Y)]; faces ¥ = +1 at ¢, 

Solutions for the dimensionless temperature profile, 7°,(.Y), 
and mean plate temperature, 7,, are given below,? along with 
the integrated dimensionless thermal-stress solutions, o,*(X), 


according to Eq. (1): 


Case I. 
T; = &(1 — X?) 2) 
7; = £1/3 (3) 
o,;* = (£;/2) [X? — (1/8)] $) 
Case IT. 
T11 = cos (vy ‘'t1,X) eos (4/ €:11) — 3 5 
Tir = (1 \/t11) tan (\/En1) — 1 6) 
o11" = (1/V/ én) tan (y E11) — cos (1/é11X) cos (YY é11 7 
Case ITI 
Tit = 1 — cosh (iV E11, X)/cosh (iv/ é1y1) 8 
Tin = 1 — (1/ivfEn) tanh (i+ é111) (9 
o11* = cosh (iy E111.V )/cosh iv én) — 
(l/ivféqir) tanh G@Y és) (10 
Case IV 
Tiv = exp [—éw(1 — X)] (11) 
liv = (1/2é1v) [1 — exp (—2é:v)] (12 
ow* = (3/2¢1v) } [1 + (1/é1v)] exp (—2é1v) + 
[1 — (1/éw)]}X — exp [—én(1 — X)] + 
(1/2é1v) [1 — exp (—2é1v)] (13) 
Case | 
Ty = (1/év?) $1 — exp [—év(1 — X)]} + (1/év) (1 + X) (14) 
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\ oO - r o r, 
6/k EBq\"'6/kR 
1 1 /2é\ 2é\ exp (—2¢ _ 1 /&y 15 
o l/év*) exp [—£év(1 \ — 
(3/2tv?) } [(1/és + 1/&, exp (—Zéy) — 16 
1/ty?) — 1/ty]j XY —(1/28 1 — exp (-2 
Case VI 
Ty1 1/Evi | — exp |—évill + \ — 
fl — exp —Véy1)] 2é\1 1+ \ iW 
Tv1 1/tvi) — (1/2év1 l/évi) + 1] X 
1 — exp (—2év; IS 
o\ 1 /tv1) exp [—évi(1 + X + 
(3/2tv1 (1 + év1) exp (—2év1) — 
1 — typ )|X — ffl — exp(—2Qévy)] /2Evi27f (19 


A full-range composite plot of the mean plate temperature and 
corresponding thermal stress for each case is presented in Fig. 1 
(see Table 1 

The well-known solution for Case J predicts a maximum tensile 
thermal stress at each face which varies directly with the heat 
generation parameter £; (@ arbitrary ) If, in Case I], the power 
production remains linear while the face temperatures are main 
tained at 4, then the maximum thermal stress is again tensile at 
the faces, and o1;* increases indefinitely as £1; approaches 7*/4. 
For the opposite situation (negative a) in Case ///, the maximum 
tensile stress at the plate faces approaches a limiting value of 
these cases, the thermal stress is 


on = | In all three of 


independent of the face temperatures ¢;.. This suggests that the 
same stress solutions apply for any set of surface temperatures, 
conductances and ambient tem- 


or any combination of surfacc 


peratures 
thermal 


Case I1V,t the maximum 


the downstream face Y = 1 


For the porous plate inl 
stress appears at The stress is 
compressive, and if the downstream face is maintained at fh, 
* approaches | for large values of the coolant parameter 


the dif- 


then o1y 
Er In this particular case, the stress is dependent o1 


ference in temperature between the downstream face of the plate 
and the upstream coolant If the porous plate generates uni- 


form heat, as in Case V, then the maximum thermal stress appears 


at the same face, but is tensile rather than compressive. The 
stress vanishes for large —y (Y = 1 not maintained at f¢,), and 
approaches a limiting value of oy* = 1/3 as éy tends to zero 


This latter limit exists because in this limiting condition the face 


temperatures increase indefinitely but approach an absolute 


difference of g’’’52/2k 

The exponential heat source in Case VJ produces a maximum 
tensile stress at the plate face Y = —1. Limiting stresses are 
zero, and oy;* reaches a peak of approximately 0.18 at a value 
of the absorption parameter éy; = 1.81. Again the maximum 
stress 1s independent of boundary conditions 

Manson 


ducing plates whos¢ 


has investigated thermal stresses in free power-pro 


thermal conductivity varies linearly with 


t Use of Eq. (1) for porous structures assumes that the stress behavior of 


porous materials is substantially the same as for solid materials Hence 
Without evidence to the contrary, the results of this extension must be con 


sidered to apply only to high-solidity plates 
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local plate temperature In this case the solution contains 


the ratio of local to face conductivity as a second parameter 
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Rapid Estimation of the Elastic-Plastic Euler 
Buckling Loads for Simply Supported 
Tapered Columns Under Varying Axial 
Loading 


Bertram Klein 

North American Aviation, Inc 
Calif 

August 3, 1955 


Aerophysics Department, Downey 


PROBLEM not given sufficient 


buckling of tapered col 


IMPORTANT STRUCTURAI 


A* ( 
attention 
umns attached to sheet with consequent variable axial loading 
the these buckling 


in the literature is the 


Various methods exist for calculation of 


loads. ! The approximate method presented here is direct, 
very rapid, yet sufficiently accurate and particularly suited for 
design. The column and its loading are shown in Fig. | The 


assumed shear flow need not be constant in magnitude 

The method concerns itself with the (?/E,/) diagram of the 
This quantity may be looked upon as representing the 
buckling intensity at E,l 
is relatively high, deflections tend to be higher, and consequently 
these values should carry more weight in determining buckling 


column 


any section In regions where (P? 


However, deflections tend to be higher around the middle of the 
column because of its end supports. Therefore it is suggested 


that the (P/E,/) diagram be weighted according to some func- 


tion which damps out at the end regions of the column. The 
buckling condition then becomes 
Lp I x? 
—~ f(x)dx f(x)dx = (1 
0 E,l an L? 


In general, the function, f(x), depends on the actual shape of the 


(P/E,1) diagram. However, it has been found that taking f(« 
sin (rx/L) gives good results in most practical cases 


Now let it be assumed that (P?/E;/) can be represented by the 


power series 


J > 
at OU + Cx* + GN rT CU T «) 


where uw is the nondimensional variable, x/L. Substitution of 


this series and the sin function into Eq. (1) and integrating yields 


a + 0.56 + 0.297c + 0.197d + 0.139e +4 r?/L? (3 
The left-hand side of Eq. (3) may be approximately but con 
veniently replaced by the expression 
1 ((P/Ed)yn0.2 + 2(P/EW), + (P/EW 

a+ 0.5b + O0.295« 0.197d + 0.1366 1) 


should count 


E,l 
Since in general the value (? 
should be 


chosen so that 


4 P 


its loading 


As mentioned previously, higher values of (7? 
more in determining buckling 


E,I + (P/E,]),~9.«, 


weighted more than the smaller one 


these two 


the larger one of 
Let u be 
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Fic. 1. Simply supported column and 
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(P/E WD) ymo.2 > (P/Ed)y 


dition takes the final form 


s( P ( P ( P x? 
1.25 + 2 + 0.75 = 
ca ; ei) .. oY ‘ Le 


When the stress at a section is plastic, it is necessary to use 


the corresponding tangent modulus in Eq. (4 
then becomes iterative in general 
steps, it may be necessary to revert back to Eq. (1) 

To illustrate the ease and rapidity of the method, several prob 
lems will be worked. Comparison with known correct solu 


tions will indicate the accuracy of the results. For the sake of 


simplicity all stresses will be assumed to remain elastic—i.e., Ey = 
E 
(1) P = Pol + a), 7 = ofl + u? 


1 Pof (1) (72 (+5) r? 
1.25 + 2 + 0.75 = (6) 
1 Ely 1.04 = 1.64 L? 


Eq. (5) becomes 


and Py = (8.46/L?)EIy (7) 
A solution obtained by another method agrees exactly with this 
result. 
(2) P = (1 — «)Po, 7 =T. Eq. (5) becomes 

1 Po O.8 0.5 0.2 r2 

. 1.25 + 2 + 0.75 = (8 

1 E] l I l L? 

and Po = (18.4/L2)EI (9) 


A more exact solution was found to be Py = (19.0/L2) EI, making 
the present answer in error on the safe side by about 3 per cent 
(3) P = Po, I = 1,(0.1 + 1.84) forO <u < 1/2,7 = Ky(1.9 — 
1.8u)for1/2 cu <1. 


1 Po (! °) (7) (722) r? 
4 } = (10 
1 FEJyjo| \0.46 1 0.46 L?2 


Po = 6.23( Eli2/L? (11) 


Substitution into Eq. (5) yields 


This result agrees within less than 4 per cent on the safe side of 


the exact solution.! 


(4)P=PiI = ui, 


Eq. (5) may be written 


a : ; 2 
— {(1.25) (5) + 2(2) + 0.75(1.25)] = (12) 
4 El, L? 
P = 3.53(E1,/L?) (13 


The error here is less than 1 per cent since the correct solution 


should be 


P = 3.50( El,/L?) (14) 
Sometimes the (P/EI/) diagram exhibits a maximum value 
around u = (1/2). In such cases f(x) = sin? 7 u gives good 


results. By reasoning similar to before, the following equation 


Gi) Jte 
= (15) 
Mies EL 


is derived 


7L(z)_..**(&), 


Then the proposed buckling con- 


The procedure 


If the column has one or more 


DECEMBER, 1955 
For example, let P/E] « sin ru. Then 
(P/ET)y2 (1/7) [0.588 + 5 + 0.588] = 2x?/L? lf 
and Pye = (11.2/L?) (El 17 


The exact solution is 
Pye = (11.4/L?) (El 


which makes the above approximate answer in error by less th 
2 per cent on the safe side 

It should be noticed that in the derivation of Eq. (5), th 
damping function does not correspond exactly to the sine func 
tion In fact, the function required to give Eq. (5) is close to 
(1.5 — uw) sin ru. This function will prove of use in problems i; 
volving columns with stepped taper. For example, consider th 
following case: P = P, I = (1/2)h for 0 < u < 0.5, an 
J =J/,for0.5 <u <1 
the above values of J into Eq. (1) to obtain 


20.5 ~ . 
} (1.5 — uw) sin 7 udu f 
t (1.5 — uw) sin w udu 
» ») 
I J0 1/2 0.5 - 


e ) 


El, *l : L 
(1.5 — uw) sin wr udu 18 
J ( 


) 


Substitute (1.5—w) sin ru for f(u) ane 


Evaluation of the integrals gives 


r 


P = 6.21E!,/L? 20) 
The correct solution is 
P = 6.40EI,/L? 


making the approximate answer in error by less than 3 per cent 


on the safe side 
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NUCLEAR DEVELOPMENT CORPORATION OF AMERICA 

PAN AMERICAN WORLD AIRWAYS, INC. . 

THE PARKER APPLIANCE COMPANY 

PESCO PRODUCTS DIVISION, BORG-WARNER COR- 
PORATION , 

PHILLIPS PETROLEUM COMPANY 

PIASECKI HELICOPTER CORPORATION 

THE PURE OIL COMPANY 

RADIOPLANE COMPANY 

J. B. REA COMPANY, INC. 

REACTION MOTORS, INC. 

REPUBLIC AVIATION CORPORATION 

E. V. ROBERTS AND ASSOCIATES 

ROHR AIRCRAFT CORPORATION 

PAUL ROSENBERG ASSOCIATES 

RYAN AERONAUTICAL COMPANY 

SANDBERG-SERRELL CORPORATION 

SCHRILLO AERO TOOL ENGINEERING COMPANY 

SHELL OIL COMPANY 

SIMMONDS AEROCESSORIES, INC. 

SOCONY MOBIL OIL COMPANY, INC. 

SOLAR AIRCRAFT COMPANY 

R. DIXON SPEAS 

SPERRY GYROSCOPE COMPANY DIVISION OF SPERRY 
RAND CORPORATION 

STANDARD OIL COMPANY OF CALIFORNIA 

STANDARD OIL COMPANY (INDIANA) 

STANLEY AVIATION CORPORATION 

THE STEEL PRODUCTS ENGINEERING COMPANY 

STROUKOFF AIRCRAFT CORPORATION 

STURGESS, INC. 

SUMMERS GYROSCOPE COMPANY 

THOMPSON PRODUCTS, INC. 

TOOLKO ENGINEERING COMPANY 

TRANS WORLD AIRLINES, INC. 

TURBO PRODUCTS, INC. 

UNION CARBIDE AND CARBON CORPORATIO 

UNITED AIR LINES, INC. 

UNITED AIRCRAFT CORPORATION 

UNITED STATES AVIATION UNDERWRITERS, INC 

UNITED STATES RUBBER COMPANY 

VARD, INC. 

VICKERS INCORPORATED 

WESTERN GEAR WORKS 

WESTINGHOUSE ELECTRIC CORPORATION 

WESTON ELECTRICAL INSTRUMENT CORPORATION 

WINDER AIRCRAFT CORPORATION 

WYMAN-GORDON COMPANY 

YOUNG RADIATOR COMPANY 














